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ABSTRACT 


Concave  curvature  has  a  relatively  large*  unpredictable*  effect  on  tur¬ 
bulent  boundary  layers.  Past  studies  of  turbulent  boundary  layers  on 
concave  walls  have  emphasized  quantitative*  single-point  measurements 
which  have  yielded  similar  information  about  the  general  features  of 
this  flow.  For  example*  skin  friction  and  turbulence  levels  increase 
when  the  boundary  layer  goes  from  a  flat  wall  to  a  concave  wall.  How¬ 
ever*  there  is  disagreement  over  the  cause  of  these  results.  Many 
recent  studies  have  supported  the  idea  that  a  large-scale  array  of 
counter-rotating  vortices  exists  within  the  concave  turbulent  boundary 
layer*  similar  to  the  Taylor-CSrtter  vortices  in  laminar  boundary  layers 
on  concave  walls.  Some  investigators  have  found  stationary*  large-scale 
spanwise  variations  in  mean  velocity  and  skin  friction.  However*  in 
other  studies  these  large-scale  variations  were  not  found.  Even  among 
those  who  did  find  stationary*  large-scale  variations*  there  are  some 
inconsistencies  in  the  relative  size  and  spacing  of  the  variations  com¬ 
pared  with  the  boundary  layer  thickness  and  the  radius  of  curvature*  the 
two  most  reasonable  length  scales.  Furthermore*  the  effects  of  concave 
curvature  on  the  basic  elements  (e.g.  streaks  and  bursts)  of  the  bound¬ 
ary  layer  structure  have  not  been  studied.  At  the  start  of  the  present 
study*  no  adequate  picture  of  the  overall  flow  field  existed. 

The  objective  of  the  present  study  was  to  obtain  an  adequate  qualita¬ 
tive  model  of  the  flow  in  order  to  increase  our  understanding  of  the 
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underlying  physics  and  thereby  guide  future  work.  A  large*  free-surface 
uater  channel  was  constructed  in  order  to  perform  visual  and  quantita¬ 
tive  studies  of  the  flou.  The  size  and  design  of  this  Hater  channel 
enabled  us  to  perform  detailed*  extensive*  visual  surveys  of  the  flou 
field  by  using  dye  injected  through  uall  slots  and  by  marking  the  flow 
with  hydrogen  bubbles.  Profiles  of  mean  velocity  and  turbulence  inten¬ 
sity  uere  also  obtained. 

Large-scale  structures  consisting  of  randomly  distributed  ejections 
and  sueeps  Here  found  to  dominate  the  boundary  layer  on  the  concave 
wall.  The  sueeps  appear  to  suppress  turbulence  production  near  the  uall 
by  inhibiting  the  normal  flat  surface  bursting  process.  The  ejections 
appear  to  enhance  turbulence  production  in  the  outer  portion  of  the 
boundary  layer  as  the  louer  speed  fluid  eonvected  from  regions  near  the 
uall  interacts  uith  the  higher  speed  fluid  farther  out. 

Interpretation  of  the  flou  visualization  Has  aided  by  a  computer 
simulation  of  a  simpler,  hypothetical  flou  field.  The  results  of  the 
simulation  helped  qualify  the  visual  observations  and  provided  a  basis 
on  uhich  to  build  a  more  realistic  picture  of  the  flou  field. 

The  streamuise  velocity  Has  measured  using  a  computer-controlled  hot 
film  anemometer.  Profiles  of  mean  velocity  and  turbulence  intensity  shorn 
the  flou  upstream  of  the  atart  of  curvature  to  be  a  normal  turbulent 
boundary  layer  uith  acceptable  spanwise  uniformity  and  a  momentum  thick¬ 
ness  Reynolds  number  of  1300.  Spantiise  surveys  in  the  curve  and  in  the 
recovery  region  dounstream  of  ourvature  shorn  no  evidenoe  that  the 
large-seale  structures  observed  visually  tend  to  occur  in  preferred 
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locations.  Boundary  layer  profiles  reveal  a  fuller  mean  velocity 
profile  and  an  increase  in  turbulence  intensity*  especially  in  the  outer 
flow.  Although  logarithmic  regions  are  evident  in  profiles  taken  in  the 
curved  section,  these  "quasi -logarithmic"  zones  appear  to  be  unrelated 
to  the  usual  Lau  of  the  Mall  and  migrate  away  from  the  wall  as  the  flou 
develops. 

Additional  insight  into  the  flou  physics  uas  provided  by  visualizing 
a  transitional  boundary  layer  along  the  concave  wall.  Taylor-65rtter 
vortices  were  observed  along  uith  a  secondary  instability  in  the  form  of 
a  horseshoe  vortex  street.  These  observations  demonstrated  a  possible 
mechanism  for  turbulence  production  in  the  outer  flou  for  the  case  of 
turbulent  boundary  layers. 

As  a  result  of  this  study,  a  qualitative  picture  of  the  flow  field  is 
proposed  uhich  accounts  for  the  observed  characteristics  of  turbulent 
boundary  layers  on  concave  ualls. 
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GLOSSARY 


As  the  study  of  the  structure  of  turbulent  flows  has  progressed *  Many 
words  and  phrases  have  evolved  to  describe  various  elements  or  charac¬ 
teristics  of  the  flow.  Obviously*  clear  definitions  are  needed  in  order 
to  communicate  the  particular  attributes  of  a  given  flow.  What  follows 
is  a  glossary  of  some  of  the  terms  used  in  the  present  work.  Many  of  the 
definitions  follow  those  of  Kline  C  Falco  (1979). 


breakdown 


bursting 


coherent  struc¬ 
ture  or  eddy 


ejection 
or  outflow 

large-scale 


longitudinal 

vortex 


An  abrupt  event  in  which  the  streak  oscillations 
terminate  in  the  formation  of  a  large  region  con¬ 
taining  a  wide  range  of  small  scales. 

The  set  of  processes  beginning  with  the  lifting  of 
streaks  and  terminating  at  the  end  of  the  mixing 
region. 

A  confined  region  in  space  and  time  in  which  defi¬ 
nite  phase  relationships  exist  among  flow  variables. 
A  coherent  structure  moving  away  from  a  wall. 

Having  a  size  of  the  order  of  the  boundary  layer 
thickness. 

A  region  of  fluid*  highly  extended  (aspect  ratio 
greater  than  10:1)  In  the  flow  direction*  with  rela¬ 
tively  high  streamwise  vorticity.  Typical  fluid  par¬ 
ticles  near  the  core  perform  several  revolutions 
during  thsir  stay  within  the  vortex. 


-  xiii  - 


K* 


medium-seal e 


Having  a  size  in  wall  layer  (y4)  coordinates  between 


t 


mixing  region 


roll  cell 


smal I -scale 


streak 


streak  lifting 


sueep 
or  inflow 


Taylor-Gortler 

vortices 


50  and  300. 

Region  after  a  breakdown  in  which  chaotic  motions 
affect  a  large  propagating  region. 

A  large-scale,  slowly  revolving  body  of  fluid  with¬ 
out  a  strong  vortex-like  core. 

Having  a  size  in  wall  layer  (y4)  coordinates  between 
1  and  10. 

A  low-speed  (relative  to  the  mean)  region  in  the 
viscous  sublayer.  highly  extended  (aspect  ratio 
greater  than  10:1)  in  the  flow  direction. 

Outward  movement  of  fluid  in  a  streak  to  a  point 
outside  viscous  sublayer. 

Large-scale  inward  (towards  the  wall)  motion  of 
faster  moving  fluid,  producing  local  acceleration  in 
the  flow  field. 

A  stationary  array  of  large-scale,  counterrotating, 
longitudinal  vortices  lying  next  to  a  wall  whose 
driving  mechanism  is  a  centrifugal  instability  due 
to  flow  curvature. 


viscous  sublayer  Region  where  y*  is  less  than  7-10. 
or  wall  layer 


Chapter  I 


INTRODUCTION 


1.1  msML  neiimi 

Surfaces  with  concave  curvature  occur  in  many  flow  situations  of  practi¬ 
cal  importance.  A  concave  surface  is  one  in  which  the  center  of  curva¬ 
ture  lies  outside  of  the  body,  in  the  flow,  as  shown  in  Fig.  1.  Some 
examples  are  shown  in  Fig.  2.  They  include  the  pressure  sides  of  tur¬ 
bine  blades  or  wings  and  portions  of  ducts  and  nozzles.  Often  the  bound¬ 
ary  layers  on  such  surfaces  are  turbulent.  Previous  studies,  including 
So  G  Mel  lor  (1975),  Meroney  C  Bradshaw  (1975),  and  Hoffmann  G  Bradshaw 
(19A1),  have  found  large  local  variations  In  boundary  layer  parameters; 
spanwise  variations  of  20X  in  skin  friction  and  100X  in  momentum  thick¬ 
ness,  even  for  cases  of  mild  curvature  (G/R  as  low  as  IX),  are  not 
uncommon.  The  more  efficient  design  of  concave  surfaces  has  been  ham¬ 
pered  by  our  inability  to  predict  the  large  effects  concave  curvature 
has  on  turbulent  boundary  layers.  The  lack  of  adequate  predictive 
methods  is  due,  in  part,  to  a  gap  in  our  understanding  of  the  physics  of 
the  flow. 

The  work  described  here  is  part  of  a  larger,  ongoing  program  at  Stan¬ 
ford  to  study  and  gain  a  better  understanding  of  the  effects  of  curva¬ 
ture  on  the  hydrodynamics  and  heat  transfer  in  turbulent  boundary  lay¬ 
ers.  The  study  of  turbulent  boundary  layers  on  convex  walls  was  well 
advanced  at  the  time  the  current  study  was  undertaken.  Simon  C  Moffat 


1 


Figure  2:  Examples  of  flows  with  concave  curvature. 


(1979)  and  01111s  C  Johnston  (1900)  had  taken  detailed  heat  transfer  and 
hydrodynamic  measurements  In  boundary  layers  on  eonvex  surfaces.  They 
had  also  successfully  modelled  the  effects  of  convex  wall  curvature  by 


suitably  Modifying  the  boundary  layer  code  STAN5.  Ho  previous  work  on 
concave  waits  had  been  done  at  Stanford.  Related  studies  had  been  done 
on  the  effects  of  system  rotation  on  turbulent  channel  flow  by  Johnston 
et  al  (1972). 

A  review  of  the  literature  revealed  a  lack  of  direct  visual  observa¬ 
tion  of  turbulent  flows  over  concave  walls.  Physical  models  of  the  flow 
structure  were  limited  to  sketchy  analogies  between  this  flow  and  the 
corresponding  laminar  flow  with  its  Taylor-Gortler  longitudinal  vor¬ 
tices.  The  effects  of  concave  curvature  on  the  streaks  and  bursts  found 
in  normal  flat  plate  boundary  layers  were  unknown.  Yet,  the  quantitative 
data  showed  a  need  for  an  overall  picture  of  the  flow  field  in  order  to 
clarify  the  interpretation  of  the  available  results.  Me  therefore 
decided  to  concentrate  on  visual  studies  in  order  to  develop  a  qualita¬ 
tive  model  of  the  flow  structure.  This  qualitative  model  can  then  be 
used  to  guide  later,  more  quantitative  studies  and  thus,  ultimately,  to 
aid  attempts  to  predict  the  flow.  These  steps:  (i)  a  qualitative  model 
from  visual  studies  as  a  basis  for  (ii)  the  design  of  quantitative 
experiments  for  use  in  (ill)  a  model  for  use  in  prediction  have  proven 
to  be  fruitful  in  the  study  of  many  complex  turbulent  flows. 

1-2  LITERATURE  SURVEY 

Until  1970,  relatively  little  research  had  been  done  on  the  effects  of 
concave  curvature  on  turbulent  boundary  layers.  However,  the  last  decade 
has  seen  a  large  Increase  in  work  in  this  field.  Table  1  presents  a  sum¬ 
mary  of  past  work,  with  an  emphasis  on  presenting  and  comparing  the 
structural  characteristics  of  the  flow  found  by  different  investigators. 


Past  studies  of  turbulent  flow  over  concave  surfaces. 
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Rayleigh  (1917)  first  deduced  that  flou  on  concave  surfaces,  in  this 
case  inviscid  flou  uith  circular  symmetry,  is  unstable  to  radial  dis¬ 


turbances.  He  assumed  conservation  of  angular  momentum  and  performed  an 
energy  calculation.  He  considered  tuo  lumps  of  fluid  uith  equal  masses, 
mi=m2::m  moving  in  circular  paths  uith  radii  n  and  rz.  The  energy  of  the 
tuo  masses  is 


E=l(miVt2+»2VjJ) 
=fm(vi2+vj2) . 


(1.1) 


He  then  interchanged  the  positions  of  the  tuo  lumps.  In  the  absence  of 
friction,  conservation  of  angular  momentum  dictates  that  vr=k,  a  con¬ 
stant.  for  each  lump  of  fluid.  The  neu  velocities  for  each  lump,  denoted 
by  O'  become 

Vi's—  (1.2) 

r  2 
*2 

V2/=— .  (1.5) 

r  1 

The  difference  in  energy  uorks  out  to  be 


E 


i]H-£  ■  ill) 


m  f  1  1 

:-(k22-kl2> - 

2  In*  r22 


(1.4) 


Let  us  arbitrarily  choose  n  to  be  less  than  rr.  If  ri<ra  then  the  sign 
of  E'-E  is  the  same  as  ka-ki.  Three  possible  eases  can  be  considered: 


1.  If  ki=k*  as  in  potential  flou,  there  is  no  change  In  the  energy. 
This  is  the  ease  of  neutral  stability. 
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2.  If  ki<kj  as  in  the  boundary  layer  on  a  convex  surface,  then  the 
energy  of  the  system  is  increased  by  the  interchange  of  the  tuo 
lumps  of  fluid.  Presumably,  a  system  tends  to  louer  its  energy 
state.  Hence,  if  ki<kj,  then  the  tuo  lumps  resist  being  inter¬ 
changed  and  the  system  is  stable. 

3.  If  ki>kj  as  in  the  boundary  layer  on  a  concave  surface,  then  sys¬ 
tem  can  louer  its  energy  by  interchanging  the  tuo  fluid  elements. 
The  system  is  therefore  unstable. 

Although  Rayleigh's  analysis  provides  a  rationale  for  anticipating  inst¬ 
abilities,  it  is  far  too  limited  to  be  useful  in  actual  flous,  uhere 
viscosity  plays  an  important  role. 

Early  experimental  studies  by  Milcken  (1930)  and  Mattendorf  (1935) 
obtained  only  mean  pressure  and  velocity  measurements  in  curved  chan¬ 
nels.  From  these  measurements  mixing  lengths  and  eddy  viscosities  uere 
calculated.  Milcken  investigated  turbulent  boundary  layers  developing  on 
concave  ualls  and  found  that  the  mixing  length  he  calculated  uas  about  4 
times  that  for  a  comparable  flat  plate.  Mattendorf,  using  fully  devel¬ 
oped  channel  flou,  also  found  an  increase  in  mixing  on  the  concave  uall 
as  uel!  as  a  midchannel  profile  uhich  reached  a  constant  angular  momen¬ 
tum  distribution,  another  sign  of  vigourous  mixing.  Higher  skin  friction 
is  also  evident  from  an  examination  of  his  velocity  profiles  close  to 
the  concave  ual . . 

After  1935,  published  experimental  uork  on  curved  turbulent  flous 
ceased  for  about  tuo  decades.  In  the  meantime,  oBrtler  (1940,  1954) 
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performed  a  viscous  stability  analysis  of  laminar  boundary  layers  on 
concave  sails.  He  superimposed  small  disturbances  in  the  form  of 
counter-rotating  longitudinal  vortices  (as  shown  in  rig.  3)  on  the  mean 
flow,  linearized  the  equations,  and  solved  the  resulting  eigenvalue 
problem.  He  found  that  vortices  with  a  wide  range  of  wavelengths  could 
be  amplified  and  that  the  parameter 


on  an 


adaptation  of  the  problem  of  flow  between  concentric  cylinders  first 
solved  by  C.  1.  Taylor.  It  was  another  10  years  before  Taylor-Gortler 
vortices  were  actually  observed  by  Gregory  £  Walker  (1950). 

Eskinazi  £  Yeh  (1956)  resumed  the  experimental  investigation  of  fully 
developed  channel  flow.  Using  hot  wire  anemometers)  they  measured  mean 
and  turbulent  quantities  as  well  as  spectra.  They  obtained  good  general 
agreement  with  Wattendorf's  results.  Their  data  and  considerations  of 
the  turbulent  energy  equations  showed  that  turbulent  fluctuations  are 
amplified  on  the  concave  wall.  As  with  previous  experimental  research¬ 
ers#  no  large-scale  disturbances  were  expected  or  looked  for. 

Tani  (1962)  is  generally  credited  with  being  the  first  to  find  evi¬ 
dence  of  large-scale  structures  in  a  turbulent  boundary  layer  on  a  con¬ 
cave  wall.  Stationary,  spanuise  periodic  variations  in  velocity  were 
detected  near  the  surface  of  a  test  wall  with  mild  curvature  inserted  in 
a  wind  tunnel.  Unfortunately  the  investigation  was  very  limited,  and 
little  was  done  to  elaborate  on  the  nature  of  the  structures.  An  attempt 
was  made  by  Tani  to  apply  the  Gortler  stability  parameter  (see  above)  to 
the  turbulent  flow  by  using  an  eddy  viscosity  rather  than  the  kinematic 
viscosity. 

Patel  (1968)  also  made  limited  measurements  on  a  concave  wall.  He 
detected  large-scale,  stationary  variations  in  the  pressure  as  he  tra¬ 
versed  the  span  with  a  Preston  tube.  Variations  as  high  as  50X 
(peak-to-peak)  were  reported  with  a  spanuise  wavelength  about  twice  the 
initial  boundary  layer  thickness.  In  addition,  the  same  spanuise  profile 
was  present  at  two  streamuise  locations,  135°  and  173°  into  the  curve. 
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Studies  by  Johnston  et  al  (1972)  investigated  the  qualitatively 
similar  situation  of  flou  in  a  rotating  channel.  A  rectangular  channel 
uith  nominally  tuo-dimensional »  fully-developed  duct  flou  uas  rotated 
about  an  axis  outside  the  channel  and  perpendicular  to  the  flou.  Flou 
visualization  using  dyes  and  hydrogen  bubbles  revealed  large-scale 
motions  on  the  destabilized  side  which  penetrated  far  into  the  core.  The 
positions  of  the  "roll  cells",  as  they  uere  called,  uere  not  fixed,  but 
varied  uith  time.  The  "short  time"  average  roll  cell  uavelength  uas  com¬ 
parable  to  the  channel  uidth. 

Bradshaw  (1973)  did  a  comprehensive  review  of  all  previous  work  per¬ 
taining  to  the  effects  of  streamline  curvature  on  turbulf  >ce.  Amon»  the 
many  topics  covered  uas  the  analogy  between  certain  configurations  of 
laminar  curved  and  buoyant  flows.  For  turbulent  flows,  Bradshaw  pre¬ 
sented  the  following  argument:  A  fluid  Clemen*  in  a  buoyant  flou  uith  a 
density  p'  higher  than  the  mean  density  p  experiences  a  vertical  accel¬ 
eration  of  -p'g/p.  A  fluid  element  in  a  curved  flow  with  a  velocity  u 
higher  than  the  mean  velocity  U  experiences  an  outward  acceleration  (in 
excess  of  the  acceleration  attributable  to  curvature  of  the  mean  flou) 
of 

[ (U+u) *-Ul 1  2Uu 


The  connection  between  the  two  situations  comes  when  one  notes  that  in  a 
variable  density  shear  layer  the  correlation  coefficient  between  densit" 
and  velocity  variations  has  been  experimentally  found  to  be  as  high  as 
0.9  (0.  S.  Johnston  (1959)1.  Hence,  in  both  cases  the  fluid  element 
experiences  an  excess  body  force  proportional  to  the  fluctuating 
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property*  and  the  correlation  data  insures  that  the  results  will  be 
comparable. 

In  the  mid-seventies*  there  was  a  dramatic  increase  in  the  amount  of 
research  being  done  on  concave  flows.  Ellis  £  Joubert  (1974)  mentioned 
the  persistence  of  some  spanwise  velocity  variations  in  their  study  of 
developing  turbulent  boundary  layers.  The  variations  could  not  be  elimi¬ 
nated  despite  "gross  duct  distortions  and  deliberate  attempts  to  induce 
asymmetric  flow."  However,  detailed  investigation  of  the  variations  was 
not  undertaken. 

So  £  Mel  lor  (1975)  performed  a  careful  study  in  which  they  controlled 
secondary  flows  due  to  end  wall  boundary  layers  by  using  wall  jets.  They 
also  tailored  the  shape  of  the  convex  wall  to  minimize  the  streamwise 
pressure  gradient  on  the  concave  wall.  A  stationary  pattern  of  velocity 
variations  was  found;  velocities  at  30°  of  turning  varied  by  up  to  35% 
(peak-to-peak) .  However,  there  was  no  well-defined  periodicity  in  the 
variations.  By  the  100°  mark  the  spanwise  variations  in  mean  velocity 
had  settled  down  to  8%.  Boundary  layer  thickness  was  reported  to  vary  by 
a  factor  of  2  over  a  spanwise  distance  of  half  the  minimum  local  bound¬ 
ary  layer  thickness  at  both  the  30°  and  1 00°  stations.  Traverses  of 
turbulence  quantities  were  made  at  the  spanwise  locations  of  mean  veloc¬ 
ity  maxima  and  minima.  Measurements  of  streamwise  turbulence  intensity 
showed  a  peak  at  y/$=.4  at  30°  for  all  spanwise  locations.  A  smaller 
peak  at  y/S=.8,  was  also  detected  but  only  at  the  mean  velocity  minima 
locations  at  110°.  Other  turbulence  quantities  also  Indicated  increased 
turbulent  activity  at  positions  away  from  the  wall.  Skin  friction  was 
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inferred  from  Clauser  plots.  The  authors  expressed  some  reservations 
regarding  the  validity  of  the  technique  in  the  regions  of  curved  flow 
since  the  usual  logarithmic  regions  in  the  mean  velocity  profiles  uere 
almost  nonexistent  in  the  curved  wall  profiles.  Cf  was  found  to 
increase  with  streamuise  distance  along  the  concave  uall  with  no  signif¬ 
icant  spanuise  variation.  A  value  of  Cf=0.0054  uas  reported  at  110° 
(compared  uith  Cf=0. 00351  upstream  of  the  curve). 

At  the  time  So  £  del  lor  uere  doing  their  studies.  Meroney  £  Bradshau 
(1975)  performed  a  similar  study  using  a  constant  uidth  channel  uith 
mild  (S/R=. 008)  curvature.  They  too  found  irregular  but  repeatable 
large-scale  variations  in  the  spanuise  velocity  profile  and  355C  varia¬ 
tions  in  boundary  layer  thickness  across  the  uavelength  of  a  typical 
"roll".  They  estimated  that  the  average  uavelength  of  the  rolls  uas 
close  to  the  local  average  boundary  layer  thickness.  The  stability  and 
irregularity  of  the  roll  positions  uas  attributed  to  upstream  flou  con¬ 
ditions.  Turbulence  intensity  measurements  again  shoued  increased  activ¬ 
ity  auay  from  the  uall  although  the  peaks  found  by  So  £  del  lor  uere 
absent. 

Ramaprian  £  Shivaprasad  (1977)  also  studied  the  case  of  mild  curva¬ 
ture.  but  they  found  no  spanuise  variations  in  flou  properties.  The 
authors  admitted  the  lou  aspect  ratio  of  their  apparatus  (2.5)  may  have 
inhibited  the  grouth  of  large-scale  structures. 

Fully-developed  curved  duct  flou  uas  investigated  by  Hunt  £  Joubert 
(1979).  The  radius  of  curvature  uas  relatively  large  in  order  to  obtain 
a  "perturbed  straight  flou".  Nevertheless.  careful  measurement  of  the 


velocities  in  the  core  revealed  a  stationary  structure  which  they 
attributed  to  longitudinal  vortices  (the  velocity  variations  were  less 
than  5X).  The  wavelength  of  the  main  structure  uas  approximately  1.614 
(M=channet  width),  but  measurements  of  friction  velocity  indicated 
weaker  structu.  es  in  between  the  main  ones,  so  the  wavelength  could  have 
been  closer  to  0.8W. 

Smits.  Young.  £  Bradshaw  (1979)  studied  the  effects  of  short  regions 
of  high  curvature  where  the  flow  was  impulsively  turned  20°  and  30°. 
They  found  a  pattern  of  spanwise  variations  in  skin  friction  and  other 
properties  which  persisted  virtually  unchanged  almost  60  boundary  layer 
thicknesses  downstream  of  the  end  curvature.  Peak-to-peak  variations  in 
velocity  and  skin  friction  were  10  to  20X  of  the  mean;  the  wavelength 
was  about  twice  the  local  boundary  layer  thickness.  They  found  that  the 
positions  of  the  spanwise  variations  appeared  to  depend  on  the  wind  tun¬ 
nel  screens.'  All  of  their  results,  including  measurements  of  Reynolds 
stresses  and  triple  products,  seemed  consistent  with  the  presence  of 
longitudinal  vortices. 

Prabhu  £  Rao  (1980)  studied  the  effects  of  strong,  extended  curvature 
for  several  different  inlet  boundary  layer  thicknesses.  Particular 
attention  was  paid  to  eliminating  secondary  flows  due  to  the  endwalls 
through  the  use  of  fences.  Stationary  spanwise  variations  as  high  as  30X 


'  Other  Investigators.  Including  Tan-atichat  (1982),  and  Jayaraman  £ 
Reynolds  in  current  work  at  Stanford  have  found  that,  in  general,  sig¬ 
nificant  nonuniformities  ean  be  introduced  by  mesh  screens,  depending 
on  how  they  are  selected  and  mounted  in  wind  and  water  tunnels.  Design 
criteria  for  the  screens  used  in  the  present  study  may  be  found  in 
Appendix  A. 
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(peak-to-peak )  of  the  potential  wall  velocity2  were  measured.  However, 
the  pattern  of  variations  did  not  change  appreciably  even  uhen  the  inlet 
boundary  layer  thickness  was  varied  by  a  factor  of  two.  The  pattern  also 
persisted  unchanged  as  the  boundary  layer  developed  downstream. 

Shizawa  et  al  (1981)  performed  the  first  (to  our  knowledge)  flow 
visualization  study  of  turbulent  boundary  layers  on  concave  walls  using 
the  smoke  wire  technique.  Large-scale  structures  were  visible  as  agglom¬ 
erations  of  smoke,  and  the  wavelengths,  for  two  different  curvatures, 
were  about  20  times  the  local  momentum  thickness.  Limitations  of  the 
visualization  technique  prevented  the  investigators  from  obtaining  more 
detailed  information  about  the  structf'e  of  the  disturbances.  However, 
momentum  thickness  was  observed  to  vary  by  a  factor  of  two  across  the 
span,  and  velocity  vectors  were  inclined  up  to  4°  away  from  the  nominal 
streamwise  direction  for  flow  near  the  concave  wall. 

In  an  attempt  to  impose  some  periodicity  on  the  spanwise  variations 
for  a  more  systematic  study  of  their  effects.  Hoffmann  £  Bradshaw  (1981) 
placed  vortex  generators  on  the  walls  of  the  settling  chamber  of  their 
wind  tunnel  (the  same  wind  tunnel  used  by  Meroney  C  Bradshaw).  Before 
installing  the  generators,  the  spanwise  skin  friction  profile  in  the 
curve  was  found  to  vary  erratically  with  variations  in  amplitude  up  to 
20X.  The  skin  friction  profile  was  virtually  identical  to  the  one  mea¬ 
sured  by  fleroney.  The  generators  were  then  installed  and  different  spac- 
ings  inveitigated.  A  spacing  of  1.48  resulted  in  the  maximum  downstream 


2  The  potential  wall  velocity  is  the  the  calculated  velocity  of  a  fluid 
particle  at  the  wall,  having  the  same  angular  momentum  as  fluid  in  the 
potential  core. 
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amplification  of  the  disturbances.  The  spanuise  profiles  obtained  under 
this  condition  of  maximum  amplification  took  on  a  well-defined  periodic¬ 
ity.  The  upstream  variations  in  skin  friction  increased  from  about  5X, 
for  the  no-generator  condition,  to  10X>  with  the  generators  set  for  max¬ 
imum  amplification.  Similarly,  skin  friction  variations  in  the  curve 
increased  from  10X  to  30X.  A  comprehensive  set  of  measurements  including 
triple  products  uas  then  taken  uith  the  generators  set  for  maximum 
amplification.  Final  discussion  of  these  results  by  the  authors  is  pend¬ 
ing  (private  communication,  P.  Bradshaw,  1981). 

Later  studies  by  Bippes  £  Gortler  (1972),  Aihara  £  Sonoda  (1981),  and 
Aihara  £  Koyama  (1981)  on  laminar  boundary  layers  concentrated  on  the 
development  of  the  Taylor-Gortler  vortices  along  concave  walls.  These 
references  also  contain  other  information  and  bibliographies  on  laminar 
flows.  Of  relevance  to  the  present  study  are  observations,  especially 
detailed  by  Aihara  £  Koyama  (1981),  that  the  upwash  region  between  pairs 
of  vortices  develops  a  secondary  instability  in  the  form  of  a  horseshoe 
vortex  street.  Blackwelder  (1981)  has  drawn  analogies  between  this  phe¬ 
nomenon  and  the  bursting  process  in  normal  turbulent  boundary  layers. 

A  full  review  of  past  work  on  the  structure  of  turbulent  boundary 
layers  on  flat  walls  will  not  be  presented  here.  Cantwell  (1981)  gives  a 
good  summary  of  the  current  state  of  knowledge  of  turbulent  flows,  and 
Savill  (1979)  highlights  the  evolution  of  the  physical  model  for  plane 
turbulent  boundary  layers.  Although  many  of  the  questions  regarding  the 
flow  physics  remain  unanswered,  certain  structural  features  are  gener¬ 
ally  agreed  upon.  A  concise  summary  taken  from  Cantwell  (1981)  follows: 
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"There  appear  to  be  four  main  constituents  of  the  organized 
structure.  Nearest  the  wall  is  a  fluctuating  array  of  stream- 
uise  counter-rotating  vortices...  The  vortices  densely  cover 
all  parts  of  the  smooth  wall.  Slightly  above  the  streamuise 
vortices  but  still  quite  close  to  the  wall  is  a  layer  that  is 
regularly  battered  by  bursts  that  involve  very  intense  small- 
scale  motions  of  energetic  fluid.  The  outer  layer  is  also 
occupied  by  intense  small-scale  motions.  These  are  found  pri¬ 
marily  on  the  upstream-facing  portions  of  the  turbul ent-non- 
turbulent  interface;  the  backs  of  the  bulges  in  the  outer  part 
of  the  layer.  The  outer  small-scale  motions  are  part  of  an 
overall  transverse  rotation  uith  a  scale  comparable  to  the 
thickness  of  the  layer." 

A  variety  of  models  have  been  proposed  to  explain  these  features  and  hou 
they  come  about;  the  field  is  still  a  subject  of  current  research.  Per¬ 
tinent  results  Mill  be  recalled  from  the  literature  as  needed  in  the 
discussion  of  results  of  our  study. 


Several  conclusions  can  be  draun  from  the  past  uork.  First,  recent 
investigations  indicate  a  strong  tendency  for  large-scale  transverse 
motions  to  develop  in  a  turbulent  boundary  layer  on  a  concave  Mall,  even 
for  cases  of  ueak  curvature.  The  precise  nature  of  these  disturbances 
has  not  been  made  clear  by  prior  research;  houever,  many  experimenters 
tend  to  attribute  the  disturbances  to  a  kind  of  turbulent  Taylor-Gortl er 
vortex  system.  The  uavelength  of  the  disturbances  tends  to  be  on  the 
order  of  tuice  the  local  boundary  layer  thickness.  In  no  case  do  the 
naturally  occurring  disturbances  display  a  well-defined  periodicity. 
The  pattern  of  observed  disturbances  seems  to  depend  on  the  specific 
characteristics  of  the  apparatus  employed  since  no  simple  scaling  laws 
are  evident  which  relate  the  flow  properties  to  the  disturbances.  The 
failure  of  the  earlier  experiments  to  detect  large-scale  variations, 
despite  evidence  of  strong  mixing  in  mean  velocity  profiles,  and  the 
several  recent  results  demonstrating  the  importance  of  inlet  conditions. 
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raise  further  questions  concerning  the  true  nature  of  the  flow  field. 
Despite  extensive  single  point  measurements  in  this  class  of  flows*  lit¬ 
tle  has  been  done  to  survey  the  instantaneous  flow  field  on  a  global 
scale*  for  example  by  using  flow  visualization*  to  clarify  the  physics 
of  the  flow. 

1.3  SPECK  TIC  QBJEC-T-I VES 

In  order  to  answer  some  of  the  questions  raised  by  the  review  of  the 
literature  above*  the  following  objectives  for  the  study  were  formu¬ 
lated: 

1.  Design  and  construct  a  large  water  channel  facility  whose  primary 
purpose  is  to  allow  detailed  visual  study  of  turbulent  boundary 
layers  on  concave  walls  with  moderate  curvature  (6/R=0.05>. 

2.  Obtain  mean  velocity  data  to  qualify  the  apparatus  and  quantify 
some  of  the  more  commonly  used  boundary  layer  parameters  such  as 
momentum  thickness  and  shape  factor.  Obtain  turbulence  data  as 
equipment  permits. 

3.  Conduct  detailed  visual  surveys  of  the  flow  field  upstream  of 
curvature  in  the  curved  region*  and  in  the  recovery  region* 
employing  dye  injection  through  the  wall  to  visualize  the  flow 
close  to  the  wall  and  the  hydrogen  bubble  technique  for  flow  far¬ 
ther  out.  Determine  the  effects  of  concave  curvature  on  flow 
structures. 
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4.  Develop  a  qualitative  physical  model  of  the  flow  which  explains 


its  observed  characteristics. 
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Chapter  II 


EXPERIMENTAL  METHODS 

2.1  FACILITY  DESCRIPTION 

In  order  to  perform  the  investigation,  a  completely  new  flow  facility 
was  designed  and  constructed. 3  A  photo  of  the  finished  apparatus  is 
shoun  in  Fig.  4  and  a  schematic  of  the  curved  channel  in  Fig.  5. 
Because  ue  designed  the  uater  channel  with  a  great  deal  of  flexibility 
in  mind,  detailed  design  calculations  and  parts  specifications  are 
described  in  the  Appendices  for  the  benefit  of  future  investigators  who 
may  uish  to  modify  the  rig. 

2.1.1  General  Construction 

The  facility  is  a  large  free  surface  uater  channel  uith  plexiglas  ualls 
braced  by  steel  frames.  Mater  enters  and  leaves  through  steel  inlet  and 
outlet  plenums.  The  actual  test  walls  are  inserted  into  the  main  struc¬ 
ture  so  that  uater  stands  on  both  sides  of  the  test  walls.  This  design 
frees  the  test  ualls  of  the  loads  imposed  by  the  large  volumes  of  uater 
and  enables  the  ualls  to  retain  their  shape  once  they  are  placed  and 
adjusted  in  position. 


3  Design  and  construction  of  the  facility  was  initially  divided  among 
three  students:  I  designed  the  curvature  uater  channel  and  probe  tra¬ 
verse.  John  Simonich  took  care  of  sump  modifications,  the  uater  treat¬ 
ment  system,  and  the  mobile  Instrument  platform,  and  Brian  Gyles  han¬ 
dled  the  straight  uater  channel  and  the  plumbing.  The  reader  should 
consult  Simonich's  thesis  for  additional  design  details. 
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Y  LAYER  TRIP8 


Mgure  5:  Schematic  of  concave  flow  facility. 


2.1.2  mist 


Water  is  pumped  from  an  underground  sump  by  a  7.5  hp  deep-well  impel  lor 
pump.  It  enters  the  rig  through  a  vertical  perforated  pipe.  Two  disc 
valves  and  a  bypass  leg  are  used  to  control  the  flow  rate.  A  constant 
height  skimmer  overflow  in  the  inlet  plenum  removes  much  of  the  surging 
due  to  the  pump  and  makes  it  easier  to  maintain  the  water  at  a  fixed 
height.  The  flow  is  then  passed  through  a  perforated  plate,  diffused, 
and  turned.  The  pressure  drop  across  the  perforated  plate  is  of  the 
order  of  100  velocity  heads  in  order  to  remove  the  larger  eddies.  A 
honeycomb  panel  and  four  fine  mesh  screens  are  used  to  take  out  the 
smaller  eddies  and  reduce  the  turbulence  levels  in  the  flow.  A  two- 
dimensional  <1:1  contraction  ratio  noz2le  accelerates  the  flow  up  to  the 
nominal  speed  of  15  cm/s. 

In  the  present  study,  laminar  boundary  layers  leaving  the  nozzle  were 
tripped  110  cm  downstream  of  the  nozzle  exit  by  0.47  cm  (3/16")  square 
rods.  The  wall  opposing  the  test  surface  was  designed  to  minimize  the 
streamwise  pressure  gradient  which  would  otherwise  develop  due  to  the 
blockage  effect  of  the  boundary  layers. 

2.1.3  Test  Section 

488  cm  (16  ft)  downstream  of  the  nozzle,  the  flow  is  turned  in  a  90°, 
136  cm  radius  bend.  The  opposite  wall  was  designed  to  minimize  the 
streamwise  pressure  gradient  on  the  curved  test  surface  due  to  pressure 
gradients  owing  to  streamline  curvature  and  boundary  layer  blockage, 
fol lowing  the  curve  is  a  244  cm  (8  ft)  long  recovery  region. 


2.1.4  Exit 


At  the  end  of  the  recovery  section  snother  perforated  plate  is  used  to 
keep  the  effects  of  the  nonuniform  flow  in  the  outlet  plenum  from  propa¬ 
gating  upstream.  It  uas  not  necessary  for  the  present  study  to  use  a 
honeycomb  flow  straightener  prior  to  the  perforated  plate  at  the  channel 
exit.  A  specially  designed  valve  in  the  outlet  plenum  controls  the  flow 
back  into  the  sump  uhile  at  the  same  time  preventing  the  entrainment  of 
air  into  the  water  as  it  returns  to  the  sump. 

2.1.5  Filtration  and  Air  Control 

The  water  uas  continuously  filtered  by  a  diatomaceous  earth  swimming 
pool  filter,  operating  in  an  independent  side  loop  (see  Fig.  6).  A  cavi- 
tating  nozzle  and  a  vacuum  pump  deaerated  the  water,  eliminating  the 
fine  "mist"  of  bubbles  usually  present.  These  bubbles  would  have  other¬ 
wise  obscured  the  hydrogen  bubbles  used  for  flow  visualization  and  also 
prevented  the  use  of  hot  film  anemometers.  The  water  was  chlorinated  to 
a  level  about  twice  that  recommended  for  swimming  pools  to  prevent  algae 
growth.  The  pH  of  the  water  uas  not  controlled. 


2.1.6  Instrumentation 

A  mobile  instrument  platform  was  designed'  and  constructed  to  run  along 
rails  attached  to  the  tops  of  the  channel  frames.  This  platform  carries 
the  probe  traverse  in  the  streamuise  direction.  The  speed  of  the  plat¬ 
form  is  monitored  by  a  photodiode  which  senses  the  rulings  on  a  tape 
scale  glued  to  one  of  the  rails. 


'  For  design  details  see  thesis  of  J.  Simonioh. 
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A  two-axis  probe  traverse  was  developed  to  position  probes  anywhere 


in  a  crossplane  to  the  flow.  Stepper  motors  drive  the  two  axes  and  are 
remotely  controlled  by  an  indexer.  Provisions  in  the  construction  of  the 
indexer  allow  computer  control.  An  interface  for  a  DEC  MINC-11  micro¬ 
computer  was  designed  and  built  for  the  indexer.  Details  concerning  the 
hardware,  electronics,  and  software  of  the  traverse  may  be  found  in 
Appendix  B.  A  picture  of  the  instrument  platform  with  the  traverse  is 
shown  in  Fig.  7. 
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2.2  FLOW 


TECHNIQUE 


Two  techniques  were  used  to  visualize  the  flow:  dye  injection  and  hydro¬ 
gen  bubble  generation. 


In  the  design  of  the  test  wall,  ue  included  several  pairs  of  dye 
slots  in  order  to  Inject  dyes  into  the  viscous  sublayer  (see  Fig.  8). 
Dyes  injected  into  the  sublayer  through  wall  slots  to  visualize  the 
boundary  layer  structure  were  first  employed  by  Rundstadler,  Kline,  £ 
Reynolds  (1963)  in  their  pioneering  work  on  turbulent  boundary  layer 
structure.  In  the  present  study,  each  dye  slot  uas  30  cm  longs  the  width 
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of  the  slot  mss  adjustable  and  typically  set  to  about  0.02  cm.  The  slots 
were  arranged  in  pairs.  with  about  20  cri  between  each  pair,  in  an 
attempt  to  track  flow  structures  over  a  greater  streamuise  extent.  Each 
slot  in  a  pair  uas  fed  a  different  color  of  dye  from  a  p>astic  bottle 
(ordinarily  used  for  medical  procedures)  suspended  about  half  a  meter 
above  the  water's  surface.  Small  aquarium  valves  controlled  the  feed 
rate  of  dye  to  the  slots.  The  dye  consisted  of  food  coloring  diluted 
with  water  along  with  a  small  amount  of  isopropyl  alcohol  to  make  it 
neutrally  buoyant  and  algae  resistant.  Tests  were  run  with  the  slots 
turned  on  and  then  off.  These  tests  indicated  that  for  modest  feed  rates 
(0.15  ml/s,  0.3  cm/s  through  the  slot),  the  injection  of  dye  liad  no 
qualitative  effect  on  the  flow.  Oyes  contaminating  the  water  supply  were 
gradually  bleached  out  by  the  chlorine  added  to  the  water  for  algae  pre¬ 
vention. 

Me  also  used  a  modified  syringe  with  a  20  gauge  hypodermic  spinal 
needle  to  inject  dyes  at  other  locations  for  quick  checks  of  the  flow 
field. 

The  dyes  were  illuminated  from  the  back,  through  the  transparent 
Plexiglas  test  surface,  by  banks  of  high  intensity  fluorescent  tubes.5 
High-speed  (ASA  400)  black  and  white,  and  color  reversal  films  (Kodak  4X 
and  VNF  7250)  were  used  to  photograph  the  flow. 


5  The  60  cycle  fluctuation  of  the  light  intensity  made  it  necessary  to 
adjust  the  shutter  speeds  on  our  cameras  to  1/60  s  or  slower  to  avoid 
shadows  and  inconsistent  exposures. 
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To  visualize  the  flou  far  from  the  uall,  ue  used  the  hydrogen  bubble 
technique.  This  method  of  visualization  was  invented  by  Clutter  £  Smith 
(1961)  and  subsequently  used  and  refined  by  several  Stanford  investiga¬ 
tors  including  Schraub  et  al  (1965).  A  platinum  wire  electrode  76 
microns  (3  mils)  in  diameter  and  30  cm  long  uas  mounted  on  a  suitable 
holder  and  placed  in  the  flou.  When  energized  to  a  potential  of  about 
125  volts  relative  to  another  electrode  placed  downstream,  fine  hydrogen 
bubbles  were  generated  on  the  uire  and  swept  away  by  the  flou.  Sodium 
sulfate  (NazS0«)  uas  added  to  the  water  at  a  concentration  of  about  150 
gm/m3.  as  recommended  by  Thompson  (1973).  to  increase  its  conductivity. 
The  bubbles  were  small  enough  to  follow  the  flou  sufficiently  close  for 
qualitative  purposes.  The  uire  holder  and  suction  cup-mounted  mirror 
used  to  vieu  the  bubbles  uere  found  not  to  significantly  affect  the  flou 
in  the  region  under  observation.  An  electronic  pulse  generator  (see  the¬ 
sis  of  J.  Simonich  for  details)  enabled  us  to  generate  discrete  time 
lines  of  bubbles  of  variable  uidth  and  spacing. 

The  bubbles  were  illuminated  by  a  500  watt  spotlight  with  a  tungsten 
halogen  bulb.  suspended  above  the  channel  and  shining  down  into  the 
water  (see  Fig.  9).  Best  results  uere  obtained  uhen  the  angle  between 
the  line  of  sight  and  the  light  beam  was  as  small  as  possible.  Hence, 
end  views  uere  preferred  for  photographic  purposes.  End  vieus  also 
accentuated  transverse  motions  of  the  flow  which  were  difficult  to  per¬ 
ceive  in  plan  views.  High-speed  black  and  white  reversal  film  (Kodak  4X) 
forced  two  stops  to  ASA  1600  was  used  to  photograph  the  hydrogen  bub¬ 
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bles. 


2.3  Ml  UlU  KEAmEMENI  JEfiimmiE 


A  0.015  cm  (6  mil)  diameter  TSI  quartz-coated  cylindrical  hot  film  probe 
uas  used  to  measure  velocities  in  the  flou.  The  length  of  the  exposed 
sensor  uas  0.2  cm.  The  probe  uas  run  by  a  TSI  101 0A  anemometer  unit  at 
an  overheat  ratio  of  1.08.  The  anemometers  bridge  output  uas  connected 
to  an  a/d  converter  in  a  DEC  MINC-11  microcomputer,  (see  Fig.  10  for 
1 ayout) . 


2.3.1  Mi  Ella  Calibration 

The  channel  uas  filled  and  set  for  minimum  flourate  by  closing  the  out¬ 
let  valve.  Leakage  through  the  valve  limited  the  minimum  mean  velocity 
to  about  1 . 010.2  cm/s,  determined  by  measuring  the  time-of-f 1 ight  of 
injected  dye.  The  velocity  fluctuated  betueen  0  and  2  cm/s  uith  a  period 
of  about  2  s  due  to  sloshing  in  the  channel.  The  amplitude  of  this 
fluctuation  uas  small  enough  compared  uith  the  range  of  velocities  that 
uould  be  measured  (1  to  20  cm/s)  so  that  the  difference  betueen  the 
average  bridge  voltage  and  the  voltage  that  uould  correspond  to  the 
average  speed  uas  negligeable.  Therefore,  the  velocity  fluctuations  due 
to  the  sloshing  uere  largely  eliminated  by  the  averaging  of  the  bridge 
voltage. 

The  mobile  instrument  cart  toued  the  probe  in  the  upstream  direction. 
The  probe  sensor  uas  located  about  20  cm  belou  the  uater's  surface.  The 
photodiode  speed  indicator  provided  the  cart's  speed  uith  an  uncertainty 
of  ±0.2  cm/s  for  speeds  ranging  from  0  to  20  cm/s.  A  computer  program 
used  this  speed  along  uith  the  computed  average  bridge  voltage  to  create 
a  calibration  table.  (See  Appendix  C  for  a  program  listing.) 


-  30  - 


Figure  10s  Block  diagram  of  data  acquisition  system. 


In  subsequent  calibration  checks,  the  probe  calibration  remained 


steady  within  0.2  cm/s  as  long  as  the  temperature  remained  constant.  3 
or  A  spot  checks  were  always  performed  before  commencing  a  run. 

2.3.2  Velocity  Measurements 

the  difficulty  of  using  hot  films  in  water  is  well -documented  by  Morrow 
£  Kline  (1971).  Before  computer  control  was  instituted,  we  were  unable 
to  control  probe  drift  due  to  dirt  contamination.  In  addition,  the  water 
temperature  generally  rose  about  0.1  degree  Celsius  per  hour  due  to 
heating  from  the  pump  bearings  and  flow  system  losses.  Automating  the 
data  acquisition  process  permitted  a  measurement  technique  which  compen¬ 
sated  for  some  of  the  shortcomings  of  hot  film  probes.  The  main  feature 
of  this  technique  was  using  the  bridge  voltage  obtained  in  the  potential 
core  region  at  a  fixed  reference  location  to  normalize  all  measurements: 
«r 

*ns*»  (2.1) 

*r/ 

where 

en  normalized  voltage  at  measurement  point 
ea  =actual  voltage  at  measurement  point 

er  =voltage  at  reference  point  during  original  calibration 
er'  =voltage  at  reference  point  at  time  of  measurement 

This  procedure  had  the  effect  of  minimizing  changes  in  the  measured 
velocity  due  to  changes  in  the  water  temperature,  total  channel  flow 
rate,  drift  of  the  anemometer  electronics,  and  drift  of  the  probe  out¬ 
put  due  to  dirt  contamination.  Justification  for  this  normalization 
procedure  is  given  in  Appendix  D.  Using  this  procedure,  we  were  able  to 


obtain  results  uith  good  repeatability.  One  beneficial  side  effect  of 
this  measurement  technique  was  the  almost  total  elimination  of  probe 
contamination  due  to  dirt.  Evidentally  the  frequent  movement  of  the 
probe  helped  dislodge  dirt  which  might  otherwise  accumulate  at  stagna¬ 
tion  points.  As  a  result,  probe  drift  due  to  dirt  was  unnoticeable; 
bridge  output  at  the  reference  point  remained  constant  except  for 
changes  due  to  drift  in  the  electronics  and  water  temperature,  which 
were  compensated  for  by  the  normalization  procedure. 

The  normalization  procedure  described  above  requires  the  establish¬ 
ment  of  a  reference  point  at  a  fixed  location  in  the  flow.  Since  the 
cart  could  not  be  moved  by  the  computer,  secondary  reference  points  were 
established  at  each  streamwise  location  where  velocity  profiles  were 
desired.  The  velocities  at  the  secondary  reference  points  were  measured 
uith  respect  to  the  primary  reference  point,  which  was  chosen  to  be  at 
x=400  cm,  y=12  cm,  and  z=2  cm.  Velocity  profiles  at  different  streamwise 
locations  were  then  taken  by  using  the  secondary  reference  points  for 
the  normalization  procedure. 

The  typical  sequence  of  steps  taken  to  acquire  a  velocity  profile  are 
outlined  belou. 

1.  The  time-of-f I ight  of  a  dye  blob  over  a  244  cm  (8  ft)  stretch  of 
the  inlet  was  used  to  determine  the  nominal  velocity  at  the  pri¬ 
mary  reference  point.  Flow  controls  for  the  channel  were  sensi¬ 
tive  enough  to  enable  us  to  repeat  core  velocities  to  within  less 
than  IX  of  the  nominal  14.8  cm/s  on  any  given  day. 


2.  The  velocity  at  the  secondary  reference  point  to  be  used  was 

determined  by  measuring  the  mean  bridge  voltage  at  the  secondary 

reference  point*  normalizing  it  by  the  voltage  at  the  primary 
reference  point*  and  interpolating  between  entries  in  the  cali¬ 
bration  table. 

3.  The  wall  location  was  determined  by  moving  the  probe  against  the 
w?  I  then  gradually  moving  it  away  while  observing  the  bridge 
output.  The  low  velocities  and  high  fluctuations  near  the  wall 
made  it  difficult  to  determine  the  wall  location  to  better  than 
0.05  cm. 

4.  The  computer  was  instructed  to  commence  taking  velocity  measure¬ 
ments.  (See  program  listing  in  Appendix  C.)  The  probe  was  posi¬ 

tioned  by  the  computer  according  to  a  predetermined  set  of  coor¬ 
dinates. 

5.  Prior  to  each  measurement*  the  computer  moved  the  probe  out  to 
the  secondary  reference  point  and  sampled  the  anemometer  at  20  hz 
for  20  s.  The  mean  bridge  voltage  obtained  there  and  the  voltage 
corresponding  to  the  velocity  at  the  secondary  reference  point 
were  used  to  normalize  all  readings  taken  at  the  measurement 
point  according  to  Eq.  2.1. 

6.  In  turbulent  flows*  the  probe  remained  at  the  measurement  loca¬ 
tion  for  5  min  and  was  sampled  at  6  hz  for  a  total  of  1800  sam¬ 
ples.  In  laminar  flows,  a  20  s  average  was  sufficient.  The  pro¬ 
gram  kept  running  tallies  of  the  sum  of  the  velocities  and  sum  of 
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Chapter  III 


RESULTS 

3.1  mflsm  MEAsuammis 

Profiles  of  mean  velocity  and  turbulence  intensity  (u'J)  were  measured 
in  four  cross-planes.  For  purposes  of  qualifying  the  measurement  tech¬ 
nique  and  the  apparatus,  we  also  obtained  some  laminar  profiles. 

Fig.  11  shows  the  coordinate  system  employed  in  describing  positions 
within  the  flow.  The  origin  for  x  was  arbitrarily  chosen  to  coincide 
with  the  beginning  of  the  flat  test  wall,  at  the  exit  of  the  nozzle. 
Curvature  begins  at  x-488  cm  and  ends  at  x=701  cm.  One  cross-plane  is 
located  in  the  inlet  at  x=400  cm,  two  are  located  in  the  curve  at  30° 
and  75°  of  turning  <x=559  cm  and  x=66G  cm,  respectively),  and  the  fourth 
in  the  recovery  region  at  x=76Q  cm.  The  coordinate  y  is  measured  from 
the  wall.  z  is  positive  in  the  upward  direction  and  its  origin  is 
located  near  midspan,  58  cm  from  the  bottom  of  the  channel. 


3.1.1  Laminar  Profiles 

Fig.  12  shows  the  mean  velocity  and  turbulence  intensity  as  measured 
across  the  channel  at  xMOO  cm  at  low  channel  velocities  (6  cm/s)  with 
the  boundary  i^yer  trips  removed.  Cross-channel  profiles  were  taken  by 
traversing  out  from  each  wall  with  some  overlap  of  the  profiles  in  the 
center;  the  degree  to  which  the  readings  from  separate  profiles  at  the 
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Figure  12s  Laminar  profiles. 
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same  location  agree  with  each  other  is  a  good  measure  of  the  repeatabil¬ 
ity  (or  first  order  uncertainty,  see  Moffat  (1978)1  of  the  measurements. 

Several  sources  could  have  contributed  to  the  relatively  high,  appar¬ 
ent  freestream  turbulence  intensity  (u' r«»/Uo=1 . 350 . 

1.  The  reference  voltage4  in  the  anemometer  uas  observed  to  contain 
two  variable  components.  There  uas  a  very  slouly  varying  part 
uhich  seemed  to  be  due  to  changes  in  the  room  temperature.  There 
was  also  a  more  rapidly  varying  component,  uith  a  frequency  on 
the  order  of  a  feu  hertz,  uhich  uas  probably  due  to  noise  in  the 
pouer  supply.  The  slouly  varying  component  uas  removed  by  the 
normalization  procedure.  The  more  rapidly  varying  component  is 
estimated  to  have  added  about  0.4%  to  the  turbulence  intensity. 

2.  The  finite  resolution  of  the  A/D  converter  could  have  introduced 
some  spurious  fluctuations  in  the  signal.  The  resolution  of  the 
A/D  uas  2.5  mV  (12  bits  for  a  range  of  10.24  V).  Since  the  pream¬ 
plifier  feeding  the  A/D  had  a  gain  of  0.5.  the  effect  of  having 
the  least  significant  bit  toggle  betueen  0  and  1  uould  be  to 
cause  a  spurious  turbulence  intensity  of  about  0.5%  to  appear. 

3.  The  sloshing  motion  mentioned  earlier  could  have  been  present  to 
some  degree.  However,  no  changes  in  the  water  level  uere  visible. 


4  This  voltage  is  the  voltage  imposed  across  the  resistance  bridge  in 
the  anemomter  circuit  and  should  not  be  confused  uith  the  voltages 
measured  at  the  reference  points  in  the  flou. 
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Because  the  emphasis  of  this  study  uas  more  on  the  qualitative  features 


i 

| 

of  the  flow,  further  refinements  of  the  measuring  technique  or  the  appa-  I 

ratus  were  not  attempted. 


Fig.  13  shows  laminar  results  normalized  and  compared  with  a  standard 
Blasius  profile.  Agreement  is  fairly  good  except  near  the  wall,  where, 
because  of  the  low  flow  speed,  calibration  errors  were  greater.  Probe 
blockage  effects  might  also  have  been  present  and  affected  the  near-wall 
regions  of  the  profiles. 


3.1.2  Flat  Plate  Turbulent  Boundary  Laver.  Up-Stream  of.  Curved  Region 
Fig.  14  depicts  boundary  layer  profiles  at  the  inlet  to  the  curve, 
x=400,  cm  at  various  transverse  (spanwise)  locations.  In  each  case  the 
origin  of  the  y  coordinate  was  adjusted  to  give  a  best  guess  for  the 
actual  origin  based  on  close-up  examination  of  the  velocity  profile  near 
the  wall.  In  any  case  the  required  shift,  yo,  was  less  than  the  uncer¬ 
tainty  in  the  wall  position.  At  each  location  the  mean  velocity  profile 
is  plotted  in  raw  form  and  in  inner,  wall  coordinates.  The  value  of  the 
friction  velocity  used  to  normalize  the  data  was  obtained  using  a  method 
suggested  by  Coles  (1968)  in  which  a  least-squares  fit  is  performed  to 
fit  the  data  to  Coles'  standard  velocity  profile  by  varying  u*  and  6  for 
y*>50  and  y/6<.75: 

1  n  piy» 

u+=-Xn(y*)  +  C  +  -2sin* —  (3.1) 

x  x  v2sJ 

where  x=0.41  and  C=5.0.  II  is  eliminated  by  evaluating  Eq.  3.1  at  y=6: 

U0  1  fSUo'l  211 

— =-An  -  +  C  +  — .  (3.2) 

U*  X  v  V  j  X 

Eq.  3.2  is  used  with  Eq.  3.1  to  eliminate  II  while  the  least-squares  fit 
is  performed  to  determine  u*  and  6.  The  program  used  was  written  by 
Rangarajah  Jayaraman  at  Stanford.  The  program  also  calculates  the 
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Figure  14*  Profiles  on  inlet  fletplste. 


standard  integral  parameters  by  integrating  a  standard  sublayer  profile 


from  the  wall  to  y*=50  then  integrating  parabolic  segments  fitted  to  the 
data  for  the  remainder  of  the  boundary  layer.  The  values  used  for  the 
sublayer  integrals  are  given  by  Coles  as: 


50 


50 


ju4dy4=540.6» 

0 


Ju42dy4=6546. 

0 


(3.3) 


TABLE  2 


Spanuise  variation  ot  boundary  layer  parameters  at  inlet.  x=400 

cm. 


z 

cm 

6 

cm 

6* 

cm 

0 

cm 

H 

u* 

cm/s 

Re* 

Cf 

.001 

n 

20 

8.39 

1.176 

0.826 

1.423 

0.675 

1349 

4.  164 

0.259 

10 

7.32 

1.122 

0.782 

1.434 

0.676 

1282 

4.146 

0.380 

0 

7.56 

1.160 

0.814 

1.426 

0.686 

1360 

4.101 

0.390 

0 

6.93 

1.057 

0.734 

1.440 

0.684 

1206 

4.222 

0.361 

0 

7.75 

1.170 

0.815 

1.437 

0.666 

1428 

3.979 

0.387 

-10 

6.94 

1.046 

0.729 

1.433 

0.688 

1200 

4.259 

0.338 

-20 

6. 14 

0.931 

0.648 

1.436 

0.702 

1066 

4.434 

0.321 

design 

8.64 

1.198 

0.811 

1.477 

0.655 

1240 

3.706 

0.322 

Some  of  the  more  commonly  used  parameters  appear  in  Table  2  along  with 
their  design  values.  First-order  uncertainty  in  the  measurements  is 
reflected  by  the  scatter  of  the  results  at  2=0  (data  sets  FLAT1,  2.  and 
5).  In  general  the  boundary  layer  was  slightly  thinner  towards  the  bot¬ 
tom  of  the  channel;  this  is  further  supported  by  the  spanuise  measure¬ 
ments  discussed  below.  Because  of  the  emphasis  of  this  study  on  qualita¬ 
tive  results,  we  considered  the  variations  in  boundary  layer  parameters 


acceptable.  For  example*  over  the  central  20  cm  of  the  span  from  z=- 10 
to  z=+10,  0  and  6  had  standard  deviations  of  53  about  their  repective 
means  of  0.775  cm  and  7.30  cm  and  u*  had  a  standard  deviation  of  1.35J 
about  its  mean  of  0.68  cm/s.  In  any  case*  at  all  5  locations  the  bound¬ 
ary  layer  appeared  to  be  typical  of  a  flat  plate  turbulent  layer  grown 
in  a  zero  pressure  gradient  to  a  modest  momentum  thickness  Reynolds  num¬ 
ber  of  approximately  1300. 

Turbulence  intensity  profiles  at  the  inlet  also  appear  normal.  Fig. 
15  compares  turbulence  intensity  profiles  taken  at  x=400  cm  and  z=0  with 
those  obtained  by  Klebanoff  (1954).  Note  that  Klebanoff's  data  were 
taken  at  a  much  higher  Reynolds  number.  Considering  the  uncertainty  in 
our  data,  reflected  by  the  scatter,  agreement  is  adequate. 

3.1.3  Surveys  of  the  Flow  Field 

A  series  of  cross-channel  mean  velocity  and  turbulence  intensity  pro¬ 
files  were  taken  at  the  four  locations  described  earlier:  x=400  cm,  30°, 
75°,  and  x=760  cm.  The  integral  parameters  for  the  test  wall  boundary 
layers  at  these  locations  appear  in  Table  3.  These  parameters  were  cal¬ 
culated  by  fitting  parabolic  segments  to  the  data  and  integrating.  The 
standard  Coles  profile  (Eq.  3.1)  doesn't  appear  suitable  in  the  curved 
and  recovery  regions,  so  cf  was  not  calculated. 

To  simplify  graphing  and  interpreting  the  cross-channel  profiles,  the 
norma)  coordinate  y,  measured  from  the  concave  wall,  has  been  normalized 
by  the  local  channel  width  at  each  successive  streamwise  location. 
Velocities  measured  at  locations  in  the  curve  have  been  normalized 
according  to  the  relation 


TABLE  3 


Streamuise  development  of  integral  parameters. 


Location 

6 

cm 

6* 

cm 

8 

cm 

H 

Data 

Set 

X=400  cm 

7.56 

1.160 

.814 

1.426 

FLAT  1 

6.93 

1.057 

.734 

1.440 

FLAT  2 

7.75 

1.170 

.815 

1.437 

FLAT5 

30° 

10 

1.164 

.874 

1.354 

CRV30 

75° 

merged 

1.141 

.935 

1.221 

CRV75 

X=760  cm 

merged 

.758 

.567 

1.338 

FL760 

Note:  Integral  parameters  at  locations  in  curve 

based  on  radius-normalized  velocities.  S  at  30° 
estimated  from  velocity  profile. 


Mail  (uhere  y=0).  The  potential  wall  velocities  used  to  normalize  the 
turbulence  intensity,  u'.  uere  obtained  in  this  way  by  examining  the 
cross-channel  velocity  profiles.  For  measurements  taken  at  75°  and  x=760 
cm.  uhere  no  potential  regions  mere  evident,  the  maximum  value  of  Un  has 
been  used. 

Fig.  16  shows  the  prof i leg  at  x=400  cm.  about  10  boundary  layer 
thicknesses  upstream  of  the  start  of  curvature;  the  flow  appears  fairly 
symmetrical.  30°  into  the  curve.  (Fig.  17).  the  mean  velocity  and  turbu¬ 
lence  intensity  shoM  some  pronounced  changes.  A  potential  core  is  still 
evident  across  the  middle  fifth  of  the  channel.  However,  the  region  of 
turbulent  activity  on  the  concave  wall  has  enlarged  and  extends  further 
out  into  the  flow  white  a  drop  In  turbulent  activity  Is  evident  near  the 
eonvex  wall.  This  trend  continues  as  shown  by  the  data  at  75*  CFfg.  18). 
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Figure  16s  Cross-ehennel  prof 1  Its  ot  xMOO  cw,  Inlet  section 
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CROSS-CHANNEL  VELOCITY  PROFILE 


Figure  17*  Cross-ehannel  profiles  30*  Into  the  curve  (x=559  e») 


CROSS-CHANNEL  VELOCITY  PROFILE 


The  velocity  profile  is  considerably  fuller  than  that  for  a  flatptate, 
evidence  of  substantial  mixing.  Both  the  velocity  and  turbulence  inten¬ 
sity  profiles  suggest  that  the  boundary  layers  on  the  side  walls  have 
merged  by  this  point. 

In  the  flat  wall  region  at  x=760  cm,  59  cm  downstream  of  the  end  of 
curvature,  the  profiles  (Fig.  19)  show  no  sign  of  recovery.  The  convex 
wall  side  of  the  profile  shows  the  effects  of  the  expansion  in  channel 
width  at  the  90°  location.  As  mentioned  earlier,  the  convex  wall  was 
shaped  to  minimize  the  streamwise  pressure  gradient  on  the  test  surface. 
This  design  requires  an  expansion  of  the  channel  at  the  end  of  curva¬ 
ture.  This  expansion  produces  a  strong  adverse  pressure  gradient  on  the 
convex  wall,  and  thereby  creates  a  small  separation  bubble  or  recircula¬ 
tion  zone.  Injected  dye  indicated  flow  separation  a  feu  centimeters 
upstream  of  90°  and  reattachment  about  20  cm  later,  40  cm  upstream  of 
the  measurement  station  at  x=760  cm.  The  effects  of  the  separation  and 
reattachment  are  evident  in  the  mean  velocity  and  turbulence  intensity 
profiles  in  the  vicinity  of  the  convex  wall  at  x=760  cm.  These  profiles 
strongly  resemble  those  obtained  by  Eaton  (1980)  downstream  of  reattach¬ 
ment  behind  a  backward-facing  step.  The  effects  of  the  separation  and 
reattachment  on  the  convex  wall  appear  not  to  have  any  effect  on  the 
concave  test  wall  profiles. 

Figs.  20-22  show  the  results  of  core  surveys  at  x=400  cm,  30°,  and 
75°.  Profiles  were  taken  by  traversing  in  the  y  direction  at  each  of  21 
spanwise  (z)  locations  spaced  2  cm  apart.  Since  the  ordinate  for  the  U 
profiles  begins  at  10  cm/s,  the  scatter  of  the  data  appears  larger  than 
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Figure  20:  Spanwise  profiles  at  x=400  cm. 
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re  22:  Spanuise  profiles  at  75®. 


the  scatter  in  previous  graphs.  Fig.  20  confirms  the  trend  mentioned 
above,  namely  the  boundary  layer  is  thinner  towards  the  bottom  of  the 
channel;  the  profiles  at  y=2  and  4  cm  show  slightly  higher  velocities 
for  smaller  values  of  z.  The  turbulent  profiles,  however,  appear  fairly 
uniform.  Fig.  21  also  shows  the  boundary  layer  to  be  somewhat  thinner 
towards  the  bottom  of  the  channel.  The  skewing  of  the  turbulence  inten¬ 
sity  profile  at  y=8  cm  is  consistent  with  this  trend.  The  turbulence 
intensity  profiles  at  the  other  y  locations  are  not  as  skewed  because 
the  cross-channel  profile  (Fig.  17)  is  relatively  flat  in  those  regions. 
Fig.  21  reveals  that,  aside  from  the  fluctuations  due  to  scatter  in  the 
data,  no  regular  variation  in  mean  velocity  is  apparent  across  the  span 
of  the  channel  covered  by  the  traverse.  At  75°  (Fig.  22)  the  skewing  of 
the  turbulence  intensity  profiles  has  propagated  to  y-12  cm  as  well, 
although  the  inner  profile-  are  still  uniform.  The  mean  velocities  also 
show  less  uniformity,  and  scatter  in  the  data  is  larger  due  to  the 
higher  turbulence  intensities  and  the  lack  of  a  potential  core  for  con¬ 
sistent  normalization  of  the  anemometer  output. 

Fig.  23  directly  compares  centerline  profiles  along  the  concave  test 
surface  as  the  flow  develops.  The  velocity  profiles  are  seen  to  become 
fuller,  and  the  region  of  high  turbulent  intensity  extends  further  away 
from  the  wall  as  the  flow  proceeds  downstream. 

3.1.4  .EJjuatl  Plots 

In  order  to  determine  if  a  log  region  existed  in  the  boundary  layer  on  a 
concave  wall,  the  mean  velocity  was  plotted  against  the  logarithm  of  the 
distance  from  the  wall  (Fig.  24).  Uncertainty  in  the  probe  position 
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Figure  23:  Streamuise  development  of  velocity  and  turbulence  intensity 
profiles. 


studies  retain  the  usual  logarithmic  region  near  the  wall*  (due  in  part 
to  the  higher  Reynolds  numbers  at  which  their  tests  were  done),  the  wake 
regions  of  their  profiles  all  show  a  tendency  to  undershoot  the  usual 
Lau  of  the  Mall.  Fig.  25  shows  mean  velocity  profiles  taken  from  the 
study  of  impulsive  curvature  of  Smite  et  al  (1979).  (Their  channel  had 
a  sharp  30°  bend  in  it  at  x=0  mm.  Profiles  in  (a)  were  taken  at  a  span- 
wise  location  where  cf  was  a  maximum:  profiles  in  (b)  were  taken  where 
Cf  was  a  minimum.  The  scale  on  the  left  applies  only  to  the  lowest 
curves.  The  left-hand  end  of  the  log-law  line  marks  U0/u*s15  for  the 
other  curves,  which  are  staggered.)  In  tlie  30°  profile  of  the  present 
study,  it  appears  as  though  the  region  from  y=0 . 5  to  2.0  cm  may  actually 
be  part  of  this  undershooting  and  should  not  be  thought  of  as  the  usual 
log  region.  Further  discussion  of  this  behavior  is  given  in  the  next 
chapter. 

3.2  VISUAL  RESULTS 

Extensive  visualization  surveys  of  the  flow  were  taken  using  dyes  and 
hydrogen  bubbles.  In  addition  to  surveying  the  turbulent  flow  field,  we 
also  visualized  a  laminar  boundary  layer  on  the  concave  wall  to  see  if 
we  could  obtain  the  classical  Taylor-Gortler  vortex  flow. 

3.2.1  Laminar  Flow 

In  order  to  achieve  the  proper  conditions,  the  velocity  had  to  be  low 
enough  so  that  transition  did  not  ooour.  However,  too  low  a  velocity 
resulted  in  the  development  of  a  strong  secondary  flow  pattern  due  to 
the  boundary  layer  on  the  channel  floor.  This  seoondary  flow  dominated 
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the  entire  curved  portion  of  the  channel*  resulting  in  flow  angles 
approaching  45°  on  the  test  wall.  By  running  the  channel  fast  enough  to 
maintain  a  turbulent  boundary  layer  on  the  floor  but  slow  enough  to  pre¬ 
vent  transition  on  the  test  wall/  we  were  able  to  produce  a  nominally 
two-dimensional  laminar  boundary  layer  on  the  test  wall.  The  mean 
velocity  under  these  conditions  was  about  5  cm/s,  and  the  boundary  layer 
was  about  4  cm  thick  at  the  start  of  curvature. 

The  flow  was  visualized  by  injecting  dye  through  the  wall  slots.  Oye 
coming  out  of  the  first  dye  slot,  chosen  to  be  upstream  of  curvature, 
formed  a  uniform  transparent  sheet.  As  the  sheet  entered  the  curve,  it 
coalesced  into  several  very  narrow  discrete  streamers,  the  result  of 
Taylor-Gortler  vortices.  Oye  injected  through  a  second  slot  30°  into  the 
curve  coalesced  within  a  few  centimeters.  The  first  time  this  test  was 
performed,  the  streamers  were  evenly  spaced  about  10  cm  apart  and  tended 
to  be  quite  stable  (see  Fig.  26a).  In  a  subsequent  experiment,  the  spac¬ 
ing  was  not  nearly  as  uniform.  The  reasons  for  this  nonrepeatability  are 
not  clear  at  this  time,  but  the  transitional  nature  of  this  flow  would 
seem  to  make  it  very  sensitive  to  upstream  flow  conditions.  During  the 
second  experiment,  the  channel  had  to  be  narrowed  3  to  4  centimeters  in 
order  to  impose  a  favorable  pressure  gradient  to  delay  transition. 

An  interesting  feature  revealed  by  the  dye  was  a  secondary  instabil¬ 
ity  which  sporadically  developed  out  of  the  primary  Taylor-Gortler 
cells.  Occasionally,  the  low-speed  fluid  brought  away  from  the  wall  by 
the  roll  cells  would  Interact  with  the  outer  flow,  forming  a  horseshoe 
vortex  street  (see  Fig.  26b>.  Similar  structures  were  observed  in  the 
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studies  of  Bippes  £  Gortler  (1972).  Aihara  £  Sonoda  (1981).  and  Aihara  £ 
Koyama  (1981).  The  vortex  streets  also  resemble  the  coherent  structures 
studied  by  Perry  8  Urn  (1978)  in  their  work  on  eoflouing  jets  and  wakes. 
The  spanuise  width  of  the  horseshoe  vortices  was  roughty  5-7  cm  and  the 
spacing  between  successive  vortices  was  roughly  4-5  cm.  Evidentally, 
the  tube-like  shear  layer  between  the  high-  and  low-speed  fluids  was 
rolling  up  as  a  result  of  a  Helmholtz-type  instability.  The  resulting 
vortices  would  drift  downstream  at  velocities  close  to  that  of  the  free- 
stream.  occasionally  pairing  before  breaking  down  into  turbulence. 
Although  these  vortex  streets  usually  formed  within  a  feu  centimeters  of 
the  wall,  some  of  them  appeared  almost  at  midchannel. 


The  vortex  streets  appeared  at  irregular  intervals,  forming  for 
several  seconds  then  washing  away.  Usually,  the  streamers  appeared  to 
develop  a  transverse  oscillation  of  increasing  magnitude  which  led  to 
breakdown.  Only  "quiet  looking"  streamers  seemed  to  develop  into  vortex 
streets.  Since  the  primary  aim  of  this  work  was  the  study  of  turbulent 


boundary  layers,  we  did  not  pursue  the  investigation  of  laminar  insta¬ 
bility  any  further. 


3.2.2  ivrfrulgnt  Dan 

Flow  through  the  nozzle  appeared  to  be  acceptably  rectilinear  and  lami¬ 
nar.  However,  visualization  of  the  flow  over  the  boundary  layer  trips 
revealed  a  curious  phenomena.  The  recirculation  region  on  the  upstream 
side  of  the  trip  was  apparently  unstable  in  the  spanwise  direction. 
Every  feu  centimeters,  the  fluid  trapped  in  the  recirculation  vortex 
would  find  its  way  over  the  trip  as  shown  in  Fig.  27.  The  recirculation 


-  62  - 


Figure  27>  Schematic  of  flow  over  boundary  layer  trip. 


region  therefore  appeared  segmented  and  three  dimensional.  The  pattern 
of  "breakthroughs"  was  very  steady*  although  large  disturbances  uould 
cause  the  pattern  to  change  and  "lock  in"  to  a  new  arrangement.  More 
detailed  flou  visualization  revealed  that  the  actual  volume  of  fluid 
involved  in  the  breakthroughs  uas  very  small  and  that  the  effects  of 
transition  dounstream  of  the  trip  obliterated  all  visual  evidence  of  the 
nonuniformity.  No  signs  of  the  nonuniformity  were  apparent  in  the  span- 
wise  velocity  profiles  taken  farther  downstream. 


Because  turbulent  flow  is  strongly  three-dimensional  and  time- 
dependent,  it  is  almost  impossible  to  completely  document  it  for  future 
analysis.  The  best  we  could  do  was  to  take  motion  pictures  of  the  flow, 
selecting  the  camera  angle  to  optimize  the  visual  effects  of  the  three- 
dimensionality.  A  summary  of  our  results  is  contained  in  the  film,  "Tur¬ 
bulent  Boundary  Layers  on  Concave  Mails,"  available  on  loan.  The  film 
contains  footage  of  dye  and  hydrogen  bubble  visualization  of  the  flow  in 
the  inlet,  curved,  and  recovery  regions.  A  detailed  scene-by-scene  nar¬ 
rative  may  be  found  in  Appendix  F.  Representative  still  photos  are  shown 
in  Figs.  2B,  and  30-33.  The  geometry  of  the  hydrogen  bubble  photos  is 
depicted  in  Fig.  29.  The  solid  horizontal  lines  represent  the  locations 
of  the  hydrogen  bubble  wires,  which  are  invisible  in  the  actual  photos. 
The  dotted  lines  are  located  at  the  same  y  values  as  the  wires,  but  at 
the  mean  x  location  of  the  bubble  lines.  Each  group  of  hydrogen  bubble 
pictures  shows  the  wire(s)  at  different  distances  from  the  wall.  The 
numbers  to  the  right  give  the  y  values  of  the  wires  in  centimeters.  In 
the  case  of  the  discrete  bubble  lines,  the  photos  show  the  bubbles  about 
0.57  sec  after  the  wire  was  pulsed.  Due  to  the  mirror  arrangement  used, 
the  flow  is  predominantly  out  of  the  page  but  slightly  canted  towards 
the  lower  left  hand  corner.  Uhat  follows  is  a  description  of  the  flow 
field  based  on  still  photos,  movies,  and  actual  observations. 

The  flat  plate  boundary  layer  upstream  of  the  curve  at  x=460  cm 
exhibits  the  usual  pattern  of  streaks  and  bursts  discovered  by  Runsta- 
dler  £  Kline  (1959).  With  the  bubble  wire  at  y=0.5  cm,  (y4=30)  the 
streaks  appear  as  narrow  regions  of  low  velocity  fluid  which  oscillate 
transversely.  The  mean  streak  spacing  was  determined  by  manually 


analyzing  dye  pictures  and  found  to  be  1 . 51 ±0 .16  cm.  Using  the  friction 
velocity  calculated  at  midspan  of  0.68  cm/s.  this  corresponds  to  a  of 
103i11,  which  agrees  well  with  the  generally  accepted  value  of  100 
[Cantwell  (1981)].  In  addition,  longitudinal  vortices  near  the  wall  are 
occasionally  seen.  The  typical  length  scale  of  these  structures  and  the 
other  eddies  present  in  the  flow  is  about  1  to  2  cm.  The  intensity  of 
turbulent  activity  seems  to  diminish  monotonical ly  as  one  moves  away 
from  the  wall;  the  bubbles  at  y- cm  appear  to  be  steady  and  almost 
laminar  except  for  occasional  small  travelling  disturbances  associated 
with  the  intermi ttency  at  the  edge  of  the  boundary  layer. 


Flgurt  30:  Hydrogen  bubbles  st  x=400  cm.  inlet. 
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Figure  31:  Hydrogen  bubbles 


Figure  33«  Hydrogen  bubbles  st  x=730  cm.  recovery. 


By  about  20°  into  the  curve,  the  boundary  layer  appears  markedly  dif¬ 
ferent  from  that  on  the  flat  plate.  The  boundary  layer  is  about  10  cm 
thick  at  this  station.  Near  the  wall,  streaks  are  still  present,  but 
they  are  frequently  swept  aside  or  squashed  by  large-scale  inflows  or 
sweeps,  with  spanwise  dimensions  on  the  order  of  4  cm  and  streamwise 
lengths  of  about  10  to  IS  cm.  The  streaks  and  their  surrounding  fluid 
are  pushed  into  piles  and  the  usual  near-wall  bursting  process  appears 
to  be  mostly  suppressed.  Small  longitudinal  vortices  are  occasionally 
seen  at  the  interface  between  the  high-speed  inflows  and  the  low-speed 
wall  fluid.  Bubbles  farther  from  the  wall  show  that  the  low-speed  col¬ 
lections  of  fluid  are  squeezed  or  ejected  away  from  the  wall  by  the 
sweeps  and  frequently  make  their  way  out  to  the  potential  core.  Dual 
bubble  lines  4  cm  apart  in  the  y-direction  show  that  vertical  motions 
are  often  correlated  across  the  boundary  layer.  Planar  visualization 
with  one  bubble  wire  pulsed  continuously  near  the  wall  or  near  the  edge 
of  the  boundary  layer  (y=9  cm)  shows  that  the  ejections  and  sweeps  occur 
randomly  in  space  and  time.  The  impression  one  gets  from  viewing  the 
actual  flow  is  that  the  sweeps  are  somewhat  more  coherent  and  appear  to 
be  responsible  for  the  ejections,  which  seem  to  form  as  a  result  of  the 
sweeps.  Small  vortices  are  also  present  at  the  edges  of  the  outflowing 
bulges  of  the  ejections  as  they  penetrate  the  core.  Careful  observation 
i  shows  gentle  streamwise  circulations  which  extend  over  the  entire  bound¬ 

ary  layer  taking  place.  However,  these  weak  roll  cells,  reminiscent  of 
but  weaker  than  the  roll  cells  indentified  in  the  related  work  by  John¬ 
ston  et  al  (1972)  in  a  rotating  channel,  appear  to  be  more  a  result  of 
the  sweeps  and  ejections  than  Independent  flow  modules  themselves. 
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As  the  flow  proceeds  downstream,  its  qualitative  features  as 
described  above  remain  the  same.  However,  the  boundary  layer  gets 
thicker,  and  the  crossplane  movements  of  fluid  increase  in  speed  and 
extent.  Intermittent  values  of  v  and  u  are  often  on  the  order  of  A  cm/s 
or  25JC  of  the  mean  channel  velocity. 

Near  the  downstream  end  of  curvature  and  in  the  recovery  region  fol¬ 
lowing  the  concave  wall  an  additional  flow  characteric  becomes  evident. 
Large  (10  cm)  patches  or  globs  of  laminar-like  fluid  appear  throughout 
the  boundary  layer,  interspersed  within  the  large-scale  motions  remain¬ 
ing  from  the  preceeding  curvature. 

To  assist  in  interpreting  these  visual  results  a  computer  simulation 
of  streaklines  and  timelines  was  performed.  This  simulation  is  discussed 
in  the  following  chapter.  Interpretations  of  these  observations  are  then 
given. 


3.3  SUMMARY  HE  flAIM  RESULTS 

Although  there  was  minor  spanwise  nonuniformity  in  the  inlet  section 
boundary  layer,  the  flow  was  considered  acceptable  for  the  present  qual¬ 
itative  study.  The  boundary  layer  upstream  of  the  curved  section  appears 
to  be  a  normal,  zero-pressure-gradient  flat  plate  boundary  layer  with  a 
momentum  thickness  Reynolds  number  of  about  1300.  The  streak  spacing 
agrees  with  the  generally  accepted  value  of  100  wall  units,  and  the 
boundary  layer  appears  normal  in  the  flow  visualizations. 

Once  in  the  curved  section,  the  appearance  of  the  boundary  layer  is 
changed  by  the  presence  of  large-scale  sweeps  and  ejections.  The  sweeps 
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arc  typically  about  half  a  boundary  layer  thickness  in  span  and  appear 
to  be  coherent  over  the  height  of  the  boundary  layer.  The  streamuise 
extent  of  the  sweeps  is  difficult  to  ascertain  using  the  present  Methods 
of  visualization  but  is  estimated  to  be  several  boundary  layer  thick¬ 
nesses.  The  sweeps  appear  randomly  in  space  and  time  over  the  concave 
wall.  Although  no  stationary  pattern  is  observed,  the  sweeps  are  spaced 
about  one  boundary  layer  thickness  apart,  on  the  average,  in  the  span- 
wise  direction.  Fluid  near  the  wall  is  displaced  by  the  sweeps  and 
ejected  away  from  the  wall.  Bursting  appears  to  be  inhibited  in  the 
regions  underlying  the  sweeps. 

The  random  occurrence  of  the  sweeps  and  ejections  is  supported  by 
velocity  and  turbulence  intensity  profiles  obtained  by  traversing  in  the 
spanuise  direction;  nothing  in  these  profiles  suggests  that  the  large- 
scale  sweeps  and  ejections  occur  in  preferred  locations.  Traverses 
across  the  height  of  the  boundary  layer  reveal  fuller  velocity  profiles 
and  higher  values  of  turbulence  intensity,  especially  in  the  outer 
layer.  These  results  confirm  the  visual  observations,  which  indicate  an 
increase  in  mixing  due  to  the  large-scale  sweeps  and  ejections. 

The  modification  of  the  velocity  profiles  by  concave  curvature  fol¬ 
lows  a  trend  apparent  in  the  studies  of  earlier  investigators.  Concave 
curvature  seems  to  cause  the  outer  layer  profile  to  undershoot  the  usual 
Law  of  the  Hall,  resulting  in  an  attenuation  or  elimination  of  the  usual 
logarithmic  region. 
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Chapter  IV 


DISCUSSION 

4. 1  SIMULATION  AND  INTERPRETATION  HE  FLOW  VISUALIZATION 
Flow  visualization  of  an  unsteady  ftou  by  fluid  marking  has  a  major 
drawback,  namely,  the  fluid  markers  show  a  time  integrated  picture  of 
the  flowfield  since  the  current  position  of  each  marker  depends  on  its 
prior  history.  In  steady  flows  streaklines,  streamlines  and  pathlines 
are  identical,  but  this  cannot  be  said  of  unsteady  flows.  The  human  eye 
and  mind  are  able  to  integrate  and  interpret  streak  and  pathline  pic¬ 
tures  to  a  limited  extent.  Thus  when  viewing  a  flow,  some  useful  inform¬ 
ation  regarding  the  instantaneous  velocity  distribution  can  be 
extracted.  However,  there  is  still  a  tendency  for  the  eye  to  make  erro¬ 
neous  conclusions  based  on  some  of  the  shapes  formed  by  the  markers.  In 
order  to  help  interpret  the  visualization  data  in  the  present  study  and 
to  try  out  some  simple  models  of  the  flow  field,  a  computer  program  was 
developed  which  is  capable  of  tracing  fluid  markers  (e.g.  streak  and 
timelines)  in  a  prescribed  three-dimensional,  time-dependent  velocity 
field.  The  results,  when  displayed  graphically  and  compared  with  actual 
visualizations  of  the  flow,  give  an  idea  of  what  the  actual  velocity 
fields  may  be  1  ike. 

Given  a  known  Eulerian  velocity  field,  for  example  u(x,t),  and  a 
fluid  marker  which  is  identified  by  it  intial  position  at  time  t-0,  X0, 

the  position  of  the  marker,  X(t),  is  given  by  solving  the  differential 


equation i 


— =u(X,t), 
dt 


(4.1) 


subject  to  the  initial  conditions  X(0)=&o> 

In  the  present  study,  a  velocity  field  resembling  a  Taylor-Gortler 
velocity  field  was  prescribed: 

u=U0-v/y  (4.2) 

-2V0hy 

v= - cos(kzz)cosIkx(x-ut) ]  (4.3) 

h2+y2 

2Voh(h2-y2 ) 

w= - sin(kzz)cos[k*(x-ut)  1  (4.4) 

k*(h2+y2)2 

where 

2n 

k i= — ,  i=x,y,  or  z.  (4.5) 

Ai 

This  flow  field  satisfies  continuity,  but  it  is  not  irrotat .onal .  The 
intent  here  is  not  to  propose  a  realistic  model  of  an  actual  flow  field 
but  to  see  how  fluid  markers  behave  in  a  hypothetical  flow  field  which 
has  some  similarities  to  an  actual  flow. 

The  flow  described  b>  Eqs.  4. 2-4. 4  consists  of  an  array  of  counterro- 
tating  roll  cells  adjacent  to  a  wall,  each  of  width  Az/2.  with  solid 
body  rotation  in  the  center  of  each  cell  at  height  h.  At  the  corners  of 
each  cell  there  is  stagnation  point  flow  (see  Fig.  34).  The  maximum 
vertical  velocity.  Vo  and  occurs  at  a  height  h  on  the  vertical  line 
separating  adjacent  cells.  The  streamwise  (u)  velocity  consists  of  a 
constant  velooity.  Uo»  plus  a  shear,  v'.  The  shear  is  used  in  an  attempt 
to  reproduce  the  fact  that  flow  moving  towards  the  wall  has  a  higher  u 
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Figure  34:  Cross-plane  streamlines  for  simulated  velocity  field. 


velocity  than  flou  moving  auay  from  the  wall.  The  cosine  terms  at  the 
ends  of  the  expressions  for  v  and  h  modulate  the  amplitudes  of  the 
cross-stream  velocities  so  that  the  cells  are  of  limited  streamuise 
extent  and  travel  downstream  as  coherent  entities,  like  travelling 
waves. 

The  differential  equations  represented  by  Eq.  4.1  were  solved  using  a 
fourth  order  Runge-Kutta  ODE  solver  given  jy  White  (1974).  Subroutines 
were  written  to  display  the  results  in  true  perspective. 

The  resulting  velocity  field  causes  computed  time-streak  markers  to 
behave  as  shown  in  Figs.  35  and  36.  Each  cell  is  5  cm  wide  (A*=10  cm) 
and  10  cm  long  (Ax=20  cm)  with  a  core  set  at  h=2.5  cm  from  the  will 
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Figure  35:  Evolution  of  simulated  time/streak  markers  with  time. 


(units  are  arbitrary).  The  base  plane  is  20  cm  square.  Uo>  the  velocity 
at  y=0,  is  11  cm/s,  and  the  shear  v'  is  chosen  so  that  the  velocity  at 
the  height  of  the  marker  generator  is  10  cm/s  in  all  cases.  In  this  way. 
there  is  less  of  a  vertical  velocity  gradient  for  markers  placed  farther 
from  the  wall,  as  is  the  case  in  the  actual  flow.  The  maximum  vertical 
velocity,  V0,  is  4.0  cm/s.  The  values  for  these  parameters  roughly 
approximate  the  actual  observed  characteristics  of  the  flow  uhile  at  the 
same  time  simplifying  the  display  of  the  results.  Fig.  35  shows  a 
sequence  of  views  at  consecutive  times  with  the  generator  at  y=2.5  cm. 
Note  the  development  of  the  the  "structures"  as  they  proceed  downstream. 
Fig.  36  has  the  generator  at  y=0.5,  1.0,  5.0  and  10  cm,  all  at  the  same 
time.  If  we  imagine  some  random  turbulence  superimposed  on  the  simu¬ 
lated  flow  field  (Figs.  35  £  36),  the  result  is  qualitatively  similar  to 
the  flow  actually  observed  (Figs.  31-33). 

The  simulation  exhibits  several  of  the  features  observed  in  the 
actual  flow  field  and  helps  explain  how  some  of  those  features  may 
appear.  One  of  the  more  striking  resemblances  of  the  simulation  to  the 
actual  flow  is  the  narrowness  of  the  outflows  (ejections)  for  markers 
placed  near  the  wall  (Fig.  36a).  The  simulation  demonstrates  that  such 
outflows  are  simply  a  result  of  the  combined  movements  of  the  markers 
and  are  not  necessarily  indicative  of  narrow,  jet-like  ejections  in  the 
velocity  field.  Another  feature  the  simulation  helps  clarify  is  the 
streamuise  development  of  the  large-scale  structures.  The  markers  make 
it  appear  as  though  the  structures  are  larger  downstream  of  the  wire 
when  in  fact  we  know  the  prescribed  velocity  field  is  uniform  aside  from 
the  periodic  fluctuations.  The  assumption  of  a  slow  streamwise 
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development  of  the  actual  flou  field  Is  supported  by  the  observation 
that  similar  structures  appear  at  all  locations  in  the  curve.  The  varia¬ 
tions  in  shape  of  the  different  timelines  in  the  real  flow  is  due  to  the 
random  distribution  of  the  "cells"  and  the  different  positions,  relative 
to  these  "cells",  at  which  bubbles  are  injected.  In  fact,  "cells"  is 
probably  not  a  good  word  to  describe  the  regions  of  inflow  and  outflow 
in  the  actual  flou  since  planar  visualization  (in  the  x-z  plane)  shows 
them  to  be  multiply-connected  domains  of  irregular  shape  and  size. 


4.2  PHYSICAL  £ QMS IP £B AXIOMS 

Most  previous  investigators  [e.g.  Tani  (1962),  So  6  Mellor  (1975),  Hoff¬ 
man  £  Bradshaw  (1981)1  have  reported  stable  longitudinal  flou  structures 
in  the  boundary  layer  on  a  concave  wall  and  have  concluded  that  a  sta¬ 
tionary  array  of  counterrotating  turbulent  vortices  exist  there,  analo¬ 
gous  to  the  Taylor-Gortler  vortices  observed  in  a  laminar  boundary 
layer.  Clearly,  our  flou  is  also  dominated  by  large-scale  structures, 
the  difference  being  that  our  structures  have  less  coherence  and  are 
distributed  randomly  over  the  wall.  To  see  why  this  might  be  so,  let's 


examine  current  knowledge  about  the  fluid  dynamics. 


Various  stability  criteria  lead  one  to  expect  flou  with  streamline 
curvature  to  be  unstable  to  disturbances  in  the  plane  of  curvature  if  a 
certain  velocity  gradient  exists.  Gortler  (1954)  and  others  have  done 
detailed  analyses  to  determine  the  exact  conditions  for  instability  in  a 
laminar  boundary  layer  assuming  that  the  initial  disturbances  would  take 
the  form  of  vortices.  It  is  important  to  keep  in  mind  the  other  assump¬ 
tions  made  in  these  analyses.  The  initial  boundary  layer  is  steady, 
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laminar,  and  tuo-dimensional . 


The  initial  disturbances  are  assumed 


small  so  that  linearization  of  the  equations  is  possible.  The  theory 
says  nothing  about  how  the  flow  will  behave  as  the  disturbances  grou.  In 
fact,  our  oun  laminar  flow  experimental  studies  as  well  as  those  of  Bip- 
pes  (1972),  Aihara  £  Sonoda  (1981),  and  others,  have  demonstrated  that 
as  the  vortices  develop,  a  secondary  instability  occurs  in  the  shear 
layer  between  the  lou-speed  fluid  transported  away  from  the  wall  and  the 
high-speed  outer  layers,  leading  to  eventual  breakdown  and  turbulence. 

Before  moving  on  to  discuss  turbulent  flow,  it  may  be  useful  to  exam¬ 
ine  the  physical  basis  for  the  instability  of  curved  flows.  Newton's 
first  Law  says  that  an  object  in  motion  in  a  straight  line  will  continue 
moving  in  a  straight  line  unless  acted  upon  by  a  force.  In  the  absence 
of  forces,  a  particle  moving  initially  tangent  to  a  wall  with  concave 
curvature  will  approach  the  wall  and  eventually  collide  with  it.  The 
same  will  happen  even  if  there  are  forces  on  the  particle,  if  those 
forces  are  insufficient  to  turn  the  particle  (i.e.  F  <  mv*/R).  Imagine  a 
constant  density  fluid  flow  along  a  concave  wall  which  contains  a  fluid 
particle  with  an  abnormally  high  velocity.  Uith  its  higher  momentum,  the 
particle  will  move  in  a  path  with  less  curvature  than  its  neighbors  and 
approach  the  wall.  Continuity  requires  that  the  surrounding  lower-speed 
fluid  move  away  from  the  wall.  How  fast  need  the  particle  move  to  devi¬ 
ate  from  its  neighbors  in  an  unstable  manner?  For  inviscid,  constant 
density  flow,  Rayleigh's  stability  analysis  and  the  Navier-Stokes  equa¬ 
tions  show  that  potential  flow  is  the  condition  for  neutral  stability, 
i.e.  Ur=constant.  If  the  local  velocity  gradient  exceeds  that  for  a 
potential  flow,  the  flow  is  unstable.  For  viscous  flows,  the  conditions 
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are  not  so  simple  as  GSrtler  demonstrated#  but  the  basic  idea  is  the 
same. 

In  turbulent  boundary  layers  the  fluid  velocity  varies  randomly  in 
space  and  time#  although  the  mean  is  stationary.  If  a  point  function 
describing  the  local  stability  uere  defined#  its  value  would  also  vary 
randomly  throughout  the  flow.  Hence,  in  a  turbulent  flow  one  would 
expect  at  any  instant  to  see  randomly  distributed  regions  (in  all  three 
dimensions)  of  stability  and  instability.  The  size  of  the  regions  would 
depend  in  some  way  on  the  characteristics  of  the  fluid  and  the  statis¬ 
tics  of  the  velocity.  furthermore,  the  velocity  statistics  would  be 
continuously  altered  by  the  mixing  of  high-  and  low-speed  fluid  ele¬ 
ments.  Added  to  this  picture  is  the  structure  of  the  boundary  layer 
itself,  manifested  in  streaks  and  bursts  as  in  the  basic  case  of  flow 
over  a  flat  wall.  Nothing  in  this  picture  suggests  the  development  of  a 
stationary  pattern  of  Taylor-GSrtler  vortices. 

Mhy  have  so  many  other  investigators  discovered  stationary  flow  pat¬ 
terns,  e.g.  Taylor-Gortler-1 ike  cells  of  long  streamwise  extent?  Rig 
dependency  must  be  seriously  considered  when  dealing  with  instability 
phenomena.  Remember  that  it  is  variations  in  velocity  throughout  the 
flowfield  which  govern  the  stability.  It  is  therefore  important  to  con¬ 
sider  (a)  how  experimenters  deal  with  the  streamwise  pressure  gradient7 
associated  with  the  onset  of  curvature,  (b)  how  they  control  inlet  con¬ 
ditions  and  secondary  flows  on  endualls,  and  (c)  the  geometry. 


7  Aihara  £  Sonoda  (1981)  showed  that  the  onset  of  Taylor-GSrtler  vor¬ 
tices  in  a  laminar  boundary  layer  was  delayed  by  a  favorable  stream- 
wise  pressure  gradient. 


particularly  the  aspect  ratio,  of  the  channel  they  employ  to  investigate 
the  phenomenon.  It  is  felt  that  the  control  of  flou  nonuniformities,  by 
inlet  screens  for  example  (see  comments  in  previous  section),  can  have  a 
very  strong  influence  on  the  final  outcome.  Very  small  stationary  span- 
uise  disturbances  generated  by  pore  size  nonuniformities  in  the  inlet 
screens  or  other  elements  can  propagate  far  dounstream.  Depending  on 
their  relative  amplitudes  and  distributions.  these  disturbances  can  in 
fact  trigger  standing  Taylor-Gortler-1 ike  cells  in  a  turbulent  boundary 
layer,  as  demonstrated  by  Hoffman  C  Bradshaw  (1981).  In  conclusion,  it 
appears  that  some  of  the  previous  experimental  work  may  in  fact  be  con¬ 
taminated  by  influences  outside  the  control  of  the  experimenters.  In  a 
number  of  cases  these  influences  may  be  creating  the  stationary  Taylor- 
Gortler-like  cells  that  are  apparently  being  observed. 

4.3  £Hm£M  ryiy  mm. 

The  physical  arguments  given  above,  the  results  of  the  flou  simulation, 
and  the  results  of  actual  observation  of  the  flou  give  rise  to  the  fol¬ 
lowing  suggested  physical  model  for  the  flou.  This  model  is  advanced  as 
a  hypothesis  which  needs  further  verification. 

Pockets  of  relatively  high-speed  fluid  force  their  way  towards  the 
concave  wall  as  sweeps,  pushing  aside  low-speed  fluid.  The  sweeps  origi¬ 
nate  as  the  random  fluctuations  in  velocity  associated  with  the  large 
eddies  in  the  outer  flow.  The  spanwise  width  of  the  sweeps  is  on  the 
order  of  half  the  boundary  layer  thickness,  or  about  500  y*  units. 
Although  they  occur  randomly  throughout  the  boundary  layer,  the  spanwise 
spacing  of  the  sweeps  is  on  the  average  comparable  to  the  boundary  layer 
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thickness.  Furthermore*  although  the  sweeps  may  initially  be  localized* 
inertial  effects  and  the  presence  of  a  mean  velocity  gradient  which  is 
less  stable*  near  the  wall  would  tend  to  carry  the  sweeps  all  the  way  to 
the  wal 1 . 

Sweeps  impinging  on  the  wall  set-up  regions  of  stagnation-like  flow 
with  appreciable  transverse  velocities.  The  sweeps  push  aside  the  low- 
speed  fluid  near  the  wall  along  with  the  streaks  embedded  in  it.  The 
severe  distortion  and  acceleration  that  the  fluid  undergoes  may  inhibit 
streak  bursting.  Bursting  is  a  phenomenon  associated  with  the  presence 
of  streamwise  vorticity  and  a  pattern  of  low-speed/high-speed  streaks  in 
the  wall  layer.  The  streaks  may  not  get  a  chance  to  form  fully  before 
they  are  disrupted  by  the  large-scale  sweeps.  However,  the  turbulence 
intensity  remains  high  in  the  wall  layers  because  of  contributions  from 
the  larger  scales  of  motion  and  the  shearing  action  between  the  sweeps 
and  the  low-speed  wall  layers  which  get  caught  between  adjacent  high¬ 
speed  sweeps  and  are  ejected  away  from  the  wall. 

The  ejection  of  low-speed  fluid  into  the  higher  velocity  regions 
appears  to  increase  the  turbulence  production  there  with  a  consequent 
increase  in  turbulence  intensity  far  from  the  wall.  What  appears  to  be  a 
Helmholtz  instability  occurs  in  the  narrow  shear  layers  formed  between 
ejections  and  sweeps;  this  may  be  the  mechanism  that  enhances  turbulence 
production  throughout  the  boundary  layer  and  could  possibly  be  a  process 
that  keeps  the  overall  turbulence  levels  higher  than  in  a  flat  wall 
boundary  layer  even  though  normal  bursting  appears  to  be  suppressed  as 


*  Less  stable  in  the  sense  that  it  represents  a  larger  deviation  from 
neutral  stability,  ursoonstant. 


described  above. 


It  thus  appears  that  these  large-scale  sweeps  and  ejections  augment 
the  effective  Reynolds  stress  throughout  the  boundary  layer*  even  though 
bursting  near  the  wall  may  be  suppressed.  If  this  is  true,  one  model  or 
picture  of  the  structure  might  be  one  of  a  magnified  wal l-layer-1 ike 
sueep/e jection  mechanism.  Thus  the  boundary  layer  as  a  whole  may  have  a 
structure  similar  to  that  in  the  2one  0<y*<30  in  a  flat  plate  boundary 
layer  but  on  a  much  larger  scale.  This  would  explain  the  modification 
or  absence  of  a  logarthmic  zone  in  boundary  layers  on  concave  walls. 
This  conjecture  needs  further  study. 

The  actual  flow  is  of  course  complicated  by  the  random  variations  in 
sweep  size,  intensity,  and  location.  and  the  interaction  between  the 
high  and  low-speed  fluid  elements.  In  fact,  the  sweeps  and  ejections  are 
actually  more  like  multiply-connected.  intertangled  regions  in  space 
which  are  constantly  deforming  as  they  drift  downstream.  Because  of 
continuity,  sweeps  and  ejections  must  give  rise  to  localized  regions 
where  the  flow  is  rotating.  However,  the  extent,  shape,  and  strength  of 
these  rotating  bodies  of  fluid  make  it  difficult  to  think  of  them  as 
coherent  structures  in  themselves.  For  this  reason,  the  terms  longitu¬ 
dinal.  xarljsc?  and  rail  cells  aita  inappropriate  iat  describing  JMs 
flan-  il&K£¥££>  iba  structural  el.enen.ts  consisting  si  SMCCH5  &M  Cita¬ 
tions,  random  though  they  may  fee*  atCffl  la  £CBCfll  themselves  sn  C  gjiaBir. 
atEiadia  feaali  end  marc  adeguateiv  characterise  Ifee  llaw- 


*  The  turbulence  measurements  of  So  C  Mel  lor  (1972)  and  Hoffman  £  Brad¬ 
shaw  ( 19S1 )  support  this  remark,  but  additional  quantitative  measure¬ 
ments  are  needed  to  distinguish  between  turbulence  production  by 
bursting  and  turbulence  production  by  sueep/e jeotion  Interactions. 
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It  should  be  emphasized  that  the  flow  over  the  concave  wall  in  this 
study  is  a  developing  one.  The  imposition  of  concave  curvature  on  the 
flat  wall  boundary  layer  gradually  transforms  that  layer  into  a  boundary 
layer  with  a  new  equilibrium  state  governed  by  the  large-scale  sweeps 
and  ejections.  The  behavior  of  the  mean  velocity  profiles  in  the  Clauser 
plots  can  be  explained  in  terms  of  this  picture  of  a  developing  flow. 

The  boundary  layer  growing  on  the  flat  wall  of  the  inlet  accumulates 
a  relatively  large  amount  of  momentum-deficient  fluid  near  the  wall  by 
the  time  it  reaches  the  curved  section.  As  the  flow  enters  the  curve, 
this  low-speed  fluid  is  ejected  away  from  the  wall  by  the  large-scale 
sweeps.  The  incoming  sweeps  begin  to  establish  a  new  equilibrium  condi¬ 
tion  in  which  the  ejection  of  low-speed  wall  fluid  is  limited  by  the 
rate  at  which  fluid  is  decelerated  near  the  wall  by  viscosity.  Since  the 
sweeps  replace  fluid  near  the  wall  with  higher  speed  fluid  from  farther 
out.  the  mean  skin  friction  increases.  Meanwhile,  the  initial,  large 
body  of  low-speed  wall  fluid  migrates  outward,  away  from  the  wall,  grad¬ 
ually  mixing  with  the  surrounding  fluid.  In  mean  velocity  profiles  this 
mixing  region  is  believed  to  correspond  to  a  quasi-logarithmic  zone  with 
a  slope  less  than  that  corresponding  to  the  Law  of  the  Mall.  Fig.  37 
shows  the  streamwise  development  of  the  mean  velocity  profiles  with  the 
quasi-log  zones  marked  by  line  segments.  (In  the  profile  at  75°.  a  new 
logarithmic  region  may  actually  be  present  close  to  the  wall.)  Conse¬ 
quently.  all  conclusions  from  the  present  study  may  be  said  to  represent 
only  developing  concave  wall  flow.  In  practice.  regions  of  extended 
curvature  seldom  occur.  Almost  all  real  flows  are  therefore  developing, 
and  the  present  results  and  this  flow  model  are  of  practical  signifi¬ 
cance. 
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favorable  pressure  gradient  since  separation  on  the  opposing  wall  nulli¬ 
fies  the  effects  of  the  channel  expansion  to  some  extent.  The  result  is 
a  new,  thinner  boundary  layer  growing  on  the  flat,  recovery  region  test 
wall.  This  layer  has  large-scale  motions  in  the  outer  flow  whose  inertia 
makes  them  slow  to  decay  as  they  are  convected  downstream  from  the 
curved  section.  The  large  pockets  of  laminar-like  fluid  seen  near  the 
wall  may  indeed  be  the  results  of  local  relaminarization  in  various 
parts  of  the  boundary  layer.  The  laminar-like  pockets  in  the  outer  flow 
are  probably  chunks  of  the  potential  core  which  were  entrained  by  the 
large-scale  motions. 

The  flow  model  for  the  concave  wall  boundary  layer  points  to  large 
changes  in  the  distributions  of  the  sizes  of  the  eddies  which  contribute 
to  the  turbulence  as  the  flow  goes  from  flat  to  concave  wall  conditions. 
The  basic  concave  instability  amplifies  the  outer  layer  scale  of  motion 
to  create  the  sweeps  which  then  interact  with  the  wall  layers  and  change 
their  structure  by  pushing  them  into  the  outward  ejections.  Simple 
measurements  of  u'  cannot  reveal  these  changes.  Although  fiow  visualiza¬ 
tion  gives  a  good  qualitative  picture,  quantitative  data  of  a  much  more 
detailed  kind  are  needed.  Space-time  correlations  of  several  velocity 
components  would  be  particularly  useful  in  determining  the  relative  con¬ 
tributions  of  the  large  and  small  scales  and  how  the  distribution  of 
these  scales  of  motion  changes  as  we  go  from  a  flat  wall  to  a  concave 
wall,  and  back  again.  None  of  the  studies  obtained  to  date,  including 
our  own,  satisfy  this  requirement. 
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Chapter  V 


CONCLUSIONS 

1.  The  results  of  the  present  experiment  provide  a  much  clearer  pic¬ 
ture  of  the  nature  of  the  flow  in  turbulent  boundary  layers  on 
concave  ualls.  This  picture  aids  understanding  of  past  experi¬ 
ments  and  should  help  guide  the  course  of  future  studies. 

The  following  description  of  a  turbulent  boundary  layer  on  a 
concave  wall  is  proposed: 

a)  The  structure  of  the  boundary  layer  appears  to  be  dominated  by 
large-scale  sweeps  and  ejections.  These  sweeps  and  ejections 
are  coherent  over  the  height  of  the  boundary  layer  normal  to 
the  wall  and  are  typically  about  half  the  boundary  layer 
thickness  in  span.  The  streamwise  extent  of  the  structures  is 
of  the  order  of  a  few  boundary  layer  thicknesses.  The  struc¬ 
tures  are  distributed  randomly  over  the  wall  and  essentially 
travel  with  the  mean  flow  with  the  sweeps  moving  faster  than 
the  ejections  since  they  originate  farther  from  the  wall.  The 
spanwise  locations  of  the  structures  are  not  stationary  in  our 
apparatus*’0  nor  do  the  structures  give  rise  to  well-defined 
vortices  or  roll  cells. 


10  Some  past  experiments  by  others  did  find  stationary  structures  per¬ 
haps  due  to  apparatus  irregularities.  The  question  of  what  determines 
fixed  versus  moving  structures  is  not  settled.  See  below. 
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The  effect  of  the  sweeps  is  to  bring  high-speed  fluid  from 


the  outer  flow  closer  to  the  wall.  This  action  displaces  the 
low-speed  fluid  near  the  wall  which  is  ejected  away  from  the 
wall.  These  large  movements  of  fluid  appear  to  inhibit  the 
near-wall  bursting  process,  the  major  turbulence  production 
mechanism  for  flat  wall  flow.  However,  the  interaction  between 
the  ejections  and  sweeps  appears  to  create  a  Helmholtz  insta¬ 
bility  of  the  shear  layer  between  the  high-speed  sweeps  and 
the  low-speed  ejections,  and  creates  additional  turbulence 
production  throughout  the  rest  of  the  boundary  layer. 

b)  An  effect  of  the  large-scale  sweeps  and  ejections  is  a  fuller 
mean  velocity  profile  which  gradually  develops  through  the 
curved  region  as  the  flow  moves  toward  a  new  state  of  equili¬ 
brium.  The  usual  Law  of  the  Mall  is  either  attenuated  or  elim¬ 
inated  by  the  onset  of  concave  curvature;  instead,  a  quasi-lo- 
garithmic  zone  appears  which  seems  to  correspond  to  the 
outward  migration  of  the  ejections  created  at  the  start  of 
curvature. 

c)  The  large-scale  fluid  motions  induced  in  the  curved  boundary 
layer  persist  long  distances  after  wall  curvature  has  ended. 
The  removal  of  low-speed  fluid  from  the  near-wall  region  in 
the  curve  causes  the  boundary  layer  in  the  recovery  region  to 
appear  somewhat  like  a  newly  developing,  laminar-like  boundary 
layer  with  little  near-wall  bursting. 
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The  present  results  lend  credibility  to  those  earlier  investigations 
which  failed  to  find  evidence  of  large-scale  structures  in  their 
measurements  of  mean  quantities.  This  suggests  the  need  for  further 
work  to  determine  under  uliat  circumstances  the  sweeps  and  ejections 
remain  fixed  in  position.  In  particular,  weak  upstream  disturbances, 
such  as  those  introduced  by  screens,  honeycombs,  and  trips  may  have  con¬ 
taminated  several  of  the  more  recent  studies,  causing  their  large-scale 
structures  to  occur  at  fixed  locations. 

2.  The  present  facility  provides  an  acceptably  uniform  flow  which 
can  be  adequately  visualized.  The  flexibility  inherent  in  the 
water  channel's  design  should  make  the  channel  a  useful  tool  for 
future  investigations  into  a  variety  of  flows.  Much  of  this  flex¬ 
ibility  lies  in  the  channel's  ability  to  accept  different  test 
wall  inserts.  The  large  size  of  the  outer  channel  allows  the  use 
of  a  wide  variety  of  inner  channel  geometries.  Another  useful 
feature  is  the  ability  to  produce  laminar  or  transitional  flows 
by  varying  the  flow  rate.  More  drastic  modifications  are  also 
possible:  the  outer  channel  modules  can  be  unbolted  and  recon¬ 
nected  in  different  configurations.  Furthermore,  should  any  of 
the  modules  become  worn  or  damaged,  they  can  be  replaced  rela¬ 
tively  easily. 

3.  The  computer  controlled  hot  film  probe  enabled  us  to  acquire 
velocity  data  with  good  repeatability  in  what  would  otherwise  be 
an  unfavorable  environment  for  such  an  instrument.  The  largest 
uncertainty  in  the  results  was  due  to  calibration  uncertainty; 
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further  refinements  in  the  speed  sensor  of  the  instrument  plat¬ 
form  could  reduce  this  uncertainty. 


Chapter  V! 


SUGGESTIONS  FOR  FUTURE  RESEARCH 

Future  investigations  of  turbulent  boundary  layers  on  concave  walls 
might  profitably  proceed  along  various  paths*  two  of  which  are: 

1.  Obtain  more  detailed  quantitative  data  on  the  flow.  Suggestions 
include: 

a)  Measurements  of  turbulent  quantities*  especially  turbulent 
shear  stress,  obtained  by  conditional  sampling  employing 
simultaneous  visualization  to  identify  structural  features. 

b)  Measurements  of  space-time  correlations  to  quantify  the  size 
of  the  large-scales  of  motion. 

c)  Quantitative  analysis  of  hydrogen  bubble  pictures  by  image 
processing  to  more  accurately  determine  the  shapes,  sizes,  and 
motions  of  the  large-scale  structures. 

2.  Attempt  to  resolve  the  role  of  upstream  or  other  disturbances  in 
determining  the  flow  in  the  curve.  A  systematic  study  of  the 
effects  of  various  commonly  encountered  wind  tunnel  or  water 
channel  flow  abnormalities  (e.g.  streamwise  vortices  or  wakes 
from  the  inlet  screens)  can  perhaps  answer  the  question  of  why 
some  of  the  flows  contain  stationary  structures  while  others  do 
not. 
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Appendix  A 


MATER  CHANNEL  DESIGN 


A. 1  OUTER  CHANNEL  DESIGN 

The  basic  design  criteria  for  the  outer  channel  structure  were: 

1.  that  it  safely  hold  the  large  volumes  of  water  required, 

2.  that  it  be  as  large  as  practically  and  economically  possible, 

3.  that  it  permit  visual  access  to  as  much  of  the  test  wall  as  pos¬ 
sible,  and 

4.  that  it  be  flexible  enough  in  design  to  allow  for  future  modifi¬ 
cations  or  repairs. 

The  last  criterion  and  the  desire  to  facilitate  the  construction  led  to 
a  modular  approach  where  the  channel  could  be  assembled  from  parts  of  an 
easily  manageable  size.  Consequently,  the  fabrication  of  the  parts  could 
be  done  elsewhere  and  the  channel  assembled  on  site  with  little  or  no 
additional  modification  to  the  parts. 

Visualization  requirements  dictated  that  most  of  the  channel  be 
transparent.  To  support  the  weight  and  pressure  of  the  large  volumes  of 
water,  a  steel  reinforcing  frame  was  necessary.  In  order  to  facilitate 
fabrication,  the  channels  and  frames  were  designed  as  separate  units. 
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Adjustable  braces  (see  arrows)  were  used  to  fit  each  frame  to  Its  chan¬ 
nel*  as  shown  in  Fig.  39.  In  addition*  this  design  feature  increased  the 
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modularity  of  the  apparatus.  Me  chose  plexiglas  over  glass  for  the  chan¬ 
nels  because  of  the  relative  ease  with  which  it  can  be  machined  and 
joined.  This  was  considered  more  important  than  the  scratch  resistance 
of  the  material.  The  standard  4'x#'  sheets  plexiglas  comes  in,  led  to  a 
channel  module  size  of  4'  long.  A'  high,  and  2'  wide.  The  size  and 
weight  of  these  modules  permitted  easy  handling  by  two  people  without 
using  power  equipment. 

Each  channel  module  had  flanges  at  either  end  which  were  used  to  bolt 
adjacent  channels  to  each  other.  A  soft  neoprene  rod,  sandwiched 
between  flange  faces  and  held  in  place  by  a  groove  in  one  face,  sealed 
each  joint.  Design  data  for  the  groove  was  obtained  from  0-ring  manu¬ 
als,  but  the  best  bolt  spacing  was  ultimately  determined  by  constructing 
a  small,  prototype  flange. 

The  inlet  and  outlet  plenums  were  fabricated  out  of  sheet  steel  rein¬ 
forced  with  channel  bar  (small  section  channel)  because  we  did  not 
expect  to  have  to  visualize  the  flow  in  these  parts.  A  special  neoprene- 
asphalt  paint  was  painted  on  the  inside  surfaces  to  prevent  rusting. 
This  paint  has  the  advantages  of  being  very  flexible  and  easy  to  apply. 
In  addition  it  possesses  outstanding  corrosion  resistance  properties. 
The  plenum  exteriors  and  the  frames  were  painted  with  an  industrial 
grade  waterproof  epoxy  paint. 

The  joints  between  the  channels  and  the  plenums  were  fabricated  on 
site  since  it  was  difficult  within  cost  constraints  to  construct  flanges 
of  sufficient  flatness  on  the  plenums.  At  first,  closed  cell  neoprene 
foam  rubber  was  used  to  bridge  the  gap  between  the  plexiglas  and  the 


100  - 


steel,  but  this  material  was  found  to  shrink  considerably  when  exposed 


to  our  water  supply,  due  perhaps  to  the  high  chlorine  levels.  Finally, 
solid  neoprene  sheet  was  employed  with  better  success.  Contact  cement 
was  used  to  bond  the  neoprene  to  the  plexiglas  and  steel,  followed  by 
silicone  rubber  caulk  to  make  the  joints  water  tight. 

The  stress/strain  calculations  needed  to  specify  the  sizes  of  the 
various  members  of  the  frames  and  channels  were  done  by  breaking  the 
structures  into  their  individual  components  and  applying  two-dimensional 
strength  of  materials  formulae  for  beam  deflections,  assuming  that  each 
member  was  independent  of  the  rest.  In  most  cases  the  results  were  con¬ 
servative  estimates  of  the  deflection  to  be  expected.  Most  of  the  beams 
were  sized  to  keep  deflections  under  0.050".  In  the  actual  structures, 
three-dimensionality  would  make  the  structures  stronger  and  the  deflec¬ 
tions  less.  Beam  sections  (i.e.  angle,  channel,  or  I-beam)  were  selected 
to  minimize  weight  and  cost  and  to  simplify  the  design  and  fabrication 
of  the  beam  joints. 

The  size  and/or  quantity  of  the  parts  made  fabrication  of  most  of 
them  in  our  own  shop  unfeasible  so  vendors  were  contracted  to  do  the 
initial  construction.  The  steel  frames  were  assembled  by  Bill  Arbogast 
of  Precision  Melding,  Mountain  View,  who  also  painted  the  parts.  The 
outer  channel  walls  were  fabricated  by  Star  Products  in  Sunnyvale,  while 
the  inner  channel  walls  (described  in  more  detail  below)  were  con¬ 
structed  by  the  Hughes  Plastio  Co.  In  Redwood  City.  Me  received  good, 
friendly  service  from  these  vendors  and  recommend  them  to  future  rig 
builders. 
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The  costs  of  the  various  parts  were  roughly  $13,000  for  the  steel 
frames  and  plenums,  $12,000  for  the  outer  channel  walls,  and  $5,000  for 
the  inner  channel  walls. 

A. 2  JBlSHtM  m  SPEC.iriJATIQJiS  fll  JLLQJ4  C9K.TRPL  ELEBEHI.S 
The  various  flow  control  elements  in  the  water  channel  were  sized  assum¬ 
ing  a  nominal  flow  speed  of  0.5  ft/sec  through  the  test  section,  whose 
cross-sectional  area  was  assumed  to  be  3  to  4  ft2.  However,  the  pump 
obtained  is  capable  of  delivering  1200  gal/min.  Me  urre  unable  to  find 
much  literature  with  design  information  (mainly  pressure  drop  coeffi¬ 
cients)  for  low-speed  flows  such  as  ours.  However  Miller's  Internal  Flow 
Systems1 1  served  as  an  excellent  guide  to  design  considerations  for 
higher  speed  flows,  and  we  were  able  to  extrapolate  from  those  figures 
and  obtain  reasonable  guesses  for  our  system.  The  reader  should  refer 
back  to  Fig.  5  for  a  diagram  of  the  channel  and  its  flow  control  ele¬ 
ments. 

A. 2.1  SiflJLgja  flmmUjttM 

The  overflow  in  the  inlet  plenum  was  designed  to  help  regulate  the  water 
level  in  the  channel  as  well  as  serve  as  a  safety  feature  to  prevent  the 
water  from  overflowing  the  channel  walls  should  something  obstruct  the 
flow  downstream.  Me  required  that  the  overflow  pipe  pass  at  least  a 
quarter  of  the  maximum  flow  of  1200  gal/min.  A  A"  diameter  pipe  with  a 
A'  vertical  drop,  capable  of  passing  about  600  gal/min,  was  considered 
adequate. 


"  Miller,  Donald  S.  (1960)  Internal  Flow  Systems.  Vol.  5  in  the  British 
Hydromechanics  Research  Assoc.  Fluid  engineering  Series. 


A. 2. 2  Perforated  Pipe 

In  order  to  reduce  the  effects  of  pump  surging  and  to  break  up  the  large 
eddies  generated  by  the  valves  and  elbous  in  the  delivery  piping,  ue 
forced  the  flow  to  pass  through  a  perforated  pipe  uith  a  large  pressure 
drop  across  its  surface.  The  pipe  uas  simply  a  vertical,  detachable 
extension  of  the  delivery  pipe  coming  through  the  bottom  of  the  inlet 
plenum.  Although  some  simple  calculations  using  Bernoulli's  equation 
provided  a  rough  estimate  of  the  porosity  of  the  pipe  required  for.  say. 
a  head  loss  of  one  foot  of  water,  the  final  design  was  determined  by 
trial.  The  final  design  incorporates  9  columns  of  0.5"  diameter  holes 
placed  on  the  side  of  the  pipe  facing  upstream,  future  workers  may  wish 
to  increase  the  number  of  holes  to  increase  the  channel's  flowrate, 
since  this  pressure  drop  is  a  major  resistance  in  the  flow  system. 

a. 2. 3  perforated  Plate 

The  requirement  for  the  perforated  plate  was  that  it  cause  a  loss  of 
about  100  velocity  heads  and  have  a  porosity  less  than  20%  to  prevent 
jet  coalescing.  For  plate  porosities  on  the  order  of  10%.  such  pressure 
drops  occur  for  a  wide  range  of  pitch  Reynolds  numbers  (Reynolds  numbers 
based  on  the  mean  velocity  upstream  of  the  plate  and  the  distance 
between  holes).  The  final  design.  Influenced  to  some  extent  by  the 
availability  of  stainless  steel  perforated  plates,  uses  3/16"  diameter 
holes  spaced  0.5"  apart  in  a  hexagonal  pattern,  with  an  open  area  of 
13%.  The  plate  was  punched  for  us  by  the  Duus  Perforating  Co.  in  San 
Jose. 
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A. 2. 4  HflUfeVLEfimh 


Pressure  losses  across  honeycomb  flow  straighteners  are  usually  neglige-  ^ 

able.  Me  selected  a  honeycomb  panel  (Kexeet  HRP  glass-reinforced  phe-  j 

nolic)  3"  thick  Kith  1/4"  cells.  the  cell  diameter  Reynolds  number  is 

then  about  300.  | 

i 

A. 2. 5  Screens  j 

i 

The  screens  were  selected  to  have  the  lowest  possible  wire  diameter  key-  ' 

nolds  number  and  an  open  area  of  at  least  60%.  Stainless  steel  wire 

I 

] 

cloths  available  from  Tyler  Industrial  Products,  Hayward,  a  large  sup- 

t 

i  plier  of  industrial  wire  cloths  and  screens>  use  wires  with  diameters  as 

>  small  as  0.004".  Using  wire  of  this  size  in  a  screen  with  60%  open  area 

(56  mesh)  at  the  nominal  flowrate,  the  wire  diameter  Reynolds  number  is 
less  than  10.  Me  estimated  the  pressure  loss  across  each  screen  to  be 
about  3  velocity  heads  or  about  0.015".  The  screens  were  mounted  in 
plexiglas  frames  in  a  manner  much  like  conventional  window  screens.  A 
flexible  rubber  rod  was  squeezed  into  a  groove  along  with  the  screen, 
holding  it  in  place.  Me  were  able  to  pre-tension  the  frames  so  that  the 
screens  remained  fairly  flat  and  taut  when  placed  in  the  flow. 

a. 2. 6  ttaazlt 

Following  the  advice  of  Bradshaw  £  Pankhurst12  we  shaped  our  nozzle 
using  cubic  splines.  The  geometrical  requirements  for  the  nozzle  are 

shown  in  Fig.  39.  where  the  solid  lines  represent  the  walls  of  the  i 


12  Bradshaw.  P.  £  Pankhurst.  R.  C.  (1964)  "The  design  of  low-speed  wind 
tunnels."  Progress  in  Aeronautical  Scienoes.  Vol.  5.  Pergamon  Press. 
I  PP-  1-69. 
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plenum,  outer  channel,  and  test  walls,  and  the  dotted  lines  represent 
the  approximate  shape  of  the  nozzle  to  be  designed.  Because  of  the 
desire  to  allow  more  space  between  the  test  wall  and  the  outer  channel, 
the  nozzle  had  to  be  asymetrical. 

The  matching  point  for  each  pair  of  cubic  arcs  was  chosen  to  lie  on 
the  line  joining  the  two  endpoints,  about  2/3  of  the  distance  downstream 
of  the  inlet  (refer  to  Fig.  40).  Instead  of  performing  two  calculations, 
one  for  a  width  of  B.5"  and  another  for  a  width  of  13. 5",  one  calcula¬ 
tion  was  done  for  a  width  of  11.0"  and  the  results  scaled  accordingly. 
The  6  constants  in  the  two  cubic  arcs: 
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y i=aj+bjX+ciX2+d jx3»  1=1,2  (A.1) 

were  determined  by  imposing  the  following  boundary  conditions: 

1.  The  ends  of  the  arcs  must  be  fixed  at  their  proper  locations: 

y,(0)=11,  yi(21)=y2(21)=3,  y2C30>=0  (A. 2) 

2.  The  first  derivatives  must  vanish  at  the  inlet  and  outlet  points: 

dy!  dy2 

- <0)  = - ( 30 ) =0  (A. 3) 

dx  dx 

3.  The  slopes  must  match  at  the  matching  point: 

dy*  dy2 

- (21)= - ( 2 1 ) =m  (A. 4) 

dx  dx 


4.  The  second  derivative  at  the  outlet  must  vanish: 
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dy22 

- (30)=0 

dx2 


(A.S) 


The  solution  to  the  above  equations  and  boundary  conditions  is: 

pi. 00  -  6.S00x10‘5x2  -  5. 396xl0**x*  C01xi21) 

y=|  (A. 6) 

U11.1  +  ll.llx  +  0.3703X2  -  4.  1 1Sx10*3x3  C2Hxi30) 

The  y-values  from  these  equations  were  multiplied  by  A. 5/11  and  13.5/11 

to  obtain  the  final  nozzle  configuration. 

A. 2. 7  Drain  Val ve 

The  initial  design  of  the  channel  incorporated  an  adjustable  weir  at  the 
end  of  the  channel  to  regulate  the  level  of  the  uater  and  its  flourate. 
However,  the  resulting  waterfall  entrained  unacceptably  large  amounts  of 
air  and  this  design  was  scrapped  when  we  wanted  to  start  using  hot  film 
probes  and  hydrogen  bubbles.  The  problem  was  to  regulate  the  flow  going 
back  into  the  sump  without  entraining  any  air.  There  was  no  room  to 
place  conventional  valves  in  the  piping  that  was  above  floor  level,  and 
valves  placed  in  the  sump  would  be  difficult  to  get  to.  Me  tried  using 
homemade  disc  valves  which  were  controlled  by  cables  running  back  out  of 
the  plenum,  but  the  hydrodynamic  forces  were  too  large  to  make  fine 
adjustments  possible.  Finally  a  special  valve,  shown  in  Fig.  41,  was 
designed  and  fabricated  .  Once  the  channel  was  filled,  the  depth  of  the 
water  in  the  plenum  prevented  strong  vortices  from  forming  and  any  air 
remaining  in  the  pipe  below  the  valve  was  quickly  entrained  and  removed 
by  the  flow.  Although  the  valve  had  a  tendency  to  vibrate  some.  It  was 
rugged  enough  to  withstand  the  large  hydrodynamic  forces  and  performed 
satisfactorily. 
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A. 3  INNER  CHANNEL  DESIGN 


In  order  to  minimize  the  effects  of  the  stresmuise  pressure  gradient  on 
the  turbulent  structure  at  the  test  wall,  a  design  calculation  was  per¬ 
formed  to  determine  the  geometry  of  the  opposite  uall .  First  the  poten¬ 
tial  flou  velocities  at  the  tuo  walls  (test  wall  with  concave  surface 
and  opposite*  flow  control  wall)  were  computed  for  a  variety  of  opposite 
wall  geometries.  The  potential  uall  velocities  were  then  used  in  a 
flat-wall  boundary  layer  calculation  to  compute  the  streamuise  develop¬ 
ment  of  the  displacement  thicknesses  on  both  walls.  The  sum  of  the  dis¬ 
placement  thicknesses  at  each  streamuise  location  was  then  added  to  the 
local  channel  width  to  determine  the  shape  of  the  opposite  uall.  A  more 
detailed  description  follows. 

The  geometry  of  the  inner  walls  was  to  a  large  extent  determined  by 
the  size  constraints  imposed  by  the  outer  structure,  which  in  turn  was 
made  as  large  as  economically  and  practically  feasible.  The  radius  of 
the  test  uall  was  made  to  fit  the  corner  module  with  a  minimum  clearance 
between  the  inner  and  outer  uall  of  about  9"  for  access  to  the  dye  slots 
or  other  devices.  The  opposite  wall  was  then  initially  placed  so  that 
there  was  a  clearance  of  about  4"  between  it  and  the  outer  wall.  This 
then  determined  the  inlet  geometry  and  a  nozzle  was  designed  to  carry 
the  flou  from  the  settling  chamber  to  the  channel. 

In  order  to  calculate  the  potential  flow  through  the  channel,  a  fast 
Laplace  solver  developed  by  Wooley,  White,  Bardina,  C  Kline13  was  used 


13  A  description  of  the  program  may  be  found  in  "A  Boundary  Method  for 
Potential  Flows"  by  Juan  G.  Bardina.  Report  IL-17,  Thermosciences 
Div.,  Dept,  of  Mechanical  Engineering,  Stanford  Univ.,  1977. 
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which  determines  the  tangential  velocities  at  specified  points  along  the 
boundaries  given  the  normal  velocities.  The  program*  PHIVN,  calls  on 
several  large  library  subroutines  and  is  not  listed  here.  In  this  case 
the  coordinates  of  points  describing  the  channel  geometry  were  entered 
in  a  counterclockwise  fashion  along  with  the  normal  velocity  at  the  mid¬ 
dle  of  the  next  segment.  This  normal  velocity  was  zero  for  all  points 
except  two  at  the  inlet  and  outlet.  The  initial  geometry  was  a  90° 
annular  sector*  outer  radius  136  cm  (53.6")*  width  24.1  cm  (9.5'/)>  with 
straight  lead-in  and  trailing  sections.  Because  of  the  symmetry  of  the 
channel*  only  half  of  the  coordinates  had  to  be  entered*  the  rest  were 
generated  by  the  computer.  The  coordinates  were  stored  in  separate  files 
for  the  inner  and  outer  walls.  Small  auxiliary  programs  combined  the 
files  to  create  an  input  for  PHIVN.  By  running  PHIVN  and  gradually 
adjusting  the  shape  of  the  flow  control  wall*  the  velocity  along  the 
entire  test  wall  was  made  uniform  to  within  3%  of  the  mean.  Finer 
adjustments  were  not  undertaken  because  it  was  known  that  the  walls 
could  not  be  manufactured  to  closer  tolerances.  The  final  potential 
flow  geometry  is  described  by  the  coordinates  in  Table  4.  Note  that  X 
and  Y  refer  to  standard*  cartesian  coordinates.  Fig.  42  shows  these 
points  plotted  along  with  a  circular  arc*  which  is  shown  by  the  dashed 
line.  Fig.  43  graphs  the  computed  wall  velocities  corresponding  to  this 
channel  geometry. 

The  output  of  PHIVN  was  then  modified  by  another  small  program  to 
serve  as  the  input  to  a  flatplate  boundary  layer  code.  The  data  for  each 
wall  was  isolated*  converted  into  the  form  of  streamwise  displacement 
versus  velocity*  and  appended  to  additional  points  needed  to  make  up  the 
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TABLE  4 


Channel  configuration,  cartesian  coordinates. 


Concave  (test)  wall 

X,  cm  Y,  cm 


-0. 1219E+03 

-0. 1361E+03 

-0. 1524E+02 

-0. 1361E+03 

-0. 1016E+02 

-0. 1361E+03 

-0 . 5080E+0 1 

-0. 1361E+03 

0 . 0000E+00 

-0. 1361E+03 

0 . 50  80E+0 1 

-0. 1360E+03 

0. 1016E+02 

-0. 1 358E+03 

0. 1524E+02 

-0. 1353E+03 

0. 2032E+02 

-0. 1346E+03 

0 . 2540E+02 

-0. 1338E+03 

0.3302E+02 

-0. 1321E+03 

0. 4572E+02 

-0.  1 282E+03 

0 . 6096E+02 

-0. 1217E+03 

0 . 7620E+02 

-0. 1128E+03 

0 . 9627E+02 

-0.9627E+02 

Convex  (control)  wall 

X#  cm  Y ,  cm 

0.11 20E+03 

0.  1219E+03 

0. 1120E+03 

0. 1524E+02 

0.1 121E+03 

0.  1270E+02 

0. 1121E+03 

0. 1016E+02 

0. 1122E+03 

0. 7620E+01 

0.  1124E+03 

0. 5080E+01 

0. 1 125E+03 

0. 2540E+01 

0.  1127E+03 

0 . 0000E+00 

0. 1130E+03 

-0 . 2540E+0 1 

0. 1132E+03 

-0. S080E+01 

0. 1133E+03 

-0.7620E+01 

0. 1133E+03 

-0. 1016E+02 

0. 1133E+03 

-0. 1270E+02 

0. 1130E+03 

-0. 1524E+02 

0. 1126E+03 

-0. 1778E+02 

0. 1122E+03 

-0. 2032E+02 

0.11 12E+03 

-0 . 2540E+02 

0.  1092E+03 

-0.3175E+02 

0. 1075E+03 

-0.3810E+02 

0. 1033E+03 

-0.4826E+02 

0.9639E+02 

-0.6096E+02 

0.891 5E+02 

-0.7112E+02 

0 . 8064E+02 

-0 . 8064E+02 

X,  cm  Y,  cm 


0 . 1 1 28E+03 

-0 . 7620E+02 

0.  12 17E+03 

-0 . 6096E+02 

0.  1282E+03 

-0.4572E+02 

0. 1321E+03 

-0.3302E+02 

0. 1338E+03 

-0. 2540E+02 

0. 1346E+03 

-0 . 2032E+02 

0. 1 353E+03 

-0. 1 524E+02 

0. 1 358E+03 

-0. 1016E+02 

0. 1360E+03 

-0. 5080E+01 

0. 1361E+03 

0 . 0000E+00 

0. 1361E+03 

0 . 5080E+0 1 

0. 1361E+03 

0.  1016E+02 

0. 1361E+03 

0. 1524E+02 

0. 1361E+03 

0. 1219E+03 

X,  cm 

Y ,  cm 

0.71 12E+02 

-0.8915E+02 

0. 6096E+02 

-0. 9639E+02 

0 . 4826E+02 

-0. 1033E+03 

0.3810E+02 

-0. 1075E+03 

0. 3175E+02 

-0. 1092E+03 

0 . 2540E+02 

-0. 1 1 12E+03 

0.2032E+02 

-0. 1122E+03 

0. 1778E+02 

-0. 1 126E+03 

0. 1524E+02 

-0.11 30E+03 

0.  1270E+02 

-0. 1133E+03 

0. 1016E+02 

-0. 1133E+03 

0 . 7620E+0 1 

-0. 1133E+03 

0 . 5080E+0 1 

-0. 1132E+03 

0 . 2540E+0 1 

-0. 1 130E+03 

0 . 0000E+00 

-0. 1127E+03 

0 . 2540E+0 1 

-0. 1125E+03 

0 . 5080E+0 1 

-0 . 1 124E+03 

0.7620E+01 

-0. 1 122E+03 

0. 1016E+02 

-0. 1121E+03 

0. 1270E+02 

-0. 1 121E+03 

0. 1524E+02 

-0. 1120E+03 

0. 1219E+Q3 

-0. 1 120E+03 
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entire  wall.  A  code  written  by  Cebeci  £  Bradshaw'*  used  the  velocity 
data  to  compute  boundary  layer  parameters.  The  resulting  displacement 
thioknesses  are  shown  in  Fig.  44.  For  the  case  of  the  variable  width 
channel*  separation  occurs  when  the  flow  is  expanded  downstream  of  the 
bond.  The  sum  of  the  displacement  thicknesses  for  the  two  walls  at  a 


Cebeci,  Tuncer  C  Bradshaw,  Peter  (1B77)  Momentum  Transfer  ifl  ifiuarimt 
Lavers.  Hemisphere  Publishing  Corp.  fi  McGraw-Hill  Book  Co.  pp. 
258-26*. 


(cm/s)  LU  (cm/s) 


POTENTIAL  WALL  VELOCITIES 


X  (cm) 


POTENTIAL  WALL  VELOCITIES 


X  (cm) 


Figure  43«  Computed  potential  well  velocities,  oonetent  end  variable 
width  channels. 


DELTA*  (cm)  DELTA*  (cm) 


given  x,  or  angle  (in  the  bend),  uas  then  added  to  the  uidth  of  the 
channel  to  obtain  the  final  opposite  uall  geometry.  This  correction  uas 
at  most  about  55C  of  the  channel  uidth.  The  channel  uidth  including  the 
displacement  thickness  correction  is  tabulated  in  Table  5.  For  the  part 
of  the  channel  dounstream  of  the  calculated  separation  point,  no  addi¬ 
tional  corrections  uere  made. 


Channel  uidth  as 

TABLE 

a  function 

5 

of  streamuise  distance. 

x.  cm 

u,  cm 

x,  cm 

u,  cm 

122 

24.92 

533 

24.54 

244 

25.66 

559 

24.67 

366 

26.32 

584 

24.68 

419 

26.61 

610 

24.82 

490 

25.52 

635 

25.32 

508 

24.21 

660 

25.36 

Admittedly,  a  flatplate  code  uas  used  for  curved  flou;  houever,  ue 
felt  that  for  a  first  approximation,  the  results  uould  greatly  reduce 
the  streamuise  pressure  gradient  on  the  test  uall  and  uould  certainly  be 
an  improvement  over  no  corrections  at  all. 

The  finished  ualls  uere  anchored  to  the  channel  bottom  by  adjustable 
brass  braces  and  to  the  top  of  the  channel  frame  by  threaded  rods  as 
shoun  in  Fig.  45.  This  mounting  system  enabled  us  to  position  the  ualls 
very  accurately.  The  various  sections  uere  attached  to  each  other  by 
screus. 
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a)  top  and  (b)  bottom  of  the  Inner 


Appendix  B 


TRAVERSE  DESIGN 


B.  1  HARDWARE 

In  order  to  accurately  position  probes  uithin  the  channel  a  tuo-axis 
traverse  had  to  be  designed  which  could  move  probes  about  in  the  channel 
cross-plane  (y-z).  The  large  size  of  the  rig  and  the  potential  availa¬ 
bility  of  a  computer  led  to  a  design  incorporating  digitally  controlled 
stepper  motors.  The  motors  chosen  for  the  job  were  capable  of  rotating 
their  shafts  in  0.9°  increments  with  an  accuracy  of  i5%.  Since  position¬ 
ing  on  the  horizontal  (y)  axis  had  to  be  mere  accurate  than  on  the  ver¬ 
tical  (z)  axis,  a  lead-screw  was  used  on  the  former  and  a  rack  and  pin¬ 
ion  on  the  latter.  The  wide  range  of  gearing  available  made  it  possible 
to  convert  the  rotation  of  the  motor  shafts  into  precise  metric  steps 
while  still  using  American  Standard  gears,  screw  threads,  etc.  In  the 
final  design.  20  full  (1.8°)  steps  of  the  horizontal  motor  caused  the 
traverse  to  move  0.1  mm  with  an  error  of  0.05%.  and  4  steps  of  the  ver¬ 
tical  motor  resulted  in  a  1.0  mm  vertical  displacement  with  an  error  of 
0.25%.  The  parts  specifications  are  shown  in  Table  6.  Initially,  the 
vertical  portion  of  the  traverse  rode  directly  on  the  lead-screw.  How¬ 
ever,  eccentricity  of  the  screw  threads  and  the  shaft  ends  resulted  in 
objectionable  oscillations  of  the  vertical  axis,  so  linear  bearings  were 
added  to  take  the  vertical  load  off  of  the  lead-screw.  A  close-up  photo 
of  the  traverse  appears  in  Fig.  AH. 
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TABLE  6 


Traverse  Parts  Specifications. 


Horizontal  Mi • 


Gear 

Teeth 

Pi  tch 

Tooth  Angle 

Pitch  Oia. 

drive 

17 

24 

14.5° 

0.708" 

driven 

27 

24 

14.5® 

1.125" 

Center 

distance 

=0.917". 

Lead  screw  thread:  3/8-16  stainless  steel  (specially  rolled). 


Vertical  Ah  1.8  = 

Gear  Teeth  Pitch  Tooth  Angle  Pitch  Oia. 
pinion  20  32  14.5®  0.625" 

rack  —  32  14.5° 


B. 2  MOTOR  CONTROLS 

In  order  for  a  stepper  motor  to  turn,  its  4  leads  must  be  energized  in  a 
particular  sequence.  To  do  this  manually  would  be  tedious  indeed*  so 
various  electronic  controllers  are  available  to  simplify  the  task. 
Because  we  wanted  to  be  able  to  position  the  probe  without  necessarily 
relying  on  a  computer*  we  purchased  what  are  known  as  preset  indexer 
modules  (PIM's).  These  devices  send  out  the  sequence  of  pulses  necessary 
to  move  the  motor  a  user-specified  number  of  steps.  The  indexers  were 
available  as  unmounted  po  boards;  To  save  space  and  money,  both  indexers 
were  mounted  in  a  common  ohassis  along  with  the  necessary  power  supplies 
and  switches.  The  various  leads  controlling  the  number  of  steps*  the 
initiation  of  a  move,  ete.#  while  being  wired  internally  to  the  switches 


on  the  front  panel,  were  also  made  accessible  through  multi-pin  connec¬ 
tors  on  the  back  panel  for  optional  computer  control. 

B. 3  COMPUTER  CONTROLS 

Shortly  after  the  motor  controller  uas  built,  we  were  able  to  obtain  a 
DEC  MINC-11  microcomputer  for  our  use.  To  enable  the  computer  to  control 
the  motors,  an  interface  was  needed  between  the  indexers  and  the  MINC. 
In  order  to  understand  the  workings  of  the  interface,  one  must  know 
something  about  the  requirements  of  the  PIM's  and  the  MINC. 
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The  PIM's  require  inputs  to  determine  5  decimal  digits  and  a  sign 
(for  the  number  of  steps  and  the  direction  of  motion),  and  an  input  to 
start  the  indexed  move.  The  PIM's  provide  a  "done"  signal  on  another 
line.  Each  digit  is  represented  in  binary  coded  decimal  by  the  states  of 
4  leads  (8-4-2-1).  Rather  than  use  4  leads  for  each  digit,  all  the  dig¬ 
its  use  the  same  BCD  lines  (the  sign  uses  line  8)  and  6  other  control 
(COM)  lines  determine  the  place  of  the  digit  (1's,  10's,  etc.  and  the 
sign).  All  lines  are  active  low  except  the  sign  (SN)  line  which  is 
active  high.  A  PIM  operates  in  the  following  fashions  The  user  signals 
the  PIM  by  pulling  the  INDEX  line  low.  This  causes  the  SN  and  DONE  lines 
to  go  high  (see  timing  diagram  in  Fig.  47).  When  the  INDEX  line  is 
allowed  to  go  high  again,  the  PIM  samples  the  8  line  for  the  sign.  A 
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short  time  later,  the  SN  line  returns  to  its  tow  state.  The  Pin  then 
sequentially  strobes  each  of  the  5  other  COM  lines  and  samples  the  BCD 
lines  on  the  rising  edge  of  the  con  signal  to  determine  each  digit. 
After  the  last  digit  is  acquired,  the  Pin  starts  the  motor  and  causes  it 
to  turn  the  required  number  of  steps.  Uhen  the  move  is  completed  the 
DONE  line  is  pulled  low  again. 

The  niNC  Digital  Out  unit  allows  the  user  to  output  16  bits  under 
computer  control.  In  addition,  two  lines  are  available  which  signal  uhen 
the  16  bits  are  available  on  the  lines,  and  a  reply  line  is  available  to 
signal  the  computer.  The  design  problem  facing  me  was  to  use  these  lines 
to  send  the  number  of  steps  to  the  Pin's  and  initiate  a  move. 

It  was  clear  that  somehow  the  BCD  values  for  each  digit  would  have  to 
be  made  available  to  the  Pin's  on  demand.  This  could  be  done  either  by 
having  the  niNC  respond  to  each  of  the  COn  signals  with  the  correct  val¬ 
ues  or  by  having  these  values  stored  in  an  auxiliary  memory  addressed  by 
the  Con  lines.  Because  the  Pin's  strobed  the  con  lines  at  a  fixed  rate 
with  no  dependence  on  any  sort  of  reply  line,  it  was  better  to  use  an 
auxiliary  random  access  memory  (RAM  to  minimize  the  possibility  of  tim¬ 
ing  problems.  The  niNC  could  then  write  data  into  the  RAn  at  any  rate, 
and  the  Pin's  could  read  the  data  at  their  own  rate. 

A  block  diagram  and  schematic  of  the  final  interface  circuit  are 
shown  in  Pigs.  48  and  49.  The  niNC  writes  the  4  bits  of  data  for  each  of 
the  6  digits  into  RAn  locations  addressed  by  4  of  the  lines.  The  niNC 
then  signals  the  appropriate  Pin  using  one  of  two  INDEX  lines.  That  Pin 
strobes  its  COn  lines,  which  feed  an  encoder  which  produces  the  address 
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corresponding  to  the  given  COM  line.  Part  of  the  address  places  the  RAM 
in  read  mode,  and  the  4  bits  addressed  appear  on  the  RAN's  output  lines, 
which  are  connected  to  the  BCD  lines  of  both  Pin's.  The  DONE  signal 
produced  by  the  PIN's  goes  back  to  the  NINC.  telling  the  program  the 
move  has  been  completed.  In  the  actual  circuit  many  of  the  lines  had  to 
be  buffered  and  the  duration  of  the  "data  ready"  strobe  produced  by  the 
niNC  lengthened. 
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Appendix  C 


OATA  ACQUISITION  PROGRAMS 


What  follows  are  the  source  listings  for  ths  BASIC  programs  written  to 
oslibrsts  ths  hot  film  probe  end  acquire  velocity  profiles. 

20  REfl*  Program  to  calibrate  hot  films  under  partial  computer  control.  User 
30  REM*  must  move  probe  and  monitor  speed;  computer  samples  anemometer* 

40  REM*  requests  speed,  and  stores  calibration  date.  Computer  e*ope  sampling 
50  REM*  if  anemometer  output  falls  belou  a  specified  percentage  of  the  mean. 
60  REM*  If  user  inputs  a  speed  of  aero*  that  calibration  point  is  discarded. 
70  REM*  Flow  speed  is  the  nominal  velocity  in  the  channel  during  calibration. 
SO  REM*  Results  are  graphed;  olesr  soreen  by  typing  any  character. 

Mas  P  —  ««“  Hmmm  “MmaM  **  mmmmM»m«m«MmmmmMMMi«i«mmmMmM«mmm»mM»mm»i«M».*mgnMM  SMMiiMBMMMMMIS 
H ^ ^WWWWWW WWWWW WWWWWWW WWWWWWW WWWWWWWWWWWWW WWW WWW WWWWWWWWWWWWWW WWW WWWWIrWWWW 

100  Din  V(20).U(20) 

no  oin  xciooo) 

120  PRINT  "OUTPUT  FILE"; 

130  LINPUT  FS 

140  OPEN  FS  FOR  OUTPUT  AS  FILE  1 
ISO  PRINT  "FLOM  SPEE0"; 

160  INPUT  U0 

170  PRINT  "SAMPLING  FREQUENCY  AND  DROPOUT  LIMIT"; 

ISO  INPUT  F.L 
190  K=0  \  S1=0 

200  PRINT  "BEGIN  BY  TYPING  RETURN"; 

210  LINPUT  S0.AS 
220  SET.GAIN(.I.S) 

230  A1N('C0NT1NU0US'.X().F.1/F.S) 

240  WAIT_F0R_DATA(X() »  I ) 

250  FOR  J=I  TO  I-1+F/2 
260  SisSHXU) 

270  K*K*1 
2S0  M=S1/K 

290  IF  X( J)/M<L  SO  TO  330 
300  NEXT  J 

310  PRINT  USING  "••.••".M 
320  GO  TO  240 
330  PRINT  "SPEED"; 

340  TERMINATE!' IMMEDIATE' «XO) 

350  INPUT  S  \  IF  S*0  00  TO  190 
360  S*S+U0 

370  PRINT  II. USING  "•».»», II. It", S.M 
300  PRINT  "ANOTHER  SAMPLE"; 

390  LINPUT  AS 


400  IF  SEC*(A»,  1,  60  TO  190 

410  PRINT  61. USING  "M.  »»,  M.  M",  0, 0 
420  CLOSE  «1 

430  OPEN  FS  FOR  INPUT  AS  FILE  2 

440  FOR  1=1  TO  20 

4S0  INPUT  «2,U(I),V(I> 

460  IF  U(I)=0  60  TO  480 
470  NEXT  I 
480  N=I- 1 

490  GRAPH(~-HLINES',’#N,U<t).V<n} 

500  LINPUT  60. A6 
510  0 I SPLAY_CLEAR 
520  END 
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20  REM  *  Program  to  take  mean  and  turbulenoe  Intensity  profiles*  to  be  used 
30  REM  •  with  the  DEC  MINC-11  with  PREAMP*  A/0*  DIGITAL  OUT.  and  CLOCK  LAB 
40  REM  *  MODULES,  along  with  the  speolally  designed  stepper  motor  oontrol 
50  REM  *  and  Interface  units.  U  refers  to  velocities,  V  to  voltages. 

70  DIM  L(2) >T(  1 ) ,X( 1000) 

60  DIM  Y C 1  > > U ( 50 > , V ( 50 > 

90  REM  initialize  axes  indices*  step-distanoe  conversion  factors*  and 
100  REM  backlash  constant. 

110  READ  Z0,Y0,C1(1)*C1(2)*B1 

120  DATA  1,2,80,-4000,. 1 

t30  LS»"««.«»,  «•«.«•*  ••.«#,  I.IM" 

140  REM  get  lists  of  coordinates  for  traverse 
150  PRINT  "Y  AND  Z  COORDINATE  FILES'*) 

160  INPUT  AS, BS 

170  OPEN  AS  FOR  INPUT  AS  FILE  1 

180  OPEN  BS  FOR  INPUT  AS  FILE  2 

190  REM  get  name  of  output  file  for  data 

200  PRINT  "OUTPUT  FILE"; 

210  LINPUT  AS 

220  OPEN  AS  FOR  OUTPUT  AS  FILE  4 
230  REM  Initialize  probe  looation 
240  PRINT  "MOVE  PROBE") 

250  LINPUT  AS 

260  IF  SEGSCAS, 1* 1)<)"Y"  60  TO  350 
270  PRINT  "Y.Z  COORDINATES") 

280  INPUT  Y.Z 

290  C=Y0  \  P=Y  \  OOSUB  3000 
300  C=Z0  \  P=Z  \  GOSUB  3000 
310  PRINT  "SET  Y  ORIGIN") 

320  LINPUT  AS 

330  IF  SE6S(AS,1,1)="Y"  THEN  U2)*0 

340  GO  TO  240 

350  PRINT  "TAKE  PROFILE"; 

360  LINPUT  AS 

370  IF  SE6S ( AS ,  1 , 1 )  ="Y"  CO  TO  1000 
380  REM  terminetion  sequence 
390  PRINT  "HALT"; 

400  LINPUT  AS 

410  IF  SEGSCAS,  1 , 1  )0"Y"  60  TO  240 

420  REM  data  file  terminated  with  special  sequence 

430  PRINT  14, USING  LS, -1,0, 0,0 

440  CLOSE  14 

450  REM  reposition  probe  at  knoun  looation 
460  CsY0  \  P*1 2  \  60SUB  3000 
470  C=Z0  \  P=0  \  80SUB  3000 
480  STOP 

1000  REM  set-up  calibration 
1010  OOSUB  6000 

1020  REM  initialize  graphlos  to  display  data 
1030  HIND0H("EXACT", . . 12, 16) 

1040  GRAPH("-HLIN£S", 1>X(Q)*Y(0)) 

1050  F»6 
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1060 
1070 
1080 
1090 
1100 
1110 
1120 
1130 
1140 
1150 
1160 
1170 
1180 
1190 
1 200 
1210 
1220 
1230 
1240 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1 510 
3000 
3010 
3020 
3030 
3040 
3050 
3060 
3070 


PRINT  "SAMPLING  FREQUENCY  IS  6  HZ,  WHAT  IS  DURATION*! 
INPUT  T 
PRINT  ** 

PRINT  "  Y  Z  U  <U'*U'>* 

PRINT  * - * 

REM  get  coordinates  to  Move  probe  to 
INPUT  82, Z 

REM  coordinate  list  terminator  Is  -99 
IF  Z=-99  00  TO  390 
INPUT  #1 , Y 
IF  Y>-99  80  TO  1200 

REM  coordinates  Itersted  on  Y  for  eaoh  Z 
RESTORE  »1 
GO  TO  1120 

REM  get  freestream  condition 
GOSUB  5000 

C=Z0  \  P=Z  \  GOSUB  3000 
C=Y0  \  P=Y  \  GOSUB  3000 
REM  sample  a/d 
S1=0  \  S2=0 
SET_GAIN(, 1,8) 

AIN ("CONTINUOUS*, X( ) »2*F , 1/F, 8) 

FOR  1=1  TO  T 
NAIT_FOR_DATA(XO,K1> 

REM  process  data 
FOR  J=K1  TO  K1+F-1 
X=X(J)*« 

IF  X>=V(0)  GO  TO  1360 
PRINT  *0UT  OF  RANGE* 

GO  TO  1400 
FOR  K=t  TO  N 
IF  X<=V(K)  GO  TO  1390 
HEXT  K 

U=U(K-1)+(U(K)-U(K-1)>*(X-V(K-1>)/(Y(K>-Y(K-1>) 

S1=S1+U 


S2=S2+U*U 
NEXT  J 
NEXT  I 

TERMINATE (*IMME0I ATE*) 

REM  calculate  mean  and  Mean  square,  type,  and  plot 
M=S1/(F*T) 

U2=S2/(F*T)-M*M 
PRINT  USING  LG, Y,Z,M,U2 
PRINT  »4, USING  LG,Y.Z,M,U2 
POINT ( , Y,M) 

GO  TO  1150 

REn***»»**« 

REM*  Subroutine  to  Move  channel  C  (*ChannelG)t«P(=Plaoe). 


REM  figure  out  displacement  from  current  location  then  update  position 
0=P-L(C) 


L(C)*P 
S=0*C1(C) 
GOSUB  4000 
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30 >0  IF  D>=0  THEN  RETURN 

3090  REN  correct  for  backlash  If  displacement  is  nagativa 
3100  S=-B1*C1 (C) 

3110  60SU0  4000 
3120  S=-S 
3130  COSUB  4000 
3140  RETURN 

4000  REM********************************************************************** 
4010  REM*  Subroutine  to  index  channel  C(  S  steps. 

4020  ren********************************************************************** 
4030  REM  sat  up  addresses  for  step  nusbers 
4040  FOR  1=0  TO  7  \  Ad >  =  10*d+4)  \  NEXT  I 
4050  0(51=0 

4060  REn  decompose  number  of  steps  into  separate  digits  and  sign 
4070  S1=ABS(S) 

4080  IF  SIOS  THEN  0(5)=9 
4090  FOR  1=4  TO  0  STEP  -1 
4100  D(I)=JNTCS1/10AI) 

4110  S1=S1-DC I )*( 10A1 ) 

4120  NEXT  I 

4130  REN  add  on  address  bits  and  convert  to  binary  for  output 
4140  FOR  1=0  TO  5  \  Ad)=A<I)+Dd)  \  NEXT  1 
4150  A(6)=A(6)*C 
4160  A(7)=A(6) 

4170  FOR  1=0  TO  7  \  NAKE_BCO(A(1 1 .  Ad  1 )  \  NEXT  I 
4180  OOUTCf  AO ,  8, . 0 1 , ) 

4190  REN  send  out  dummy  argument  to  start  indexing  and  wait  for  controller's 
4200  REN  00NE  signal 
4210  D0UT(,T(),2.0) 

4220  RETURN 

5000  REn******************************** ************************************** 
5010  REN*  Subroutine  to  get  freestream  eorreotion  factor. 

5020  REN****************** **************************************************** 

5030  REN  move  probe  to  freestream  conditions 

5040  C=Y0  \  P=12  \  60SUB  3000 

5050  C=Z0  \  P=2  \  60SUB  3000 

5060  V1=0 

5070  REN  get  mean  bridge  output 
5080  SET_GAIN(,1,8) 

5090  AIN(,X().20»F,1/F.8) 

5100  FOR  1=0  TO  20*F-1 
5110  V1=»1+Xd) 

5120  NEXT  I 
5130  V1=V1/(20*F) 

5140  q=vo/vi 

5150  PRINT  USINB  »V1M.*l«'r,V1 
5160  RETURN 

5170  REN********************************************************************** 
5180  REN*  Subroutine  to  read  velocity  off  calibration  curve  using  X. 

5 1 90  REN********************************************************************** 

5200  IF  X>V(N)  80  TO  5270 

5210  IF  X<V(0)  60  TO  5270 

5220  FOR  K=1  TO  N 

5230  IF  X<=V(K)  00  TO  5250 


5240  NEXT  K 

5250  U=U<K-1)«-<U(K)-U<K-1))*(X-V<K-1))/{V<lO-V(IC-1)) 

5260  RETURN 

5270  PRINT  "OUT  OF  RANGE" 

5200  RETURN 

6000  REM************** ******* ************************************************* 
6010  REM*  Subroutine  to  sot-up  calibration  ourv*. 

6020  REM********************************************************************** 
6030  PRINT  "CALIBRATION  FILE"; 

6040  LINPUT  AS 

6050  OPEN  AS  FOR  INPUT  AS  FILE  3 
6060  FOR  1=0  TO  50 
6070  INPUT  #3,U(I),V(I) 

6000  IF  UCI )-0  GO  TO  6100 
6090  NEXT  I 
6100  H=I- 1 

6110  PRINT  "VELOCITY  AT  REFERENCE  POINT"; 

6120  INPUT  UO 

6130  FOR  1=1  TO  50 

6140  IF  U0<=U(I)  GO  TO  6160 

6150  NEXT  I 

6160  V0=V(I-1)+(U0-U(I-1))*(vm-V(I-1))/(UCI)-U(I-l)) 

6170  RETURN 
9000  ENO 
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Appendix  0 

HOT  FILM  NORMALIZATION  PROCEDURE 

Nomenclature 

d  sensor  diameter 

e  bridge  voltage 

ea  bridge  voltage  at  measurement  station 

*r 

en  normalized  bridge  voltage*  ea- - 

*r' 

er  bridge  voltage  at  reference  point 

H  rate  of  heat  transfer  from  sensor  to  fluid 

kaff  an  effective  thermal  conductivity  of  the  fluid*  including  the 
effects  of  dirt  on  the  sensor 

jt  sensor  length 

Nu  Nusselt  number  (defined  below) 

Re  Reynolds  number  (defined  below) 

Ra  sensor  resistance  at  ambient  temperature 

R«  sensor  cold  resistance 

Rop  sensor  operating  resistance 

u  fluid  velocity 

ua  fluid  velocity  at  measurement  station 

ur  fluid  velocity  at  reference  point 

«  temperature  coefficient  of  resistivity  of  sensor 
r  ratio  of  heat  loss  due  to  oonvection  to  total  heat  loss 
v  kinematic  viscosity  of  tha  fluid 
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9a  ambient  temperature  of  fluid 

9C  ambient  temperature  during  cold  resistance  measurement 
9op  operating  temperature  of  sensor 
O'  quantity  at  some  later  time 

Hot  film  anemometers  have  been  widely  used  as  an  inexpensive  means  for 
acquiring  velocity  measurements  in  fluid  flows.  Because  of  their  dura¬ 
bility.  these  probes  are  often  used  in  water,  where  the  relatively  low 
frequency  response  of  the  large-area  sensor  is  not  a  problem.  A  common 
problem  one  encounters  when  using  a  hot  film  probe  is  the  drift  of  the 
probe's  calibration  over  time,  due  to  the  probe's  acute  sensitivity  to 
dirt  on  its  surface  and  because  of  changes  in  the  fluid  temperature.  In 
the  past,  experimenters  have  sought  to  overcome  these  difficulties  by 
carefully  filtering  their  water  and  by  regulating  the  water  temperature 
or  employing  a  temperature  compensation  probe.  Uith  the  wide  availabil¬ 
ity  of  laboratory  microcomputers,  however,  another  approach  is  possible. 
Here,  an  alternate  method  is  described  which  enabled  us  to  obtain 
repeatable  measurements  in  a  low-speed  water  flow  (u<15  cm/s). 

Preliminary  measurements  showed  that  no  reliable  data  could  be 
obtaired  from  the  hot  film  anemometer  we  were  using  unless  something  was 
done  to  account  for  changes  in  the  probe  environment  due  to  dirt  and 
temperature  fluctuations.  Fig.  50  shows  the  anemometer  bridge  output 
over  an  hour's  time  uith  the  probe  operating  at  an  overher'i  ratio  of 
1.08.  The  change  in  water  temperature  over  this  period  was  less  than 
0. 1°C»  so  the  drift  of  the  probe  output  is  attributed  to  the  accumula¬ 
tion  of  dirt.  The  probe  was  held  fixed  in  a  constant  velocity,  10  caw's, 
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Figure  51s  Variation  of  calibration  with  temperature  (simulated 
by  change  in  overheat  ratio). 


For  the  long  running  times  needed  to  acquire  stable,  mean  velocity 
profiles  in  our  lou  speed  flou,  changes  in  the  uater  temperature  might 
also  have  a  significant  effect  on  the  calibration,  especially  uith  the 
lou  overheat  ratios  that  had  to  be  used  to  prevent  bubble  formation  on 
the  sensor.  Fig.  51  shous  calibration  curves  taken  on  different  days. 
The  cold  resistance  of  the  probe  uas  about  0.05  ft  higher  on  the  8th  than 
on  the  13th,  but  the  operating  resistance  of  the  probe  uas  set  at  5.27  ft 
on  the  8th  and  5.25  ft  on  the  13th.  Consequently,  the  probe  uas  operated 
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at  an  ovarheat  ratio  of  1.07  on  th*  0th  and  1.00  on  the  13th.  The 


results  therefore  simulate  uhat  would  happen  to  the  calibration  curve  of 
a  probe  if  the  fluid  temperature  were  to  decrease  and  emphasize  the  sen¬ 
sitivity  of  the  probe  to  changes  in  the  overheat  ratio. 

The  similarity  of  the  curves  in  Fig.  51,  and  the  idea  that  small 
changes  in  the  probe/s  response  due  to  other  factors,  such  as  dirt  build 
up,  might  also  cause  the  calibration  curve  to  shift  slightly,  suggested 
that  a  normalization  procedure  might  be  able  to  compensate  for  changes 
in  the  probe's  response  and  thereby  improve  its  reliability  and  utility. 

The  normalization  procedure  adopted  corrects  the  anemometer's  output 
using  information  about  how  the  probe's  response  has  changed  at  a  loca¬ 
tion  where  the  velocity  is  assumed  not  to  change.  Specifically,  the  nor¬ 
malized  voltage,  en>  used  to  look  up  the  flow  speed  from  a  calibration 
table,  is  obtained  as  follows: 

*  r 

e „=e a  .  ( D . 1 ) 
*r/ 

The  justification  for  this  normalization  procedure  follows. 

Hills’*  shows  that  for  most  situations  of  practical  interest,  the 
heat  transfer  from  a  cylindrical  hot  wire  or  hot  film,  expressed  as  a 
dimensionless  Mussel t  number,  is  a  function  of  only  the  Reynold  number 
viz. 

Nu-f (Re) .  (0.2) 
This  was  the  case  in  our  flow,  Nu  is  defined  as 


’•  Hills,  J.  A.  B.  (1900)  "Hot-wire  snd  hot-film  snemometry*,  Nstionsl 
Maritime  Institute,  May  1980. 
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H 

Nu= -  (0.3) 

Rk«f  fA  (9  op-0  ») 

and 

ud 

R#=— .  (0.4) 

v 

The  sensor  is  heated  by  the  current  which  flows  through  it  due  to  a 
voltage,  e.  Part  of  the  heat  is  lost  through  conduction  to  the  sensor 
supports.  For  hot  films,  this  heat  is  probably  not  negligeable.  The  heat 
lost  to  the  fluid  by  convection  is  therefore  some  fraction  of  the  total 
heat  generated  in  the  sensor » 
e* 

H=y - .  (D.  5) 

Rop 

Since  the  heat  lost  through  conduction  ultimately  ends  up  back  in  the 
fluid,  and  the  temperature  gradient  over  which  the  conduction  occurs  is 
the  same  as  for  convection  from  the  sensor,  y  is  not  expected  to  vary 
significantly.  Substituting  Eqs.  0.3-0. 5  into  Eq.  D.2  gives 
ye1  rudi 

- =f  —  .  (0.6) 

R  opflk  •  f  (8  n  )  U  . 

Note  that  everything  on  the  left  hand  side  of  Eq.  0.6  except  for  e  does 
not  depend  on  u.  We  can  then  write 

where  f  depends  on  the  operating  conditions  of  the  probe  but  not  the 
fluid  velocity.  In  the  discussion  which  follows,  the  variation  of  v  with 
temperature  (about  IK  per  °F)  will  be  assumed  negligeablei  this  was  the 
ease  for  our  flow.  Oirt  accumulating  on  the  sensor  affects  the  effective 
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*Sr-**'  j»****tf,  v,1  »'f' ; 


*/ 


diameter  of  the  sensor,  but  if  we  assume  most  of  the  larger  dirt  parti 
eles  can  be  shaken  off,  then  d  should  remain  fairly  constant. 


Ue  can  now  demonstrate  how  the  normalization  procedure  works.  Sup¬ 
pose  ue  calibrate  the  probe  under  a  certain  set  of  conditions.  The  var¬ 
iation  of  bridge  voltage  uith  fluid  velocity  can  be  represented  by  Eq. 
0.7.  In  particular,  the  voltage  at  a  reference  point  can  be  expressed  as 


— Frl- 


(0.8) 


"Frl 


(0.9) 


Nou  suppose  at  some  later  time,  ue  measure  a  voltage  under  different 
operating  conditions  at  some  point  in  the  flow.  The  voltage  and  velocity 
there  are  related  by 
fu.dl 
v 

where  ♦'  is  unknown  as  is  ua.  t'  can  be  determined  by  returning  to  the 
reference  point  and  measuring  the  voltage  corresponding  to  ur  under  the 
new  operating  conditions: 
furdi 
v 

Combining  Eqs.  0.8  and  0.10  ue  obtain 
.2 


♦'eP'* 


"It) 


(0.10) 


♦'=*- 


*r 


(0.11) 


•  t 


•  r 


Substituting  Eq.  0.11  into  Eq.  0.9  yields 


te 


*r*  fuadT 

'* — =f—l 
er'*  l  v  J 


(0.12) 


Comparing  Eq.  0.12  with  Eq.  0.7,  we  see  that  everything  in  Eq.  0.12  is 
known  except  for  uM,  which  can  then  be  found  by  Inverting  the  funotion 
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f.  In  practice  what  ua  do  is  look  up  the  velocity  corresponding  to  the 
normalized  voltage  using  the  old  calibration  curve.  For  example,  by 
applying  this  procedure  to  one  of  the  curves  in  Fig.  SI  the  third  curve 
results.  Consequently,  the  2  cm/s  error,  which  would  have  occurred  had 
the  unnormalized  voltages  been  used  with  the  older  calibration  curve,  is 
substantially  reduced. 

It  is  instructive  to  examine  the  sensitivity  of  •  to  various  parame¬ 
ters.  Recall  that 
V 

*= - .  (D. 13) 

Ropfk •f  f X ( 0 0p-6 1 ) 

Rop  is  generally  not  varied  by  more  that  IX  in  order  to  achieve  the  same 
overheat  ratio  on  any  given  day.  ksff  can  have  a  large  effect  on  ♦  .  as 
demonstrated  by  Fig.  50  if  dirt  is  allowed  to  accumulate  on  the  sensor. 
In  the  present  tests,  the  constant  movement  of  the  probe  between  the 
measurement  stations  and  the  reference  point  appeared  to  keep  k«ff 
fairly  constant.  However,  a  thin  film  of  dirt  can  accumulate  with  time 
and  slowly  lower  k«ff. 

9 op  depends  on  the  resistance  value  at  which  the  sensor  is  operated. 
The  resistance  of  the  sensor  varies  approximately  linearly  with  tempera¬ 
ture: 

Rop=Kc(Ha(e0p-ae)].  (0.14) 

«  is  about  0.002  per  °C  for  the  probes  used  in  this  study.  Solving  for 
9«p  and  substituting  into  Cq.  0.13, 

y 

t  =  - . 

rfRop  -p  i  (0.15) 

*op»kef f* JJ -  “  lj“  +  0e”0»j 
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In  our  study  0e-6«  was  generally  less  than  t°C  if  the  cold  resistance  of 
the  probe  was  measured  at  the  ambient  temperature  at  the  beginning  of  a 
run.  The  other  term  in  the  brackets  represents  the  difference  between 
the  temperature  of  the  sensor  and  the  temperature  at  which  the  cold 
resistance  of  the  sensor  is  measured  and  is  on  the  order  of  50°C.  This 
term  turns  out  to  be  quite  sensitive  to  the  overheat  ratio.  Rop/Re.  A 
drift  or  other  small  change  in  R0|>  of  0.01  ft  1 0 . 250  results  in  a  1% 
change  in  this  temperature  difference  term.  Fortunately,  as  long  as  If 
is  independent  of  the  fluid  velocity,  its  variations  aree  compensated 
for  by  the  normalization  procedure  as  shown  above. 

It  should  be  noted  that  this  normalization  procedure  also  has  some 
nonlinear  effect  on  the  measured  velocities  if  the  actual  reference 
velocity  in  the  freestream  changes.  Because  of  the  nonlinearity  of  the 
basic  correlation,  Eq.  0.2,  velocity  ratios  are  not  scaled  linearly  with 
voltage  ratios.  For  example,  an  increase  in  the  freestream  velocity 
over  the  nominal  reference  velocity  will  cause  the  anemometer  to 
increase  the  voltage  applied  to  the  probe.  The  normalization  procedure 
scales  all  of  the  measured  voltages  so  that  the  voltage  at  the  reference 
point  after  normalization  remains  constant.  The  measured  velocity  at  the 
reference  point  therefore  remains  constant.  However,  it's  not  as  if  all 
of  the  other  actual  velocities  were  scaled  on  the  freestream  velocity 
too.  The  normalization  procedure  multiplies  the  measured  voltages  by  a 
fraction,  er/er'>  Due  to  the  nonlinearity  of  the  calibration  curve,  this 
is  not  equivalent  to  multiplying  all  the  velocities  by  Up/Ur'.  It  can  be 
shown  that,  depending  on  the  shape  of  the  calibration  curve,  velocities 
other  than  the  reference  velocity  are  either  overcorrected  or  undercor- 


rected.  For  example*  if  the  actual  freestream  velocity  increasea  by  10X* 
the  normalization  reduces  the  measured  voltages  by  an  amount*  say  3X» 
corresponding  to  a  9. IX  reduction  in  the  measured  velocity  at  the  refer¬ 
ence  point.  However*  velocities  at  other  locations  will  be  reduced  by 
more  or  less  than  9. IX  since  a  3X  change  in  voltage  corresponds  to  dif¬ 
ferent  changes  in  velocity  at  different  locations  on  the  calibration 
curve.  Since  the  freestream  velocity  in  the  present  study  was  carefully 
regulated*  this  side  effect  of  the  normalization  procedure  was  probably 
negl igeable. 
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Appendix  E 


SPECIALIZED  TOOLS 

Several  tools  uere  designed  and  fabricated  to  perform  special  tasks  in 
the  assembly  and  operation  of  the  water  channel.  Shown  in  Fig.  52  are 
(a)  a  flexible  tip*  ball  driver  for  placing  and  driving  the  socket  head 
screws  used  to  hold  the  inner  walls  in  place.  (b>  a  remote  manipulator 
for  retrieving  objects  dropped  into  the  channel*  (c)  a  hydrogen  bubble 
wire  holder,  and  (d)  a  modified  syringe. 

The  ball  driver  was  made  by  attaching  a  conventional  ball  driver  tip* 
taken  from  a  screwdriver  handle*  to  a  miniature  universal  joint  and  an 
extra  long  handle.  A  piece  of  vinyl  tubing  slipped  over  the  joint  gave 
the  joint  some  rigidity  so  that  the  tip  could  be  positioned  as  desired, 
besides  being  easier  to  use  than  slot  or  Phillips  head  screwdrivers* 
ball  drivers  also  allow  you  to  temporarily  hold  on  to  a  screw  because  of 
the  way  in  which  the  ball  fits  into  the  socket  head  of  the  screw. 

The  remote  manipulator  was  fashioned  from  a  bicycle  brake  outfitted 
with  modified  caliper  arms.  The  manipulator  possessed  enough  dexterity 
to  pick  up  thin  plastic  rulers  dropped  on  the  channel  floor. 

The  hydrogen  bubble  wire  holder  was  designed  to  allow  easy  replace¬ 
ment  of  broken  wires  while  also  affording  the  safety  needed.  The  holder 
was  fabricated  from  nonconduotive  plexiglas.  The  spacing  between  the  two 


Special ized 


prongs  was  adjustable.  The  electrical  leads  running  doun  the  center  of 

the  stem  terminated  at  small  screws  sticking  out  of  the  base  of  each 

prong.  A  pair  of  nuts  on  each  screw  were  used  to  clamp  the  wires  in 

* 

place.  Parafin  wax  applied  with  a  brush  insulated  terminals  from  the 
water.  The  vinyl  tube  wound  around  the  stem*  corkscrew  fashion*  kept  the 
stem  from  oscillating  in  the  flow  by  disrupting  the  vortices  being  shed. 

Finally,  the  modified  syringe,  used  for  injecting  dyes,  was  fabri¬ 
cated  from  two  syringe  bodies  joined  with  a  plexiglas  adapter. 


Appendix  F 


FILM  NARRATIVE 

"Turbulent  Boundary  Layers  on  Concave  Walls"  is  a  16  mm,  color*  silent 
film*  430'  long*  which  runs  for  12  minutes  if  projected  at  24  frames  per 
second  (sound  speed).  The  film  was  produced  by  Albert  H.  Jeans  and  James 
P.  Johnston  of  the  Thermoseiences  Division  of  the  Dept,  of  Mechanical 
Engineering,  Stanford  Univ.  in  1981.  The  film  gives  some  background 
information  and  then  depicts  systematic  visual  surveys  of  the  floufield 
using  dyes  and  hydrogen  bubbles.  A  scene-by-scene  narrative  follows. 

1.  Title  and  producers. 

2.  Schematic  plan  vieu  of  the  concave  flow  facility.  Uater  enters 
from  an  underground  sump  through  a  perforated  pipe.  The  flou  pas¬ 
ses  through  a  perforated  plate*  honeycomb*  and  screens  to  achieve 
uniformity  and  lou  turbulence  levels.  A  nozzle  with  a  contraction 
ratio  of  4s 1  accelerates  the  flou  up  to  the  nominal  running  speed 
of  15  cm/s.  The  laminar  boundary  layers  are  tripped  about  1  m 
downstream  of  the  nozzle  exit.  The  flow  passes  through  a  90°  bend 
and  a  2.5  m-long  recovery  region  before  returning  to  the  sump. 

3.  Color  diagram  of  the  test  section.  The  boundary  layers  are  about 
7.5  cm  thick  at  the  start  of  curvature  with  a  momentum  thickness 
Reynolds  number  of  about  1300.  The  channel  depth  is  about  1.1  m 


and  the  radius  of  curvature  1.36  m.  The  opposing,  convex  uall  was 
designed  to  minimize  the  streamuise  pressure  gradient  at  the  test 
ual  1 . 

4.  Diagram  of  dye  injection  technique.  To  visualize  the  flou  near 
the  uall.  dye  is  injected  through  narrou  slots  into  the  viscous 
sublayer  and  lit  from  behind. 

5.  Diagram  of  hydrogen  bubble  technique.  To  visualize  the  flou  auay 
from  the  uall  hydrogen  bubbles  are  used  to  mark  the  flou.  The 
diagram  ahous  the  effective  vieuing  angle  uhich  uill  be  seen  in 
the  movies.  In  reality,  this  vieuing  angle  is  accomplished  using 
a  mirror  placed  in  the  flou,  dounstream  of  the  bubble  uires. 

6.  The  dye  movies  uere  shot  at  24  frames  per  second,  the  hydrogen 
bubble  movies  at  45  frames  per  second. 

7.  A  scale  ruled  in  centimeters  appears  on  the  uall  in  many  of  the 
scenes.  Diagram  shouing  the  location  of  the  first  region  to  be 
visualized:  the  flat  plate  region  upstream  of  the  curve. 

6.  This  region  is  characterized  by  uniformly  distributed  streaks  and 
bursts  near  the  uall. 

9.  Dye  is  being  injected  through  tuo  slots  about  20  cm  apart.  He  see 
the  usual  streaky  structure  characteristic  jf  turbulent  boundary 
layers.  The  nondimensional  streak  spacing  X*  agrees  uith  the  gen¬ 


erally  quoted  value  of  100 .  Velooity  measurements  made  in  this 
region  also  indioate  a  normal  turbulent  flatplate  boundary  layer. 


10.  Single  bubble  Mire  et  y=0.6  cm  (y*=44),  pulsed  continuously.  Typ¬ 
ical  eddies  are  1  or  2  cm  in  size  although  larger  ones  occur 
occasionally.  Streaks  appear  as  lou-speed  regions  uhich  oscillate 
back  and  forth.  Bursting  disrupts  the  bubble  lines  soon  after 
they  are  generated. 

11.  Single  Hire  at  y=1.0  cm  (y*=73)»  pulsed  intermittently.  This 
pulsing  pattern  does  a  better  job  of  depicting  the  eddies  as  they 
travel  downstream.  Note  the  appearance  of  streamuise  vorticity. 

12.  Dual  bubble  wires  are  always  separated  by  4  cm  and  are  pulsed 
intermittently. 

13.  Dual  Hires,  louer  wire  at  y=1.0  cm.  Note  the  lack  of  correlation 
between  the  two  bubble  lines,  another  indication  that  the  eddies 
in  this  region  are  not  very  large.  Turbulent  activity  seems  to  be 
diminishing  farther  from  the  wall. 

14.  Dual  wires,  louer  Hire  at  y=5. 0  cm.  The  upper  wire  is  now  outside 
of  the  boundary  layer.  The  bubble  lines  from  the  upper  wire 
remain  Intact  aside  from  occasional  eddies  from  below. 

15.  Diagram  showing  the  region  in  the  curve  uhich  was  observed,  a 
region  starting  at  about  30°. 

16.  The  flow  in  this  region  appears  to  be  dominated  by  randomly  dis¬ 
tributed  large-scale  sweeps  (inflows)  and  ejections  (outflows). 

17.  Dye  visualization.  Note  hew  the  dye  eotleets  into  large  stream¬ 
ers.  whioh  are  spaced  about  10  cm  apart  on  the  average.  The 
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boundary  layer  thickness  is  also  about  10  cm  here.  The  streamers 


are  actually  formed  by  sweeps  from  the  outer  flow  which  push 
aside  the  fluid  near  the  wall.  The  sweeps  appear  randomly  in  pos¬ 
ition  and  time  and  seem  to  greatly  reduce  the  bursting  activity 
seen  in  the  flatplate  boundary  layer. 

18.  Bubble  wire  at  y=0.6  cm.  pulsed  continuous) y .  Here  again  we  see 
the  wall  fluid  being  pushed  into  piles  by  sweeps  from  the  outer 
flow.  Bursting  is  greatly  diminished  where  the  sweeps  impinge  on 
the  wal 1 . 

19.  Bubble  wire  at  y=1.0  cm,  pulsed  intermittently.  The  vertical 
movement  of  the  fluid  is  much  greater  here  than  in  the  flat  plate 
region.  Typical  eddies  are  about  4  or  5  cm  in  size. 

20.  Dual  wires,  lower  wire  at  y=1.0  cm.  (lotions  depicted  by  the  bub¬ 
ble  lines  appear  highly  correlated,  indicating  that  the  large- 
scale  eddies  extend  a  long  distance  through  the  boundary  layer. 
The  large-scale  motions  also  appear  to  be  promoting  a  lot  of  mix¬ 
ing  and  turbulence  production  far  from  the  wall. 

21.  Oual  wires,  lower  wire  at  y=5.0  cm.  The  upper  wire  is  near  the 
edge  of  the  boundary  layer.  Many  of  the  ejections  appear  quite 
strong  at  the  upper  wire  position. 

22.  Mire  at  9.0  cm,  pulsed  continuously.  This  visualization  gives  you 
a  better  Idea  of  the  planar  distribution  of  the  ejections.  The 
ejections  are  irregularly  shaped  and  randomly  distributed.  Note 
the  Intense  longitudinal  vortioity  occasionally  generated  at  the 
interfaces  between  the  ejections  and  the  outer  flow. 


146 


23.  The  flow  displays  the  same  qualitative  features  at  other  stations 
around  the  curve,  the  differences  being  that  the  large-scale 
motions  appear  to  intensify  and  grow  somewhat  larger. 

24.  In  the  recovery  region  downstream  of  the  curve,  bursting  is 
greatly  suppressed  over  large  regions. 

25.  Diagram  showing  region  visualized  in  recovery  region. 

26.  Dye  visualization.  Note  large  regions  of  low  turbulent  activity. 
Perhaps,  the  large-scale  motions  in  the  curve  have  removed  so 
much  fluid  from  the  wall  that  a  new  boundary  layer  is  growing 
here. 

27.  Bubble  wire  at  0.6  cm,  pulsed  continuousl y .  The  large  patches  of 
laminar-like  fluid  alternate  with  patches  of  more  conventional 
bursting  activity. 

28.  Dual  wires,  lower  wire  at  y=1.Q  cm.  Remnants  of  the  large-scale 
motions  still  appear  quite  strong.  Chunks  of  laminar-like  fluid 
also  appear.  These  may  be  remnants  of  the  potential  core  or  zones 
where  rel aminar ization  has  occurred;  at  this  time  no  definite 
conclusions  can  be  drawn. 

29.  Dual  wires,  lower  wire  at  y=5.0  cm.  Large  scale  motions  dominate 
the  entire  flowfield.  The  scale  and  inertia  of  these  motions  evi- 
dentatly  makes  them  very  resistant  to  decay. 

30.  Single  wire  at  y=9,0  cm,  pulsed  continuously.  Again,  we  see  large 
patches  of  laminar-tike  fluid.  Velooity  and  turbulence 


Measurements  show  no  sign  of  a  potential  core  by  the  75° 
location.  Visualization  gives  us  a  better  idea  of  uhat  is  happen¬ 
ing. 

31.  This  work  was  funded  by  the  Air  Force  Office  of  Scientific 


Research. 
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DATA  SETS 

The  data  sets  used  in  the  graphs  are  fisted  below  in  two-column  format. 
Except  where  noted*  the  anemometer  was  sampled  at  6  Hz  for  a  period 
T =300  s  to  obtain  the  mean  velocity  and  turbulence  intensity  at  each 
location.  Y0  is  a  shift  applied  to  the  y-eoordinates  in  order  to  par¬ 
tially  compensate  for  positioning  errors.  Y0  was  determined  by  examining 
the  near-wall  data  and  estimating  the  probable  location  of  the  wall  by 
extrapolation.  The  value  of  Y0  given  in  the  data  sets  is  the  amount  that 
was  added  to  the  nominal  y-coordinates  to  produce  the  coordinates  shown. 
THETA*  expressed  in  degrees*  is  the  location  in  the  curve  at  which  some 
of  the  data  sets  were  taken.  M  is  the  channel  width  at  the  streamuise 
location  at  which  the  cross-channel  profiles  were  taken.  TEMP  is  the 
water  temperature*  in  degrees  Fahrenheit*  at  the  time  the  run  was 
started.  The  integer  following  the  title  of  eaeh  data  set  is  the  number 
of  points  in  that  set. 
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S.700E+00 

4. OOOE+OO 

0. OOOE-OI 

5.74QE+00 

0.0Q0E-03 

4.900E-01 

2.000E+01 

9. 270E+00 

4.400E+00 

0. OOOE-OI 

S.7S0E+O0 

S.OOOE-03 

6. 200E-01 

2.QOOE+OI 

9.050E+00 

150 


7. 9001-01 

2.00QE+01 

1.00SE+01 

1 . 200E+00 

O.OOOE-OI 

1 . 107E+01 

1 . 373E+00 

1.010E+00 

2 . OOOE+O 1 

1.064E+01 

1 . 610E+00 

O.OOOE-OI 

1. 144E+01 

9.670E-01 

1 . 270E+00 

2. OOOE+O 1 

1.096E+01 

2.0I0E+00 

O.OOOE-OI 

1 . 191E+Q1 

1 . 002E+00 

1.600E+00 

2. OOOE+O 1 

1 . 126E+01 

2.480E+00 

O.OOOE-OI 

1 . 236E+01 

9. S40E-01 

2.000E+00 

2.000E+01 

1. 173E+01 

3. 170E+00 

0. OOOE-OI 

1.303E+01 

0.620C-II 

2. 470E+00 

2. OOOE+O 1 

1. 190E+01 

3.680E+00 

O.OOOE-OI 

1 . 347E+01 

7.490E-M 

3. teoE+oo 

2. OOOE+O 1 

1.2S0E+01 

4.290E+00 

O.OOOE-OI 

1.377E+01 

6.100E-M 

3.670E+00 

2. OOOE+O 1 

1.273E+0! 

4. 980E+00 

O.OOOE-Ot 

1.4IOE+01 

3. 090E-9I 

4.280E+00 

2.000E+01 

1.333E+01 

S.980E+00 

O.OOOE-OI 

1.461E+01 

1.S60E-0I 

4.970E+00 

2. OOOE+O 1 

1.3S1E+Q1 

7.900E+00 

O.OOOE-OI 

1.S15E+01 

1.000E-02 

5.970E+00 

2.00QE+01 

1.429E+01 

9.980E+00 

e.oooE-ot 

1.482E+01 

7.000E-03 

7. 970E+00 

2. OOOE+O 1 

1.479E+01 

1.198E+01 

O.OOOE-OI 

1.401E+01 

6.000E-03 

HU  All  5/13/81  XM00 

CM 

FUT2  S/11/81  X*400 

CM 

2=10  CM  Y0*.05  CM  TEHP*76 

2*0  YO*-. 

02  TEHI»*76 

22 

20 

Y  CM 

2  CM 

<U>  CM/O 

Y  CM 

2  CM 

<U>  CM/S 

1. SOOE-OI 

1. OOOE+O 1 

S.730E+00 

1.2Q0E-01 

O.OOOE-OI 

S.  150E+00 

3.264E+00 

1.700E-01 

1. OOOE+O 1 

6. 360E+00 

1.S00E-01 

O.OOOE-OI 

5. 880E+00 

3.273E+00 

1.900E-01 

1. OOOE+O 1 

6.630E+00 

1. 800E-01 

O.OOOE-OI 

6.690E+00 

3.433E+00 

2.200E-0I 

t.OOOC+OI 

7. 100E+00 

2.200E-0I 

O.OOOE-OI 

7. S30E+00 

3. 104E+00 

2.500E-01 

1.000E+01 

7 . 590E+00 

2.800E-01 

0. OOOE-OI 

0. 090E+00 

2.747E+00 

2.900E-01 

1. OOOE+O 1 

8.240E+00 

3.800E-01 

O.OOOE-Ot 

9.040E+00 

2.474E+00 

3.S00E-01 

1.000E+01 

8. 680E+00 

5.000E-01 

0. OOOE-OI 

9.6S0E+0O 

2. 122E+00 

4.500E-01 

1.000E+01 

9. 390E+00 

6.300E-01 

0. OOOE-OI 

1.007E+01 

1.609E+00 

S.700E-01 

1.000E+01 

9.700E+00 

8.000E-01 

0. OOOE-OI 

1.0S0E+01 

1 . 561E+00 

7.000E-01 

1.000E+01 

1.0Q7E+01 

1.020E+00 

0. OOOE-OI 

1. OOOE+O 1 

1 . 224E+00 

8.700E-01 

1.000E+01 

1.042E+01 

1 . 280E+00 

0. OOOE-OI 

1. 112E+01 

1 . 281E+00 

1 . 090E+00 

1. OOOE+O 1 

1.066E+01 

1.610E+00 

0. OOOE-OI 

1. 142E+01 

1.200E+00 

1 . 350E+00 

1.000E+01 

1. 104E+01 

2.010E+00 

0. OOOE-OI 

1.177E+01 

1.037E+00 

1.680E+00 

1. OOOE+O 1 

1. 150E+01 

2.480E+00 

0. OOOE-OI 

1.23SE+01 

9. 130E-01 

2.080E+00 

1.000E+01 

1 . 1S1E+01 

3. 170E+00 

0. OOOE-OI 

1.291E+01 

9.260E-01 

2. 550E+00 

1.000E+01 

1.224E+01 

3.080E+00 

0. OOOE-OI 

1.333E+01 

0.790E-01 

3.240E+00 

1. OOOE+O 1 

1. 207E+01 

4. 290E+00 

9.00DE-01 

1.307E+01 

6.230E-01 

3.750E+00 

1.0QQE+01 

1.314E+01 

4.900E+00 

0. OOOE-OI 

1.427E+01 

3. 460E-01 

4.360C+00 

1.000E+01 

1.366E+01 

S.980E+00 

0. OOOE-OI 

1.464E+01 

1.700E-01 

S.050E+00 

1.000E+01 

1.417E+01 

7.9S0E+00 

0. OOOE-OI 

1.489E+01 

1 . 900E-02 

6.DS0E+00 

1.000E+01 

1.440E+01 

6.0S0E+00 

1. OOOE+O 1 

1.484E+01 

FUT5  S/20/01  X*400  CM 

2*0  CM  TEHPM1 

run  s/ii/«i  xmoo 

CM 

26 

2*0  CM  Y0«-.02  CM  TEHT«76 

Y  CM 

2  CM 

<0>  CH/S 

<U'«U'> 

22 

1.000E-81 

0. OOOE-OI 

4.070E+00 

2.S32E+00 

Y  CM 

2  CM 

<U>  CM/9 

1 . 200E-01 

0. OOOE-OI 

4.920E+00 

2. 5S4E+00 

1.200E-01 

0. OOOE-OI 

S.200E+00 

2.935E+00 

1.400E-01 

0. OOOE-OI 

5. 500E+00 

3. 100E+00 

1.500C-01 

0.000E-01 

$. 140E+00 

3.435E+00 

1.700E-01 

0. OOOE-OI 

S.910E+00 

3. 113E+00 

1. OOOE-OI 

O.OOOE-OI 

6. 730E+00 

3.C11E+00 

2. 000E-01 

0. OOOE-OI 

6.600E+00 

3. 47SE+00 

2.200E-01 

0. OOOE-OI 

7.3SOE+00 

3. 379E+00 

2.400E-01 

0. OOOE-OI 

7. 180E+00 

3. 155E+00 

2.000E-01 

O.OOOE-OI 

8.320E+00 

2.I73E+O0 

3.000E-01 

0. OOOE-OI 

0.010E+00 

3.069E+00 

3. OOOE-OI 

O.OOOE-OI 

9.000E+09 

2. 423E+00 

4.000E-01 

o.eeoE-oi 

0.960E+00 

2. 5S1E+00 

5.000E-91 

O.OOOE-OI 

9.670E+00 

2.010E+00 

S.200E-0I 

0. OOOE-OI 

9.490E+00 

1.922E+00 

6. 300E-01 

O.OOOE-OI 

1.017E+01 

1.711E+00 

6. S00E-01 

0. OOOE-OI 

9. 720E+00 

1. 5S2E+00 

0.000E-01 

O.OOOE-OI 

1.039E+01 

1 . 488E+00 

0.200E-01 

0. OOOE-OI 

1.026E+01 

1 . 280E+00 

1.020E+00 

O.OOOE-OI 

1.079E+01 

1 . 332E+00 

1.040E+00 

0. OOOE-OI 

1.0S6E+01 

1 . 120E+00 

-  tsi  - 


1.300E+00 

O.OOOE-OI 

1.090E+01 

1.0S5E+00 

1 . 630E+00 

0.000E-01 

1. 120E+01 

9.070E-01 

2.030E+00 

O.OOOE-OI 

1. 173E+01 

1. OOOE+OO 

2. 500E+00 

O.OOOE-OI 

1.200E+01 

9.390E-01 

3. 190E+00 

O.OOOE-OI 

1.20SE+01 

9.240E-01 

3. 700E+00 

O.OOOE-OI 

1 . 301E+01 

0.430E-Q1 

4.310E+00 

O.OOOE-OI 

1.350E+01 

7.700E-01 

5. OOOE+OO 

O.OOOE-OI 

1.401E+01 

S.920E-01 

0. OOOE+OO 

O.OOOE-OI 

1.45SE+01 

2.7Q0E-01 

0. OOOE+OO 

O.OOOE-OI 

1.493E+01 

3.700E-02 

l.OOOE+OI 

O.OOOE-OI 

1.4SSE+01 

1.000E-02 

1.200E+01 

O.OOOE-OI 

1.473E+01 

1 . 200E-02 

1.400E+01 

O.OOOE-OI 

1.477E+01 

1.000E-02 

l.OOOE+OI 

O.OOOE-OI 

1.471E+01 

2.900E-02 

QUA12  5/13/01  XMOO  CH 

z*-io  cn  yo*o  cn  tehm70 

22 

r  cn  z  cn  <u>  cn/« 

I.OOOE-OI  -l.OOOE+OI  4.410E+00 
1.200E-01  -1 .OOOE+OI  4. 710E+00 
1.400E-01  -l.OOOE+OI  5.070E+00 
1. 700E-01  -l.OOOE+OI  5.900E+00 
2.000E-01  -l.OOOE+OI  6.670E+00 
2.400E-01  -l.OOOE+OI  7.470E+00 
3.000E-01  -l.OOOE+OI  O.SOOE+OO 
4.000E-01  -l.OOOE+OI  9.130E+00 
5.200E-01  -l.OOOE+OI  9.S20E+00 
e.SOOE-OI  -l.OOOE+OI  l.OOOE+OI 
0.200E-Q1  -l.OOOE+OI  1.049E+01 
1.040E+00  -l.OOOE+OI  1 . 091E+01 
1 . 300E+00  -l.OOOE+OI  1. 119E+01 
1 . 030E+00  -l.OOOE+OI  1.109E+0I 
2.030E+00  -l.OOOE+OI  1 . 194E+01 
2. SOOE+OO  -l.OOOE+OI  1.2S7E+01 
3. 190E+00  -l.OOOE+OI  1.290E+01 
3. 700E+00  -l.OOOE+OI  1.342E+01 
4. 310E+00  -l.OOOE+OI  1.3901+01 
S. OOOE+OO  -l.OOOE+OI  1.427E+01 
6.000E+00  -l.OOOE+OI  1.404E+01 
0. OOOE+OO  -l.OOOE+OI  1.490C+01 

QUAL7  5/10/01  XMOO  CH 

z«-20  cn  to*-. 03  on  tehmzo 
22 

t  cn  z  cn  <»>  en/t 

7.000E-02  -2.000E+01  3.I79E+00 

9.000E-02  -2.000E+01  4.150C+00 

1 . 100E-01  -2.000E+01  4.9201+00 

1.400E-01  -2.000E+01  0. 190E+00 

1.700C-01  -2.000E+01  0.090E+00 

2. 100E-01  -2.009C+01  7.430E+00 

2. 700E-01  -2.000E+01  0.3201+99 


3.700E-01 

-2.000E+01 

9. 260E+00 

4. 900E-01 

-2.000E+01 

9. 050E+Q0 

0.200E-01 

-2.000E+01 

1.021E+01 

7. 900E-01 

-2.000E+01 

1 .066E+01 

1.010E+00 

-2.000E+01 

1.119E+01 

1 . 270E+00 

•2.000E+01 

1 . 140E+01 

1. OOOE+OO 

-2.000E+0 1 

1. 109E+01 

2. OOOE+OO 

-2.000E+01 

1 . 220E+01 

2.470E+00 

-2.000E+01 

1.261E+01 

3. 100E+00 

-2.000E+01 

1 .325E+01 

3.070E+00 

-2.000E+01 

1 . 370E+01 

4.200E+00 

-2.000E+01 

1.423E+01 

4.970E+00 

-2.000E+01 

1.4S0E+01 

S.970E+00 

-2.8Q0E+01 

1 . 474E+0 1 

7. 970E+00 

-2.000E+01 

1.490E+01 

FL400*FlAT5+n.AT3  5/27-20/01 

xmoo  cn 

H*27. o  cn 

TEHOMI 

54 

y  cn 

z  cn 

<U>  CPU'S 

<U'«U'> 

I.OOOE-OI 

O.OOOE-OI 

4.070E+00 

2 

S32E+00 

1.200E-01 

O.OOOE-OI 

4.920E+00 

2 

554E+00 

1.400E-01 

O.OOOE-OI 

5. SOOE+OO 

3 

100E+00 

1.700E-01 

O.OOOE-OI 

S.910E+00 

3 

113E+00 

2.000E-01 

O.OOOE-OI 

0. OOOE+OO 

3 

475E+00 

2.400E-01 

O.OOOE-OI 

7. 100E+00 

3 

15SE+00 

3.000E-OI 

O.OOOE-OI 

O.OtOE+OO 

3 

0S9E+00 

4.000E-01 

0.090E-01 

0.960E+00 

2 

S51E+0Q 

5. 200E-01 

O.OOOE-OI 

9.490E+00 

1 

922E+00 

6. 500E-01 

O.OOOE-OI 

• . 720E+00 

1 

502E+00 

0.200E-91 

O.OOOE-OI 

1.02SE+01 

1 

200E+00 

1.040E+00 

O.OOOE-OI 

1.050C+01 

1 

I20E+00 

1 . 300E+00 

O.OOOE-OI 

1.090E+01 

1 

055E+00 

1.030E+00 

O.OOOE-OI 

1. 120E+01 

9 

•70E-01 

2.030E+00 

O.OOOE-OI 

1.173E+01 

1 

OOOE+OO 

2. SOOE+OO 

O.OOOE-OI 

1.200E+01 

9 

390E-01 

3. 190E+00 

O.OOOE-OI 

1.2C5E+01 

9 

240E-01 

3.700E+00 

O.OOOE-OI 

1 . 301E+01 

• 

430E-01 

4.310E+00 

O.OOOE-OI 

1 . 350E+01 

7 

700E-01 

5.  OOOE+OO 

O.OOOE-OI 

1 . 401E+01 

5 

920E-01 

•.OOOE+OO 

O.OOOE-OI 

1.455E+01 

2 

700E-01 

t.OOIE+OO 

O.OOOE-OI 

1.430E+01 

4 

tOOE-Ot 

0. OOOE+OO 

O.OOOE-OI 

1.493E+01 

3 

700E-02 

0.901E+00 

O.OOOE-OI 

1.477E+01 

4 

•OOE-02 

l.OOOE+OI 

O.OOOE-OI 

1.40SE+01 

1 

000E-02 

1.001C+01 

o.oeoE-oi 

1.470E+01 

2 

1 00E-02 

1.200E+01 

O.OOOE-OI 

1.473E+01 

1 

200C-02 

1 . 201E+01 

0.800E-01 

1.472E+01 

1 

000E-02 

1.400E+01 

O.OOOE-OI 

1.477E+01 

1 

OOOE-02 

1 . 401E+01 

O.OOOE-OI 

1.470E+01 

2 

000E-02 

1.090E+01 

O.OOOE-OI 

1.471E+01 

2 

9002-02 

1.001E+01 

O.OOOE-OI 

1.475E+01 

2 

50SE-02 

l.OOOE+OI 

O.OOOE-OI 

1 . 400E+01 

1 

1S0E-01 

2.0001+01 

O.OOOE-OI 

1 . 377E+0 1 

5 

070E-01 

2.1002*01 

0.0002-01 

1.3292*01 

7.1102-01 

2. 1692*01 

0.0002-01 

1.2012*01 

0.1302-01 

2. 2301*01 

0.0002-01 

1.2462*01 

9.0702-0! 

2.2012*01 

0.0002-01 

1.2162*01 

0.7702-01 

2.3502*01 

0.0002-01 

1.1752*01 

9.2002-01 

2.3972*01 

0.0002-01 

1.1252*01 

9.0102-01 

2.4372*01 

0.0002-01 

1.0922*01 

9.3902-01 

2.4702*01 

0.0002-01 

1.0002*01 

1.0152*00 

2.4962*01 

0.0002-01 

1.0422*01 

1.0262*00 

2.5102*01 

0.0002-01 

1.0172*01 

1.2642*00 

2.5352*01 

0.0002-01 

9.9002*00 

1 . 3792*00 

2.5402*01 

0.0002-01 

9.6002*00 

1.6462*00 

2.5602*01 

0.0002-01 

9.2602*00 

1.9512*00 

2.5702*01 

0.0002-01 

0.5002*00 

2.6092*00 

2.5762*01 

0.0002-01 

0.1902*00 

2.7732*00 

2.5002*01 

0.0002-01 

0.1302*00 

2.7372*00 

2.5032*01 

0.0002-01 

7.3102*00 

3.2762*00 

2.5062*01 

0.0002-01 

6.0002*00 

3.2052*00 

2.5002*01 

0.0002-01 

6.5702*00 

3. 1342*00 

2.5902*01 

0.0002-01 

6.3702*00 

3.3022*00 

atv30«cimvi*cuRV4  5/01 

TH2TA*30 

49 

r  en 

Mzts.Q  en 

T2HM01 

2  cn 

<u>  en/s 

<«'•»'> 

1.2002-01 

0.0002-01 

3.9302*00 

1.7122*00 

1.4002-01 

0.0002-01 

4.6302*00 

2.3142*00 

1.7002-01 

0.0002-01 

5.0902*00 

3.2012*00 

2.0002-01 

0.0002-01 

7.0002*00 

3.9652*00 

2.4002-01 

0.0002-01 

0.1302*00 

3.5962*00 

3.0002-01 

0.0002-01 

0.6602*00 

3.6342*00 

4.0002-0 1 

0.0002-0 1 

9.7102*00 

3.3702*00 

5.2002-01 

0.0002-01 

1.0322*01 

2.3002*00 

6.5002-01 

0.0002-01 

1.0502*01 

1.0772*00 

0.2002-01 

0.0002-01 

1.0702*01 

1.9252*00 

1.0402*00 

0.0002-01 

1.1062*01 

1.7502*00 

1.3002*00 

0.0002-01 

1.1202*01 

1 . 5352*00 

1.6302*00 

0.0002-01 

1.1472*01 

1.4592*00 

2.0302*00 

0.0002-01 

1.1012*01 

1.4702*00 

2.5002*00 

0.0002-01 

1.1902*01 

1.4702*00 

3. 1902*00 

0.0002-01 

1.2472*01 

1.4422*00 

3.7002*00 

0.0002-01 

1.2602*01 

1 . 5532*00 

4.3102*00 

0.0002-01 

1.3122*01 

1.3012*00 

5.0002*00 

0.0002-01 

1.3532*01 

1.1512*00 

5.0012*00 

0.0002-01 

1.3542*01 

1.1202*00 

6.0002*00 

0.0002-01 

1.4112*01 

7.0302-01 

7.0002*00 

0.0002-01 

1.4612*01 

1.3702-01 

0.0002*00 

0.0002-01 

1.4912*01 

2.3502-0! 

9.0002*00 

0.0002-01 

1.5102*01 

1.3602-01 

1.1002*01 

0.0002-01 

1.5532*01 

3.0002-02 

1.3002*01 

0.0002-01 

1.5602*01 

3.2002-02 

1.5002*01 

0.0002-01 

1.6022*01 

2.6002-02 

1.7002*01 

0.0002-01 

1.6022*01 

7.0002-02 

1.9002*01 

0.0002-01 

1.3732*0! 

2.2402-01 

-  Ill 


2.0002*01 

0.0002-01 

1.5592*01 

2.9002-01 

2.0692*01 

0.0002-01 

1.5172*01 

2.5202-01 

2.1302*01 

0.0002-01 

1.5132*01 

2.6702-01 

2.1012*01 

0.0002-01 

1.4002*01 

2.7702-01 

2.2502*01 

0.0002-0t 

1.4672*01 

3.0102-01 

2.2972*01 

0.0002-01 

1.4302*01 

3.5602-01 

2.3372*01 

0.0002-01 

1.4102*01 

4.2402-01 

2.3702*01 

0.0002-01 

1.3742*01 

6.5302-01 

2.3962*01 

0.0002-01 

1.3122*01 

1.0302*00 

2.4102*01 

0.0002-01 

1.2642*01 

1.4662*00 

2.4352*01 

0.0002-01 

1.1072*01 

1.0312*00 

2.4402*01 

0.0002-01 

1 . 1 10E*01 

2.1112*00 

2.4602*01 

0.0002-01 

9.9202*00 

2.7612*00 

2.4702*01 

0.0002-01 

0.7502*00 

3.0032*00 

2.4762*01 

0.0002-01 

7.6202*00 

3.4522*00 

2.4002*01 

0.0002-01 

6.7502*00 

3.6172*00 

2.4032*01 

0.0002-01 

6.3502*00 

3.4522*00 

2.4062*01 

0.0002-01 

5.6902*00 

3.3922*00 

2.4002*01 

0.0002-01 

4.5702*00 

2.6092*00 

2.4902*01 

0.0002-01 

3.0102*00 

2.2112*00 

CRV7S*eURV3*CURV5  5/01 

TN2TA*75 

49 

y  en 

m>25.7  en 

T2HP*01 

2  en 

<u>  cn/s 

1.2002-01 

0.0002-01 

6.5602*00 

4.2502*00 

1.4002-01 

0.0002-01 

7.6202*00 

4.2192*00 

1.7002-01 

0.0102*00 

4.0442*00 

2.0002-01 

0.0002-01 

0.6002*00 

4.4572*00 

2.4002-01 

0. 0002-0 1 

9.3102*00 

4.6932*00 

3.0002-01 

0.0002-01 

1.0202*01 

4.2652*00 

4.0002-01 

0.0002-01 

1.0062*01 

4.0012*00 

5.2002-01 

0.0002-01 

1.1352*01 

3.2562*00 

6.5002-01 

0.0002-01 

1.1712*01 

2.0742*00 

0.2002-01 

0.0002-01 

1.2112*01 

2.4062*00 

1.0402*00 

0.0002-01 

1.2452*01 

2.0232*00 

1.3002*00 

0.0002-01 

1.2992*01 

1.6602*00 

1.6302*00 

0.0002-01 

1.2072*01 

1.4142*00 

2.0302*00 

0.0002-01 

1.3022*01 

1.3692*00 

2.5002*00 

0.0002-01 

1.3252*01 

1.3532*00 

3.1902*00 

0.0002-01 

1.3542*01 

1.1302*00 

3.7002*00 

0.0002-01 

1.3602*01 

1.2362*00 

4.3102*00 

0.0002-01 

1.3572*01 

1.1072*00 

5.0002*00 

0.0002-01 

1.3932*01 

9.7502-01 

6.0002*00 

0  ■  Q0QE~Q 1 

1.4102*01 

9.9202-01 

0.0002*00 

Q  QQQf—Q 1 

1.4742*01 

7.2602-01 

1.0002*01 

Q  QOQt~0 1 

1.5362*01 

4.0702-01 

1.2002*01 

Q  Q0QE‘~u  1 

1.5032*01 

3.6402-01 

1.4002*01 

0.0002-01 

1.6092*01 

1.4602-01 

1.5702*01 

0.0002-01 

1.6572*01 

7.0002-02 

1.6002*01 

0.0002-01 

1.6502*01 

6.4002-02 

1.7702*01 

Q  QQQe  ~0 1 

1.6302*01 

1.0402-01 

1.0002*01 

Q  QQQf  «■  0  1 

1.6422*01 

1.1102-01 

1.9702*01 

a  0001*0 1 

1.5602*01 

1.0002-01 

2.070E+01 

O.OOOE-OI 

1 . 524E+01 

2.620E-01 

2.200E401 

O.OOOE-OI 

1.066E401 

2.064E400 

2. 139E+01 

O.OOOE-OI 

1.469E+01 

3.030E-01 

2. 269E401 

O.OOOE-OI 

1 . 039E+01 

2.311E400 

2.200E401 

O.OOOE-OI 

1.448E+01 

3.760E-01 

2.330E4Q1 

O.OOOE-OI 

9.700E400 

2.295E400 

2.251E+01 

O.OOOE-OI 

1 . 409E+01 

4. 540E-01 

2.301E401 

O.OOOE-Ot 

9. 250E+00 

2.272E+00 

2.320E+01 

O.OOOE-OI 

1 . 341E+01 

S.S90E-01 

2.4S0E401 

O.OOOE-OI 

6.540E400 

2.145E400 

2.367001 

O.OOOE-OI 

1.271E+01 

6.610E-01 

2.497E401 

0 . 000E-01 

7.950E400 

1.9S1E400 

2.407E+01 

O.OOOE-OI 

1 . 191E+01 

7.270E-01 

2.537E401 

O.OOOE-OI 

7.690E400 

1 . 793E+00 

2.440E401 

O.OOOE-OI 

1. 123E+01 

6. 150E-01 

2.570E401 

0 . 000E-01 

7.470E400 

1.697E400 

2.466001 

O.OOOE-OI 

1.063E+01 

1.0S9E400 

2.596E401 

0. 000E-01 

7.060E400 

1.624E400 

2.466001 

0. 000E-01 

1.0072*01 

1.200E400 

2.616E401 

O.OOOE-OI 

6.940E400 

1.646E400 

2. SO 5001 

O.OOOE-OI 

9. 390E+00 

1.452E400 

2.635E401 

O.OOOE-OI 

6. 500E+00 

1 . 544E400 

2.516001 

O.OOOE-OI 

6.010E+00 

1.677E+00 

2.646E401 

O.OOOE-OI 

6. 300E+00 

1 . 541E+00 

2.530001 

O.OOOE-OI 

6. 18OE.0O 

2.173E400 

2.660E401 

0. OOQE-OI 

6. 160E400 

1.690E400 

2.S40E401 

O.OOOE-OI 

7.300E+00 

2.436E+00 

2.6702401 

O.OOOE-OI 

5. 330E+00 

1.779E400 

2.546E+01 

O.OOOE-OI 

6.570E+00 

2.7532400 

2.676E401 

0 . 000E-01 

5.110E400 

1 . 906E+00 

2.5S0E+01 

O.OOOE-OI 

5. I70E+00 

2.609E400 

2.600E401 

O.OOOE-OI 

4.600E400 

1.698E400 

2.553E+01 

O.OOOE-OI 

5. S30E+00 

2.395E400 

2.663E401 

O.OOOE-OI 

3.990E400 

1.S52E400 

2. 556E+01 

O.OOOE-OI 

4.710E+00 

1.9932400 

2.666E401 

O.OOOE-OI 

3.440E400 

1.368E400 

2.5S6E+01 

O.OOOE-OI 

4 . 090E+00 

1.906E400 

2.666E401 

0. OOOE-O 1 

2.790E400 

1.188E400 

2. 560E+01 

O.OOOE-OI 

3.260E+00 

1.3622400 

CORO 4  5/20/61  X»400 

cn  TEnP*61 

FL7S0*CURV6+CURY7  5/61 

42 

X»760  Ctl 

H=27.0  cn 

TEHPM1 

y  en 

z  cn 

<u>  cn/s 

46 

2.0002400 

-2.000E401 

1 . 1 79E+0 1 

1.074E400 

r  cn 

z  cn 

<u>  cn/t 

<«'««'> 

4.000E400 

-2.000E401 

1.355E+01 

6.050E-01 

1.200E-01 

4.560C+00 

2.5672400 

2.0002400 

-1.6C0E+01 

1. 162E4Q1 

1.003E400 

1.400E-01 

O.OOOE-OI 

5.680E+09 

3.662E400 

4.000E400 

-1.600E401 

1.348E401 

7.280E-01 

1.700E-01 

O.OOOE-OI 

7.330E+00 

4.0632400 

2.0002400 

-1.600E401 

1.177E491 

9.730E-01 

I.000E-01 

O.OOOE-OI 

7.960E+00 

4.4072400 

4.0002400 

-1.600E401 

1.354E401 

7.490E-01 

2.4Q0E-01 

Q  OOOE-Ol 

9. 1S0E+00 

4.4402400 

2.9002400 

-1.400E401 

1.193E401 

9. 690E-01 

3. 000E-01 

0  OOOE-O 1 

1.03IE+01 

4.2612400 

4.0002400 

-1.400E401 

1.352E401 

6.960E-01 

4.000E-01 

n  OOOE^O 1 

1 . 1682*01 

3.6272400 

2.0002400 

-1.200E+01 

1.187E+01 

9.770E-01 

5.200E-01 

O.OOOE-OI 

1 . 2  52C-.0 1 

3.S10E400 

4.0002400 

-1.2002401 

1.3S9E4Q1 

7.660E-01 

6.500E-01 

0.Q0QC-Q1 

1.332E.01 

3.1792460 

2.0002400 

-1.0002401 

1.178E401 

1.019E400 

6.200E-01 

1 . 347E+01 

2.6612400 

4.0002400 

-1.0002401 

1.3542401 

7.430E-01 

1 . 040E+00 

O.OOOE-OI 

1 .406E+01 

2.4022400 

2.0002400 

-6.0002400 

1.161E401 

9.820E-01 

1 . 300C+00 

1.4482*01 

2.0022400 

4.0002400 

-0.0002400 

1.337E401 

7.600E-01 

1 . 630E+00 

O.OOOE-OI 

1 . 46 IE* 01 

1.7002400 

2.0002400 

-6.0002400 

1.174E401 

9.740E-01 

2.030E+00 

O.OOOE-OI 

1.4952*01 

1.3662400 

4.0002400 

-6.000E400 

1.340E+01 

7.990E-0I 

2.500E+00 

1.479E+01 

1.3432400 

2.0002400 

-4.0002400 

1.167E401 

1.035E400 

3. 190E+00 

0. 00QE-01 

1.5062*01 

1.0372400 

4.0002400 

-4.000E400 

1.341E401 

7. 100E-01 

3.700E+00 

O.OOOE-OI 

1.5172*01 

9.060E-01 

2.0002400 

-2.0002400 

1.161E401 

8.600E-01 

4. 3I0E+00 

o.oooe-oi 

1.5HE+01 

0. 3002-01 

4.0002400 

1.327E+01 

6.5702-01 

5.000E+00 

0. 000E-01 

1 . 527E+01 

9.3602-01 

2.0002400 

O.OOOE-OI 

1.170E401 

9.900E-01 

6.000E+00 

1.524E+01 

9.9302-01 

4.0002400 

0.0002-01 

1.341E401 

7.740E-01 

6.000E400 

O.OOOE-OI 

1.5562*01 

6.3362-01 

2.0002400 

2.0002400 

1.169E401 

9. 130E-01 

1.O0OE4O1 

O.OOOE-OI 

1.U4M1 

6.3102-01 

4.0002400 

2.0002400 

1.3312401 

7.570E-01 

1.200E401 

q  0bSC*0 i 

1.3»tE*01 

4.2202-01 

2.0002400 

4.0002400 

1.1602401 

9.260E-01 

I.300E+01 

q  OOOE— 0 l 

1.56214-01 

5.2102-01 

4.0002400 

4.0002400 

1.3272401 

6.200E-01 

1.400E401 

0. OOOC-OI 

1 . 562001 

2.5102-01 

2.0002400 

6.0002400 

1.164E401 

9.760E-01 

1.5002401 

O.OOOE-OI 

1 . 566C+0 1 

3.0002-01 

4.0002400 

6.0002400 

I.343E401 

7.2502-01 

1 .7002+01 

o.oooe-oi 

1.49724-01 

3.7002-01 

2.0002400 

6.0002400 

1.1602*01 

9.5502-01 

1 . 900E+0 1 

1 . 3532*01 

9.9402-01 

4.0002400 

0.0002400 

1.3262401 

7.9702-01 

2. 1002401 

0. OOOE-O 1 

1.1512401 

1.5912400 

2.0002400 

1.0002401 

1.1632401 

1.0122400 

-  t*«  - 


4. 000E+00 

1.000E+01 

1.324E+0I 

8.300E-01 

1.550E+01  2. OOOE+OO 

1.474E+01 

2. OOOE+OO 

1.200E+01 

1. 161E+01 

9. 500E-01 

I.900E+01  2. OOOE+OO 

1.4S8E+01 

4. OOOE+OO 

1.200E+01 

1.329E+01 

8.110E-01 

8. OOOE+OO  0.000E-01 

1 . 488E+01 

2.000000 

1.400E+01 

1 . 104E+01 

9. 190E-01 

1 . 150E+01  0.Q00E-01 

1.479E+01 

4. OOOE+OO 

1.400E+01 

1.321E+01 

8.030E-01 

1.550E+01  O.OOOE-OI 

1.476E+01 

2. OOOE+OO 

1.000E+01 

1.154E+01 

9.460E-01 

1.900E+01  0. 000E-01 

1.4SSE+01 

4.000000 

1.000E+01 

1.322E+01 

7.580E-01 

8.  OOOE+OO  -2. OOOE+OO 

1.490E+01 

2.000000 

1.800E+01 

1 . 161E+01 

9. 380E-OI 

I.I50E+01  -2. OOOE+OO 

1.480E+01 

4.000000 

1.800E+01 

1.328E+01 

8.780E-01 

1.5S0E+01  -2. OOOE+OO 

1.477E+01 

2.  OOOE+OO 

2.000E+01 

1. 16SE+01 

9.670E-01 

1.900E+Ot  -2. OOOE+OO 

1 . 482E+01 

4.000E+00 

2.000E+01 

1.324E+01 

7.590E-01 

8. OOOE+OO  -4. OOOE+OO 

1.486E+01 

1 . 150E+01  -4. OOOE+OO 

1.477E+01 

CORE  1  5/13/01  X*400 

CM 

1.550E+01  -4. OOOE+OO 

1.474E+01 

T*60  SEC 

TEMP*76 

1.900E+01  -4. OOOE+OO 

1.4S8E+01 

84 

8. OOOE+OO  -6. OOOE+OO 

1.489E+01 

T  CM 

2  CM 

<U>  CM/S 

I.ISOE+Ot  -6. OOOE+OO 

1.481E+01 

8.  OOOE+OO 

2.000E+01 

1.481E+Q1 

I.550E+81  -8. OOOE+OO 

1.474E+01 

1 . 1 50E+0 1 

2.000E+01 

t . 483E+0 1 

1.900E+0t  -8. OOOE+OO 

1.454E+01 

1.S50E+01 

2.000E+01 

1.470E+01 

8. OOOE+OO  -8. OOOE+OO 

1.48SE+01 

1.900E+01 

2.000E+01 

1.406E+01 

1 . 150E+01  -8. OOOE+OO 

1.473E+01 

8. OOOE+OO 

1 . 800E+0 I 

1.489E+01 

1.S50E+01  -8. OOOE+OO 

1.473E+01 

1 . 150E+01 

1.800E+01 

1.479E+01 

1.900E+01  -8. OOOE+OO 

1.450E+01 

1 . 550E+01 

1.800E+01 

1.471E+01 

8.  OOOE+OO  -1.000E+01 

1.492E+01 

1.900E+01 

1.800E+01 

1.459E+01 

1. 1S0E+01  -1.000E+01 

1 . 476E+01 

8.  OOOE+OO 

1.600E+01 

1.4S8E+01 

1.S50E+01  -1.000E+01 

1.472E+01 

1 . 15QE+01 

1.S00E+01 

1.480E+01 

1.900E+01  -1.800E+01 

1.454E+01 

1.550E+01 

1.50QE+01 

1.470E+01 

8. OOOE+OO  -I.200E+01 

1. 487E+01 

1.900E+01 

1 . 600E+0 1 

1.483E+01 

1. 150E+01  -1.200E+01 

1.479E+01 

8.  OOOE+OO 

1 . 400E+01 

1 . 485E+01 

1.550E+01  -1.200E+01 

1.472E+01 

1. 150E+01 

1.400E+01 

1.4S0E+01 

1.900E+01  -1.200E+01 

1.4S6E+01 

1.550E+01 

1.400E+01 

1.469E+01 

8.  OOOE+OO  -1.400E+01 

1 . 488E+0 1 

1.900E+01 

1.400E+01 

1.459E+01 

1 . 150E+01  -1.400E+01 

1.480E+01 

8. OOOE+OO 

1.200E+01 

1.407E+01 

1  .SSQE+01  -1.400E+01 

1.473E+01 

1. 150E+01 

1.200E+01 

1.479E+01 

1.900E+01  -1.400E+01 

1.458E+01 

1.S50E+01 

1.200E+01 

1.472E+01 

8. OOOE+OO  -1.600E+01 

1.493E+01 

1.900E+01 

1.200E+01 

1.400E+0I 

1. 150E+01  - 1 . 800E+0 1 

1 . 478E+0 1 

8. OOOE+OO 

I.000E+01 

1.483E+01 

1.550E+01  -1.500E+01 

1 . 475E+0 1 

1.150E+01 

1.000E+01 

1 . 480E+0 1 

1.900E+01  -1.600E+01 

1.448E+01 

1.550E+01 

1.000E+01 

I.474E+01 

8. OOOE+OO  -1.800E+01 

1.492E+01 

1.900E+01 

1 . 000E+01 

1 . 462E+01 

1 . 15QE+01  -1.800E+01 

1 . 477E+0 1 

8.  OOOE+OO 

8. OOOE+OO 

1.491E+01 

1.550E+01  -1.800E+01 

1.475E+01 

1 . 1S0E+01 

8. OOOE+OO 

1.480E+01 

1.900E+01  -1.800E+01 

1.444E+01 

1.550E+01 

8. OOOE+OO 

1.473E+01 

8. OOOE+OO  -2.000E+01 

1.488E+01 

1.900E+01 

8. OOOE+OO 

1.460E+01 

1. 1SOE+OI  -2.000E+01 

I.479E+01 

8. OOOE+OO 

0. OOOE+OO 

1.487E+01 

1.550E+01  -2. 000E+01 

1.474E+01 

1. I50E+01 

6. OOOE+OO 

1. 481E+01 

1.900E+01  -2.000E+01 

1.4S3E+01 

1.S50E+01 

0. OOOE+OO 

1.472E+01 

1.900E+01 

0.  OOOE+OO 

1.4S9E+01 

CORE*  5/25/01  THETA«30  TEMP* 

8. OOOE+OO 

4. OOOE+OO 

1.483E+01 

105 

1 . 1S0E+01 

4. OOOE+OO 

1.477E+01 

T  CM  *  CM 

<U>  CM/S 

1.550E+01 

4. OOOE+OO 

1.471E+01 

2. OOOE+OO  -2.000E+01 

1. 141E+01 

1.900E+01 

4. OOOE+OO 

1.4S7E+01 

4. OOOE+OO  -2.000E+01 

1.270E+01 

8.  OOOE+OO 

2. OOOE+OO 

1.488E+01 

8. OOOE+OO  -2.000E+01 

1 . 501E+01 

1.150E+01 

2. OOOE+OO 

1.479E+01 

1.200E+01  -2. 000E+01 

1.571E+01 

.52oe+oo 
I . 542E+00 

.oooe-oi 

.5001-0* 
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I . 600E+Q 1  -2.000E+01 
2. OOOE+OO  - 1 . 800E+01 
4. OOOE+OO  -1.80QE+01 
8. OOOE+OO  -1 .800E+01 
1 . 200E+Q1  -1.800E+01 
1.600E+0I  -1.800E+01 
2. OOOE+OO  -1.60QE+01 
4. OOOE+OO  -1.600E+01 
0. OOOE+OO  -1 . 600E+01 
1.200E+01  -1.600E+QI 
1.600E+01  >1 .600001 
2.000000  -1.400001 
4. OOOE+OO  -1.400001 
I. OOOE+OO  -1.400E+01 
1.200001  -1.400001 
1.600001  -1.400001 
2. OOOE+OO  -1 . 200E+01 
4. OOOE+OO  -1.200E+01 
6. OOOE+OO  -1.200E+01 
1.200E+01  — 1 . 200E+0 1 
1.600E+01  -1.200E+01 
2. OOOE+OO  -1.000E+01 
4. OOOE+OO  -1.000E+01 
0. OOOE+OO  -1.000E+01 
1.200E+01  -1 . 000E+01 
1.600E+01  -1.000E+01 
2. OOOE+OO  -6. OOOE+OO 
4. OOOE+OO  -0. OOOE+OO 
6. OOOE+OO  -0. OOOE+OO 
1.200E+01  -0. OOOE+OO 
1.600E+01  -0. OOOE+OO 
2. OOOE+OO  -6. OOOE+OO 
4. OOOE+OO  -6. OOOE+OO 
0. OOOE+OO  -6. OOOE+OO 
1.200E+01  -6. OOOE+OO 
1.600E+01  -6. OOOE+OO 
2. OOOE+OO  -4. OOOE+OO 
4. OOOE+OO  -4. OOOE+OO 
0. OOOE+OO  -4. OOOE+OO 
1.200E+01  -4. OOOE+OO 
1.600E+01  -4. OOOE+OO 
2. OOOE+OO  -2. OOOE+OO 
4. OOOE+OO  -2. OOOE+OO 
0. OOOE+OO  -2. OOOE+OO 
1.200E+01  -2. OOOE+OO 
1.600E+01  -2. OOOE+OO 
2. OOOE+OO  O.OOQE-OI 
4. OOOE+OO  O.OOOE-OI 
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